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A  computer-based  system  has  been  created  using  29  features  of  hand 
geometry  to  discriminate  successfully  278  images  of  hands  taken  from 
30  male  and  female  subjects  over  a  four-month  period.  These  features 
include  hand  perimeter  and  area  measures,  finger  lengths  and  widths, 
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subjects  from  which  the  images  were  taken.  All  images  were  correctly 


DO  1473  EDITION  OF  I  NOV  «S  IS  OBSOLETE 


UNCLASSIFIED 

s E CURlTY  CLASSIFICATION  OF  THIS  PAGE  Dum  SAM* 


yy/y  ^ 


•'■VUJSJitUW**'' 


CONTENTS 


LIST  OF  ILLUSTRATIONS 


ACKNOWLEDGMENT 


I  INTRODUCTION 


II  DATA  COLLECTION 


III  DATA  PREPARATION 


A.  Data  Generation 

B.  Data  Reduction 


IV  DATA  ANALYSIS 


V  CLASSIFICATION  RESULTS 


VI  COMPUTER  PROGRAM  SUMMARIES 


A.  HIST  . 

B.  THRESH 


C.  LINEBYLINE 


D.  SCAN 


APPENDICES 

A— COMPUTER  PROGRAM  LISTINGS 


B— DOCUMENTATION  OF  LINEBYLINE 


C— ANNOTATIONS  FOR  THE  SEPARATE  SPSS  RESULT  PACKAGE 
(Submitted  Separately)  . 


Accession  For 

NTIS  GE/iil 
VTIC  TAB 

Uuannowiood 

Justif i  c?.i  ion — 


By- - 

Distribut ion/ 

Avail.'.i  : Ity  Coles 
i  Av'ii  c.nJ/er 

Diet  Sr-ciol 


ILLUSTRATIONS 


1  Apparatus  Used  for  Data  Collection .  4 

2  Fingertip  Location  Measurement  .  10 

3  Measurement  Points  for  Lengths  of  Finger  and  Thumb  .  10 

4  Finger  and  Thumb  Width  Measurement  Points  ...  .  11 

5  Curvature  and  Distance  Measures  .  12 

6  Analysis  of  Hand  Geometry  Features  .....  .  14 

7  Analysis  of  Hand  Geometry  Features  Using  Within-Groups 

Correlation  Matrix  .  16 


v 


I 


pacsbuo  rtm  num 


ACKNOWLEDGMENT 


With  appreciation,  the  authors  gratefully  acknowledge  the  substan¬ 
tial  contributions  of  the  following  toward  the  research  reported  in 
this  document: 

•  Dr.  John  Ostrem  served  as  a  consultant  on  feature-extraction 
methods,  on  data  analysis,  and — particularly — on  the  use  of 
commercial  software  used  in  the  statistical  analysis  of  the 
feature-extraction  results. 

•  Mr.  Gregory  Myers  supplied  the  computer  program  LINEBYLINE 
and  provided  guidance  in  its  application  to  this  research. 

He  also  gave  valuable  help  and  advice  in  feature  selection 
and  analysis. 

•  Mr.  Zev  Pressman  gave  valuable  guidance  in  the  design  of  the 
data-collection  apparatus  and  in  its  initial  adjustment  and 
evaluation . 

•  Mr.  Marshall  Wilson  recruited  subjects,  photographed  the 
subjects'  hands,  coordinated  the  film  processing,  and  cataloged 
all  of  the  data  in  the  data  base. 


vii  J 

j  PKBOUam  MS  ■AHT-MOt  YUMD 


I  INTRODUCTION 


This  report  describes  a  one-year  research  and  development  effort  to 
determine  the  feasibility  of  using  hand  geometry  and  other  features  of 
the  hand  as  bases  for  the  construction  of  an  automatic  personnel  authen¬ 
tication  system  for  access  control.  The  subsequent  sections  of  the 
report  describe  the  method  and  extent  of  data  collection,  the  method  of 
data  preparation  and  reduction,  the  data  analysis,  the  experimental 
authentication  results,  and  provide  brief  descriptions  of  the  computer 
programs  that  were  written  for  this  work. 

The  appendices  are  photocopies  of  listings  of  the  computer  programs 
used  and  of  critical  parts  of  the  experimental  statistical  analysis  and 
raw  data. 
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II  DATA  COLLECTION 


Previous  access-control  systems  using  hand  geometry  or  fingerprints 
have  used  devices  that  capture  data  from  physical  contact  with  or  proximity 
to  the  hand  surfaces  being  measured.  Such  contact  or  closeness  can  dis¬ 
tort  the  measurements,  particularly  if  the  hand  can  move  or  exert  pressure. 

A  measurement  system  that  can  avoid  distorting  the  features  being  measured 
will  be  more  reliable,  which  in  turn  will  reduce  the  occurrence  of  incor¬ 
rect  authentication  decisions. 

Precise  orientation  of  the  hand  in  space  and  precise  control  of  finger 
and  thumb  spread  are  difficult  to  accomplish  without  some  mechanical  device, 
however.  To  compensate  for  the  variable  position  of  the  hand  from  trial  to 
trial  (because  there  was  no  device  to  hold  it  in  a  fixed  position)  it  was 
decided  to  concentrate  on  adjusting  the  data.  Precise  data  reduction  would 
allow  easier  and  more  natural  interaction  between  the  person  seeking  auth¬ 
entication  and  the  access-control  system. 

The  system  used  to  collect  data  is  schematically  illustrated  in  Fig¬ 
ure  1.  The  subject's  hand  is  positioned  near  the  light  diffuser;  the  open 
palm  faces  the  camera  and  the  fingers  are  comfortably  spread.  The  subject's 
name  and  session  number  are  displayed  in  the  upper  quarter  of  the  imaged 
area  and  thus  become  an  integral  part  of  the  data  base. 

The  data  base  collected  contains  50  subjects,  25  males  and  25  females, 
ranging  in  age  from  the  early  20s  to  65,  in  height  from  approximately  5' 
to  6 ’3",  and  in  weight  over  a  normally  observed  distribution.  The  total 
number  of  sessions  (five  trials  per  session)  was  1020,  approximately  evenly 
divided  between  males  and  females.  Of  these,  997  sessions  produced  4985 
high-quality  images  with  full  gray  scale.  In  the  remaining  sessions,  flash 
problems  of  failure  to  fire  or  synchronization  incurred;  the  resultant 
images  were  quite  useful  for  geometric  measurements  but  contain  no  gray¬ 
scale  features. 
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FIGURE  1  APPARATUS  USED  FOR  DATA  COLLECTION 


Negative  images  have  been  cataloged  and  stored  in  three  volumes  that 
were  designed  and  manufactured  for  such  use.  Photo  development  followed 
standard  procedures  for  the  black-and-white  high-contrast  copy  film 
(Kodak  5069).  The  film  manufacturer  states  that  if  storage  conditions 
are  correct  (cool  and  dry)  the  archival  qualities  of  this  film  are  ex¬ 
cellent — it  will  retain  quality  for  more  than  30  years. 

Image  quality  is  also  excellent.  Correctly  focused  and  exposed  nega 
tives  have  been  examined  with  a  400X  microscope,  and  minutiae  are  clearly 
visible  on  fingertips. 


Ill  DATA  PREPARATION 


A.  Data  Generation 

Two  methods  were  used  to  obtain  digital-image  data  to  be  used  for 
further  processing. 

•  An  SRI-built  Reticon  CCD  array  scanner  (512  elements  per  line) 
interfaced  to  an  LSI-11  was  used  early  in  the  project. 

•  A  commercially  available  Eikonix  digitizer  (1024-element  Reticon 
*  1500  scan  steps)  was  used  to  supply  final  evaluation  data. 

Data  from  the  first  system  were  generally  unsatisfactory  because  of 
a  combination  of  problems.  The  system  interface  included  a  hardware 
thresholder  for  the  analog  data  generated  by  the  Reticon  array.  The 
user  can  vary  the  threshold,  but  it  is  uniform  for  all  samples  coming 
from  the  CCD  array.  Such  an  arrangement  works  well  only  if  the  image  area 
is  uniformly  lighted,  the  individual  CCD  elements  are  reasonably  uniform  in 
sensitivity,  and  the  scanned  images  have  equivalent  niecom  densities. 
Unfortunately,  these  conditions  were  not  satisfied,  and  the  uniform  hard¬ 
ware  threshold  was  found  to  be  quite  unsatisfactory.  It  is  probably 
possible  to  build  a  system  that  utilizes  a  fixed-level  threshold  for  this 
type  of  application;  however,  such  a  device  with  a  sufficiently  uniform 
large  field  of  view  was  not  available,  and  the  unif orm-threshold  technique 
has  therefore  not  been  thoroughly  evaluate^ 

The  Eikonix  digitizer  is  an  excellent  instrument  for  purposes  of 
this  study.  The  digitized  images  generated  are  1024  x  1500-element  data 
each  with  8  bits  of  gray  scale  recorded  on  magnetic  tape.  The  real  reso¬ 
lution  of  the  data  when  the  film  negatives  are  sampled  with  this  spatial 
frequency  is  two  pixels  (picture  elements)  per  mm.  The  system  can 
automatically  compensate  for  lighting  nonuniformities  and  variable  CCD- 
element  sensitivities.  It  is  also  fairly  fast:  one  complete  cycle  of 
generating  and  recording  an  image  requires  about  one  minute.  All  of  the 
data  used  to  generate  the  results  appearing  in  this  report  were  collected 
through  the  Eikonix  system. 
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B.  Data  Reduction 

Image  data  on  magnetic  tapes  were  transported  to  a  PDP  11/40  computer 
system  in  SRI's  Bioengineering  Research  Center.  Each  image  was  reduced 
by  extracting  the  pertinent  512  *  512-element  field  and  by  thresholding 
to  a  binary  image.  Finally  the  binary  image  was  reduced  to  a  hand- 
perimeter  list  of  X-Y  coordinates  that  was  stored  on  a  DEC  RK  05  disk. 
Although  this  process  is  somewhat  involved,  it  does  provide  final  data 
that  are  quite  reliable  and  easy  to  work  with  in  subsequent  computations. 
(It  is  to  be  noted  that  the  final  data  set  generated  by  this  procedure 
does  not  include  any  information  about  a  subject's  hand  other  than  the 
perimeter  coordinates  in  a  512  x  512-element  field.) 

A  very  important  feature  of  the  data-reduction  process  is  the  method 
used  to  compute  the  perimeter  list  from  a  512  x  512-element  hand  image. 

The  method  involves  using  a  computer  program  (LINEBYLINE)  that  was  de¬ 
veloped  earlier  at  SRI  and  recently  rewritten  in  PASCAL  for  use  on  the 
PDP  11/40.  The  program  generates  the  perimeter  list  while  accessing  the 
image  data  one  line  at  a  time,  starting  at  one  edge  of  the  field.  This 
implies,  for  a  practical  system,  that  a  minimum  of  memory  and  data  buffer¬ 
ing  is  required;  consequently,  ‘there  is  a  significant  saving  of  money  and 
execution  time. 


IV  DATA  ANALYSIS 


The  basis  for  the  analysis  and  the  geometric  features  that  have 
been  extracted  is,  in  each  case,  a  list  of  the  X-Y  coordinates  of  the 
hand  perimeter.  Such  a  perimeter  is  described  with  a  real  resolution 
of  two  pixels  per  mm. 

All  measurements  are  based  on  landmarks  that  are  easily  identified 
on  the  perimeter;  the  tips  of  the  fingers  and  thumb  were  selected  for 
this  study.  These  landmarks  were  identified  by  calculating  a  coarse 
curvature  function  over  the  perimeter  and  finding  criterion  shifts  in 
the  angle  of  a  tangent  to  the  perimeter.  Tangents  were  estimated  by 
passing  an  18-point  window  over  the  perimeter  and  using  the  window  end 
points  to  define  a  vector  in  the  X-Y  image  plane.  The  window  was  moved 
in  nine-point  increments.  This  analysis  gives  a  rough  estimate  of  the 
tip  location. 

The  rough  estimate  is  refined  by  locating  four  additional  points 
on  the  finger,  two  found  25  mm  away  along  the  perimeter  on  each  side 
of  the  estimated  tip  and  two  more  found  35  mm  away  on  each  side.  The 
virtual  lines  connecting  each  pair  of  points  are  bisected,  and  the  inter¬ 
section  of  the  tip  and  the  line  defined  by  the  bisection  points  is  de¬ 
fined  to  be  the  actual  tip  location.  This  process  is  illustrated  in 
Figure  2.  After  a  precise  location  for  each  tip  is  established,  the 
following  features  can  be  extracted: 

•  Features  1  through  5  are  lengths  of  the  fingers  and  thumb, 
starting  with  the  little  finger.  Lengths  are  measured  from 
the  tip  to  the  depth  of  the  interdigital  space  as  illustrated 
in  Figure  3. 

•  Features  6  through  10  are  widths  of  the  fingers  and  thumb. 

A  width  is  the  distance  between  two  points  on  the  perimeter 
of  a  finger;  the  points  are  45  mm  down  the  perimeter  on  each 

f  side  of  the  tips  of  the  little  finger  and  the  thumb,  and  60  mm 

’  on  each  side  of  the  other  t  ^ertips.  Width  measurement 

points  are  illustrated  in  Figure  4.  These  features  are 
ordered  as  were  the  first  five. 


FIGURE  4  FINGER  AND  THUMB  WIDTH  MEASUREMENT 
POINTS 


Features  11  through  15  are  ratios  of  the  measurements  of  finger 
lengths  to  widths. 

Feature  16  is  the  hand  width  measured  between  two  points  on  the 
palm.  One  point  is  9  cm  on  the  perimeter  from  the  tip  of  the 
little  finger  toward  the  wrist  and  the  other  is  11  cm  on  the 
perimeter  from  the  tip  of  the  index  finger  toward  the  thumb. 
These  points  are  labeled  2  and  11,  respectively,  in  Figure  4. 

Features  17,  18,  and  19  are  the  hand  perimeter,  hand  area,  and 
the  ratio  of  the  squared-perimeter  length  to  the  area.  The 
perimeter  and  area  are  measured  from  point  1  to  point  14  in 
Figure  4.  Point  1  is  13  cm  on  the  perimeter  from  the  tip  of 
the  little  finger,  and  point  14  is  10  cm  on  the  perimeter  from 
the  tip  of  the  thumb. 

Ten  features,  20  through  29,  are  paired  values  corresponding 
to  the  curvature  of  each  fingertip  and  the  distance  over  which 
the  tip  shape  can  be  represented  by  a  circular  arc  with  less 
than  a  criterion  value  of  rms  (root  mean  square)  error.  The 
process  of  extracting  a  typical  feature  pair  is  illustrated  in 
Figure  5.  A  tangent  to  the  perimeter,  such  as  T,  has  an  angle, 
a,  as  shown.  If  a  is  plotted  as  a  function  of  distance  along 
the  perimeter  defining  a  finger,  then  a  plot  such  as  the  one 
shown  in  Figure  5  results.  The  slope  of  the  plot  between 
and  P2  is  related  to  the  curvature  of  the  tip.  The  slope  is 
recorded  as  the  curvature.  A  correlation  calculation  over  the 
region  containing  this  characteristic  region  yields  a  high  value 


PERIMETER 


FIGURE  5  CURVATURE  AND  DISTANCE  MEASURES 


until  the  boundary  of  the  region  extends  into  the  flat  segments 
on  either  side  of  the  sloping  section.  Pj  and  P2  are  the  boundary 
points  when  the  correlation  coefficient  (using  a  straight-line 
fit  to  the  sloping  region)  falls  to  0.98.  The  distance  between 
P^  and  P2  is  recorded  as  the  perimeter  length  over  which  the 
fingertip  can  be  approximated  well  with  a  circular  arc. 

In  summary,  hand  geometry  is  expressed  as  a  set  of  29  features  that 
are  measured  by  locating  precisely  five  landmarks  on  the  perimeter  of  the 
hand — the  tips  of  the  fingers  and  the  thumb.  The  feature  set  contains 
16  measures  of  lengths  and  widths  and  ratios  of  lengths  and  widths,  and 
three  measures  of  area  or  area  and  perimeter  length.  Ten  measures  of 
shape  (curvature)  are  included. 
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V  CLASSIFICATION  RESULTS 


Thirty  subjects  (15  male,  15  female)  were  selected  randomly  from 
the  subject  population  and  278  images  taken  from  the  selected  subjects 
were  digitized.  The  data  were  converted  to  perimeter  lists  and  sets  of 
29  features  were  extracted.  The  data  were  submitted  to  a  standard  sta¬ 
tistical  analysis  program*  to  determine  the  power  of  the  extracted 
features  in  discriminating  among  the  30  subjects.  Linear  discriminant 
analysis  was  performed  and  all  278  feature  sets  were  correctly  classi¬ 
fied  for  all  subjects.  The  detailed  results  of  this  analysis,  as  well 
as  the  feature  data  set,  are  submitted  as  a  separate  package.  Annotation 
for  this  package  is  included  as  Appendix  C. 

Of  critical  importance  are  indicators  of  value  and  reliability  of 
the  various  features  selected.  One  such  indicator  is  the  univariate 
F-ratio,  a  number  that  expresses  the  ratio  of  between-groups  variation 
of  a  given  measure  and  the  within-groups  variation  of  that  same  measure. 
Large  values  of  F  indicate  that  a  measure  will  be  relatively  useful  in 
separating  groups  (or  samples  from  different  subjects)  and  small  values 
predict  that  a  measure  will  do  poorly.  Figure  6  summarizes  the  F-ratios 
of  the  selected  29  features  over  the  30-subject,  278-sample  data  base. 

All  of  the  features  appear  to  be  useful  and  reliable,  especially 
the  more  global  measures  of  hand  width,  perimeter,  and  area.  Some  mea¬ 
sures,  notably  the  curvature  and  distance  measures  associated  with  the 
thumb  are  marginal  as  presently  calculated.  They  could  be  made  more 
useful  by  improving  the  algorithm  that  deals  with  the  curvature  of  the 
perimeter  in  the  vicinity  of  the  thumb  and  by  locating  the  thumb  tip 
more  precisely. 
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55.9899 

FEATUR26 
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9.3*99 

FEATUR27 
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33.1988 
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.7087 

3.51*2 

FEATUR29 
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7.3202 

FIGURE  6  ANALYSIS  OF  HAND  GEOMETRY  FEATURES 


In  general,  it  appears  that  the  measures  used  are  being  generated 
very  reliably,  which  results  in  good  (high)  F-ratios.  One  difficulty 
in  previous  systems  that  used  hand  geometry  for  authentication  was  un¬ 
reliable  measures  of  certain  features,  with  correspondingly  small 
F-ratios  and  poor  discriminating  power. 

Figure  7  gives  more  information  about  the  quality  of  the  feature 
set.  In  this  table,  the  within-groups  correlation  matrix  is  listed. 

The  correlation  matrix  is  used  to  express  the  degree  to  which  pairs  of 
features  are  correlated  or  dependent.  Circled  are  the  two  largest  mag¬ 
nitudes  of  correlation  values,  corresponding  to  the  pairs 

•  Second  finger  length  and  the  ratio  of  its  length  and  width 

•  Thumb  width  and  the  ratio  of  its  length  and  width. 

It  is  to  be  expected  that  ratios  of  features  and  the  features  themselves 
will  tend  to  have  relatively  large  correlation  magnitudes,  and  the  data 
in  Figure  7  substantiate  this  expectation.  However,  the  SPSS  results 
show  that  all  features  contribute  positively  in  the  discriminant  analysis. 

In  summary,  we  have  been  able  to  improve  on  previous  attempts  to  use 
features  of  the  hand  for  access  control  in  the  areas  of  accuracy  and  user 
comfort.  This  has  primarily  been  the  result  of  removing  the  requirements 
on  both  the  subject  and  the  mechanical  device  that  produced  either  a 
distortion  of  the  data  or  discomfort  for  the  individual.  Instead,  we 
have  required  the  software  to  do  all  the  registration  and  measurements. 

The  result  is  a  fast  and  reliable  method  that  demands  only  a  minimum  of 
effort  from  the  individual  requesting  access. 

If  this  technique  were  to  be  packaged  into  a  commercial  product, 
we  estimate  that  the  throughput  rate  could  be  less  than  two  seconds. 

The  rate-limiting  step  in  the  process  is  data  acquisition  from  the 
scanner  through  an  interface  and  into  a  computer.  The  throughput  rate 
of  a  practical  interface  could  be  in  excess  of  275K  bytes  per  second. 
Because  each  image  contains  256K  bytes,  the  transfer  would  take  about 
one  second.  This  data  rate  is  well  within  the  state  of  the  art  of 
image  production  using  a  scanning  linear-array  system.  Because  the  data 
can  be  analyzed  a  line  at  a  time,  and  therefore  analyzed  as  it  is 


TmruaeJ 

2IA1UA0* 

MAT0401 

M  ATUA0* 

7UVUA01 

1 .00004 

nn«ai - 

.u*ii 

1.00000 

MAIUAO) 

.1*141 

.11191 

1.00000 

MUMM 

.22704 

•19494 

.7*799 

1.000 JO 

MMUIOl 

,i*ni 

-.014)4 

.14*07 

•20001 

tf*  rutoo 

•omi 

.0147* 

•11940 

.094)4 

MATUA07 

—04)** 

•02944 

•11904 

.100)9 

itnwii 

.inoi 

.07944 

09400 

•11144 

.00447 

.04411 

.11000 

7141 UA10 

.00911 

•00119 

-.0171* 

-.09*70 

Toman 

'  .TtO>* 

•0*042 

.19**4 

7I4TUA1I 

.104*4 

-.00*24 

.•*000 

MAfUAIS 

•0414) 

Trim 

.♦7799 

.41004 

Tmwn - 

•10701 

.«»*! 

'  .907*2 

7IATUA11 

-.0)09) 

-.020*4 

•1*141 

.17010 

MAVUAU 

-.0*747 

.000*0 

•14014 

•29904 

TUIU.IT  ■ 

•170*4 

■««.<. 

"  .Util 

”  *01107 

HiTWlI 

.14972 

•12)44 

.90144 

•40714 

7IA7UA14 

.171)0 

.27779 

-.01740 

•11*19 

TTTlUifl - 

— TTJItt  - 

^.sm* - 

-.04401 

-.01241“’ 

fiatuaii 

.024*0 

-.10*44 

-.04907 

-.22410 

7IAIUA** 

.044)4 

•42194 

.14110 

•24749 

Mtiiail 

-.04*19 

-.09701 

-.10000 

-.24014 

MttUA!« 

-.017** 

•077*4 

•09470 

••T)*l 

MATUAtt 

.014*4 

•1044* 

.0011* 

•04940 

TTlTUfl* - 791*41 

"Tuni - 

•040*4 

-.10742 

MATUAIT 

-.024*0 

.04242 

•04491 

.•2120 

MATUAIS 

-.0244) 

-.0044) 

-.09411 

-.10001 

IhTMII - 

—  12997 

714911 

.01*01 

.00741 

2IATUA10 

2IATUAU 

72ATUA10 

1.00004 

AIAIUAII 

-.044)2 

1.00040 

7t*TUA12 - —  07**4 - 

-714441 - 

1.00040 

MAfUAll 

-.191*0 

.14*11 

.2*077 

1.00000 

AIATUA1A 

•Jl&ll 

•1417* 

•14901 

•99144 

“7TATUAT5 - 

.07124 - 

.0*400 

.19179  “ 

7<47U*1* 

••17494 

-.00990 

-.00070 

AKAIUAIT 

.14)74 

.20421 

.2444* 

.21070 

IIATUAI* - .14)21  - -.1444) - 

-.17141 

-.19400 

MAfUAlt 

-.2*414 

.14711 

.474*1 

.94401 

MATUAfO 

-.01429 

.209*2 

.00401 

-.11047 

mTwn - 

— vwwi - 

714412 - 

-.00*97 

-.14177 

7IATUA** 

-.0**4* 

.09240 

-.00919 

•02011 

flAfUAO 

.14979 

-.29411 

-.19400 

-.20010 

TIATLM24 - 

’*.11441 - 

-.14147  - 

-.019)1 

-.21440 

FtAfUA*) 

-.01149 

-.0*410 

.0)900 

-.0114* 

ff«ruA24 

.074)7 

-.0*749 

-.0*04* 

-.1)240 

H»TWIT 

-.09144  * 

—.009*0 

.07041 

-.01004 

MATUA24 

.27470 

-.1491* 

-.09491 

-.11024 

MATUAI4  —07041 

-.10201 

•00944 

•01477 

2IA7U414 

MA10A20 

M0TUA21 

M4VUA22 

21 A IUA 14 

1.00000 

FI4TUA20 

.02411 

1.00000 

_flATUA21 

_-• 10474 _ 

.92000 

1.40000 

Tuiuifl 

•00004 

.2*41* 

-.14714 

1.00000 

2CATUA21 

-.14490 

-.21010 

•41*44 

-.941*2 

FfATUA** 

-.2*441 

•02100 

•14)90 

-.00049 

2IAIUA24 

-.1*417 

.074*4 

•19)91 

•11070 

2f Af MA24 

-.19419 

.0*470 

•99012 

••24*11 

21 ATUAIT 

-.01144  _ 

•10129 

•20094 

•11100 

TUtu.il 

•.HiH 

-.21200 

.01144 

-•11000 

.0*171 

-.0497* 

-••0447 

-•00901 

Tf.tiMM 

tlATOatf 

2IIVHAII 

n<i«ui 

.mil 

1.00000 

7IAYUA0* 

M41U400 

1:  *  Tutor 

H.toH. 

Ii»Tu*»* 

1.00000 

_ _  _  _ 

•02099 

l. 00000 

•00** » 

•94400 

1.00000 

•uni 

•17019 

•90920 

n m 

•1)279 

•90100 

•97140 

•42949 

1.00000 

-.19714 

•14794 

•19172 

•14944 

•21914 

.09*04 

-.94401 

-.90914 

-.iiiir 

-.0409) 

-.10491 

-.91092 

-•*4904 

-.14114 

•01412 

-•17972 

-.20049 

-.92494 

-.9*000 

.0414* 

-.21499 

-.l*»TT - 

•nmi 

■  ••Tim 

•40202 

-.14971 

-.00271 

-.11774 

-.11071 

.29900 

•29414 

•29149 

•22)44 

•2701? 

•10709 

•04044  ' 

.lllll 

•l.ITT 

•  1*247 

•44479 

•9049) 

.44*04 

•99)20 

-.00000 

••42102 

-.49194 

-.47001 

-.4*040 

•00192 

-•19900 

-.1**01 

i llAO* 

■  •.mil 

•• 10401 

-.00*10 

•  40071 

•19190 

•117)7 

•14914 

-•02*0* 

•14100 

.00)0) 

•01)94 

-.22177 

•19919 

■  IIMl 

.1742# 

'  •2494C 

.07141 

•24004 

•  14100 

•20097 

•27)91 

.04497 

•11901 

•10*01 

.00091 

•24440 

-.11*10 

'  *12914 

•#9994  " 

.  liter* 

.iiib* 

•0424* 

-•09447 

•90014 

•49912 

•moo 

-.22107 

•14249 

•90997 

•94099 

•10940 

•07490 

•01072 

-.94497 

- =79170 - 

'  .99717 

MAIU014 

7UVU019 

M  ATI*  19 

MATUA17 

9IAVU010 

1.00099 

1.49944 

.10112 

-.94042 

-.21910 

•24917 

-.44219 

.27971 

-.12174 

1.99999 

.12919 

”  .40744 - 

1.99999 
- .f424*~ 

- nooooi— 

•4IA42 

.19)49 

-•49402 

•97492 

••9170* 

.4*0*1 

-.21409 

.17994 

.90949 
-.99999  ~ 
•19294 

—  44494 

•42947 

.417*4 

• 92999 - 

•90999 

—  11199"' 
•14194 

•1047* 

-.*4747 

-•22019 

•14909 

-•99104 

•19209 

-.10914 

••41190 

“  .17047 - 

- .11149- 

- '.11149— 

-.11*40 

•99209 

.21*77 

.94199 

.1401* 

-.*1701 

-.11027 

.10999 

•9)041 

•  1*421 

-.90704 

•42191 

-  .14114  - 

- .11171" 

- nwr 

-.14999 

-.94444 

-.94490 

-.99000 

.00100 

-.411*1 

•90490 

.49011 

-.110)2 

_ -.1*0*0 

MAf  UB29 

7IATU0I4 

MAfU929 

M*fUB20 

9KAVUAI? 

1.09440 

.1)141 

•00094 

1.99990 

•44204 

1.99999 

- - 

— 

•97441 

•99719 

•19474 

1.99999 

-.49404 

.12440 

•19999 

- X'dl— 

1.99900 

•90900 

•94094 

.41l9f" 

—  144*7 

-.91997 

•99909 

•99992 

•94491 

•94194 

FIGURE  7  ANALYSIS  OF  HAND  GEOMETRY  FEATURES  USING  WITHIN-GROUPS  CORRELATION  MATRIX 
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produced,  we  do  not  believe  that  the  data  processing  would  be  limiting. 
However,  if  line  processing  takes  longer  than  scanning,  multiple  proces¬ 
sors  could  be  used  with  parallel  multiple  lines  of  processing  working 
concurrently  to  process  more  than  one  line  at  a  time.  The  need  for 
multiple  processing  has  not  been  assessed,  since  all  work  has  been  per¬ 
formed  on  a  multi-user  system  with  relatively  slow  development-type  pro¬ 
grams.  However,  because  our  software  could  be  supported  by  multiple 
processors  working  in  parallel,  process  time  should  not  be  rate-limiting. 

Since  the  potential  of  this  technique  has  been  shown  to  be  good, 
the  next  task  is  to  refine  the  technique.  Refinements  can  be  made  in 
all  areas,  including  fine-tuning  feature  extraction,  adding  new  features 
and  removing  the  less  significant  ones,  and  developing  methods  for  making 
authentication  decisions  based  on  these  data.  In  addition,  studies 
could  be  made  of  methods  for  optimizing  the  software  for  minimum  through- 


VI  COMPUTER  PROGRAM  SUMMARIES 


The  task  of  image  analysis  has  been  broken  down  into  four  discrete 

jobs: 

•  Determine  the  correct  gray  level  for  image  thresholding  using 
the  FORTRAN  program  HIST. 

•  Threshold  the  image  to  produce  a  binary  image  using  the  FORTRAN 
program  THRESH. 

•  Extract  a  list  of  hand-perimeter  points  from  the  binary  image 
using  the  PASCAL  program  LINEBYLINE. 

•  Extract  the  feature  values  from  the  perimeter  list  using  the 
FORTRAN  program  SCAN. 

Each  of  these  jobs  is  handled  by  a  single  independent  program.  A  brief 
description  of  each  program  follows,  and  listings  of  each  program  appear 
in  Appendix  A.  The  program  LINEBYLINE  is  quite  complex,  however,  and 
Appendix  B  is  included  to  describe  its  operation  in  detail. 

A.  HIST 

This  FORTRAN  program  computes  a  histogram  of  the  frequency  of  occur¬ 
rence  of  gray  levels  in  an  8-bit  gray-scale  image.  The  program  is  capable 
of  handling  a  512  x  512  pixel  image,  which  is  the  size  used  to  represent 
the  hand  data.  The  image  format  is  two  pixels  per  16-bit  word.  The  record 
size  on  a  DEC  RK  05  disk  is  256  16-bit  words.  The  output  file  lists  the 
frequency  count  of  each  gray  level  (0-255).  A  plot  of  the  output  can  also 
be  produced  on  a  TEKTRONIX  4000-series  terminal,  using  the  PLOT  10 
(TEKTRONIX)  software  interface  installed  in  the  Bioengineering  Research 
Center's  PDP  11/40  computer.  The  output  of  this  program  clearly  identi¬ 
fies  the  sets  of  gray  levels  that  represent  either  the  background  or 
the  hand;  this  information  is  used  to  pick  a  discrete  gray  level  between 
the  two  sets  of  levels  that  will  best  separate  them.  This  level  is  used 
by  THRESH  to  produce  a  binary  image. 
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B.  THRESH 


This  FORTRAN  program  thresholds  an  8-bit  gray-level  image  to  produce 
a  binary  image.  Thresholding  is  performed  on  a  pixel-by-pixel  basis. 

The  input  file  contains  two  8-bit  pixels  per  word  with  a  total  image  size 
of  262,144  pixels.  The  output  file  contains  16  one-bit  pixels  per  word. 
Both  input  and  output  files  have  256  words  per  record.  The  discrete 
threshold  value  is  entered  at  run  time. 

C.  LINEBYLINE 

This  PASCAL  program  produces  a  list  of  perimeter  points  that  bound 
"blobs"  of  color  opposite  to  that  of  the  background  in  a  binary  image. 
Eight-point  connectivity  is  checked  for  each  point  in  the  image.  Using 
the  large  area  as  a  criterion,  the  hand  blob  is  selected  from  others, 
and  only  the  hand  perimeter  is  saved.  This  list  of  hand-perimeter  points 
is  output  as  a  file  of  ordered  points  running  around  the  perimeter  of  the 
hand.  Only  one  pass  is  made  through  the  image,  and  only  one  line  is 
accessed  at  a  time.  The  program  therefore  requires  only  slightly  more 
dynamic  memory  than  is  required  to  store  perimeter  points.  A  more  exten¬ 
sive  description  of  this  program  is  included  as  Appendix  B. 

D.  SCAN 

This  FORTRAN  program  processes  the  list  of  perimeter  points  produced 
by  LINEBYLINE  and  outputs  numeric  values  that  describe  geometric  features 
of  the  hand,  as  discussed  earlier  in  this  report.  Included  below  are 
comments  intended  to  help  the  reader  interpret  the  program  listing  in 
Appendix  A. 

•  The  preferential  landmarks  on  a  perimeter  are  the  fingertips. 
These  are  located  by  the  angle  algorithm  on  page  A-20, 
starting  with  the  comment  lines.  Landmarks  are  found  by 
moving  a  window  of  size  IDIS  (18  points)  along  the  perimeter 
and  looking  for  an  angle  criterion  shift  of  90°  in  the  tangent 
to  the  perimeter. 

•  Starting  at  line  174  on  page  A-5,  calculations  of  finger  length 
are  made.  These  require  the  precise  tip  locations  calculated 
earlier  on  page  A-4,  starting  with  the  comment  lines  in  the 
middle  of  the  page. 
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Finger  widths  are  calculated  with  the  code  on  page  A-21. 

Ratios  of  finger  lengths  to  widths  are  calculated  on  page  A-22. 

Curvature  measures  are  calculated  with  the  code  in  the  middle 
of  page  A-22,  using  subroutines  TANAG  and  LINFIT  that  process 
angle-function  data. 

Hand  width,  perimeter,  and  hand  area  are  calculated  at  the  bottom 
of  page  A-22,  using  the  functions  IDELTS  and  DSTNCE  and  the  sub¬ 
routine  ARPRIM. 


Appendix  A 

COMPUTER  PROGRAM  LISTINGS 


A- 1 


oooo  nnooo  floooonooi/i  m  .  on  on  on  o  oo  oorjnoonrionnn 


PROGRAM  HIST 


WITTEN  BY  GREG  MYERS 

THIS  ROUTINE  COMPUTES  A  HISTOGRAM  OP  AN  8-BIT  GREY-LEVEL  IMAGE.  IT 
HANDLES  IMAGES  THAT  ARE  2S6.  512,  OR  1024  POINTS/LINE.  THE  IMAGE  FORMAT  IS 
2  PIXELS  PER  WORD,  AND  THE  RECORDSIZE  ON  DISK  IS  256  16-BIT  WORDS . 

THE  OUTPUT  FILE  LISTS  THE  FREQUENCY  COUNT  OF  EACH  GREY  LEVEL  (0-2SS) . 

NOTE:  TOR  256  X  256  FILES,  THERE  IS  ONLY  ONE  PIXEL  PEI  WORD I 


LOGICAL*  1  NULL,  NAMINO0) ,  N ABOUT (301  ,  A<512),  ITITLE<72) 
INTEGEN*4  1HI256) ,  IHAKEA 

REAL  G(2S6>,  H(2S6>,  XAXISU8) ,  YAXtS<2> ,  TITLEUBI 
INTEGER  STaRTL,  STARTP,  ENDL,  EH DP 
EQUIVALENCE  UTITLE. TITLE) 

DATA  XAXISl  P/'GREY  /,XAX1S(2>/’  LEV V,XAXIS(3)/’EL  V, 
l  YAXISl 1>/’C0UN'/,YAX1S12>/’T  V 


LUNIN  •  3 
LUNOUT  •  4 
NULL  •  0 


DO  I  I  •  4,18 
TITLE) I)  •  4H 
I  XAXIS(I)  •  4H 

— GET  DATA  FROM  TERMINAL 


URITE  (5,10) 

10  FORMAT ( ’  FOR  INPUT  FILE:  TYPE  (LINES,  (POINTS/LINE,  AND  FILENAME') 
READ  (S.1S)  NL,  NP,  NCHIN ,  (NAMINl  I) ,  !•  I  ,NCHIN> 

15  FORMAT (215,Q,30A1 ) 

NAMIN (NCHIN* 1 )  •  NULL 

OPEN(UNIT*UINlN,NANE-NAMIN,TYPE«'OLD' .ACCESS* 'DIRECT' .READONLY, 

I  RECORDS IZE* I 28, £ER«90I) 


WRITE  (5,20) 

»  FORMAT)’  FOR  HISTOGRAM  AREA:  TYPE  STARTING  LIME,  STARTING  POINT,  t 
ILINES,  APOINTS/LINE'I 

read  (5,25)  STARTL,  STARTP,  ML,  MP 
25  FORMAT (415) 

ENDL  •  STARTL  *  NL  -  I 
ENDP  «  STARTP  *  MP  -  I 


WRITE  <S  30) 

»  FORMAT) ' 'TYPE  ''l”  IF  YOU  RANT  A  HISTOGRAM  OUTPUT  FILE  ONLY,’/ 

1  '  '  '2"  IF  YOU  WANT  A  PLOT  ONLY,’/ 

2  •  "3"  IF  YOU  WANT  BOTH’) 

READ  (5.3S)  ISWP 

35  FORMAT  (ID 

IF  HSVP  .EX}.  I  .OR.  ISVP  .EQ.  3)  WRITE  (5,40) 

40  FORMAT)'  TYPE  HISTOGRAM  OUTPUT  FILE’) 

IF  (ISWF  .EQ.  I  .OR.  ISVP  .EQ.  3) 

I  READ  (5,451  NCBOUT, (NAMOUT( I) ,1*1 .NCHOUT) 

IS  FORMAT (Q.36A I > 

HAMOUT (NCHOUT* 11  *  NULL 

IF  (ISVP  .EQ.  2  .OR.  ISVP  .EQ.  3)  VRITE<S,50) 

9  FORMAT ( '  TYPE  0  OF  GREY  LEVELS  DESIRED  IN  THE  PLOT  (13)’) 

IF  (ISVP  .EQ.  2  .OR.  ISVP  .EQ.  3)  READ(S.SS)  NGREYL 

5  FORMAT (13) 


— INITIALIZE  HISTOGRAM  ARRAY 


100 


TIN  •  SECNDS(0.1 
DO  100  I  •  1,256 
iH( i>  •  a 


IP  (NP  .03 
IF  (NP  .55 
IF  (NP  .EQ 
GO  TO  902 


1024)  OO  TO  200 
512)  00  TO  390 

256)  00  TO  400 


— FOR  1024  POINTS/LINE  (2  RECORDS  PER  LINE) 


A- 2 


0 

D  212 


22* 


219 


C 

C 


DO  210  JREC  •  (STARTL-l)»2*l1(Dfl)t-l>»2»l,2 
IF  <M00> JIBC.26)  .BO.  !>  WHITE  (5,212)  JREC 
FORMAT ( '  JREC  •  ’,14) 

IF  (STARTP  .OT.  SI2)  00  TO  220 

READ  (LUNIN' JREC,  EM -903)  <A(  I)  ,1-1 ,512) 

XI  ■  STARTF 

K2  .  N|N0(ENDP,SI2) 

CALL  HISTLlA.KI ,K2, Ifl) 

CONTINUE 

IF  (ENDP  LE.  512)  00  TO  210 

READ  < LUNIN 1 JREC- I , EM-904)  (Ad) , I-  1 .512) 

XI  •  NAX0(STARTP,S13)  -  SI2 

K2  •  ENDP  -  512 

CALL  HISTL(A,KI,K2, IH> 

CONTINUE 
00  TO  500 


C - FOR  512  POINTS/LINE  (ONE  LINE  PER  RECORD) 

C 

300  DO  310  JREC  •  STARTL.ENDt 

READ  ( LUN IN '  JREC  .EM-903)  (A(  I) ,  I- 1 ,512) 

CALL  HISTLU. STARTF, ENDP.iH) 

319  CONTINUE 
00  TO  500 
C 
C 

C - FOR  256  POINTS/LINE  (2  LINES/RECORD) 

C 

400  DO  410  JREC  •  STARTL.ENDt 

READ  ILUNIN' JREC.ERR-903)  (A(l ) , !•! .512) 


C - CONPRESS  DATA  IN  ARRAY  A  FRO*  DEAXZA  FORMAT 

DO  411  I  •  1,256 
411  A( I )  •  A(2*|-|i 

CALL  H|N7L(A.srARTF,ENDP,IN> 

410  CONTINUE 
C 
C 
C 
C 

C - PROMICE  OUTPUT  FILE  IN  A  FORMAT  ACCEPTABLE  TO  THE  PROGRAM  PLOTDATA 

C 

500  CONTINUE 

TINE  -  SECTfOS(TIN) 

WRITE  (5,5011  TIME 
FORMAT ( 1  TIME  •  '.PI0.3,’  SECONDS') 

CLOSE  (UNIT-LUNIN) 

DO  SOS  I  •  1,256 
fl(l)  -  1-1 
Nil)  -  IH(|) 

CONTINUE 


501 


ir  (ISVP  .NB.  I  .AND.  ISWP  .ME.  3)  GO  TO  600  _ 

OPEN  (|TN  IT-LUNOUT ,  NAME- NAMOUT ,  TYPE- '  NEW'.  ACCESS-  'SEQUENT I AL ' , 
FORM- ' FORMATTED ', EM -905, CARR I AGECONTROL- 'LIST' > 


SI9 


S29 


530 

549 


C 

C 


WRITE  (4.510)  (NAMINd), 1-1,27).  STARTL , ST ARTP . ML , MP 
FORMAT (7X, 'HISTOGRAM  OF  1 ,2751,415) 

N2S6  •  256 

WRITE  (4.S20)  N2S6,  <0d)  ,B< I) .  I- 1 ,2561 

FORMAT ( I4,/,2(EI5.8,2XI ) 

I H AREA  -  0 
DO  530  I  *  1,256 
IHANEA  -  IHAREA  *  IBd) 

WRITE  (4,540)  IHAREA 
FORMAT)/,'  HISTOGRAM  AREA 
CLOSE  (UNIT-LUNOUT) 


' ,  17) 


C — PLOT  THE  DATA  ON  THE  TEKTRONIX 
C 

600  CONTINUE 

IF  (ISWP  .LT.  2)  GO  TO  999 
ENCODE  (67.510,  IT ITLE. EM-9) 


C 

C 

C 


UTVWC  'VliJIU,  III  iniwifllll  I 

CALL  PLOT(0,H,NGREYL,XAXIS, TAXIS, TITLE) 
STOP 


(NAMINd)  ,1-1,27)  .STARTL, ST  ARTP,  NL.IIP 


90) 
91 1 


’ ,  39A I  > 


912 


913 


WRITE  (5,911)  (NAMINd)  .  I- 1  .NCHIN) 

FORMAT l '  HIST:  ERROR  OPENING  INPUT  FILE 
GO  TO  999 
WRITE  (5,912) 

FORMAT ( ’  HIST  ERROR:  •  POINTS/LINE  256.  512,  OR  1024  ONLY! ' ) 
CLOSE  (UNIT-LUNIN) 

00  TO  999 
WRITE  (5  913) 

FORMAT) ' *HIST  ERROR:  ERROR  903  READING  INFUT  FILE') 

CLOSE  (UNIT-LUNIN) 

00  TO  999 


A- 3 


oorjooonooooo  oonnononoononnoono 


964 

914 


90S 

915 

9 

916 


C 

C 


166 


999 


WRITE  (5,914) 

FORMAT l '  HIST  ERROR:  ERROR  964  READING  INPUT  FILE') 
CLOSE  (UNIT-LUN1N) 

00  TO  999 

WRITE  (S.9IS)  (NANOVT ( I ) , I* I .NCHOUT) 

FORMAT  I '  HIST  :  ERROR  CREATING  OUTPUT  PILE  \38At> 
00  TO  999 

WRITE  (5,916)  ( I TITLE* I) , 1*1 ,72) 

FORMAT (72A I) 

PAUSE  916 
GO  TO  999 

END 

SUBROUTINE  fllSTL(A,H  ,K2,  IR> 


LOGICAL*!  A(512) 
INTEGER*4  IH<256) 


DO  100  I  •  E1.X2 
I  NT  •  Ad) 

IF  (INT  .LT.  6)  INT  •  INT  ♦  2S6 
1NDX  ■  INT  ♦  I 

IF  ( I H < INDX)  .EQ.  32767)  GO  TO  166 
IH(INDX)  •  IH(INDX)  ♦  I 

CONTINUE 


RETURN 

END 


PROGRAM  THRESH 
WHITTEN  BY  ORB)  MYERS 

THIS  PROGRAM  THRESHOLDS  AN  8-BIT  GREY-LEVEL  IMAGE  AND  PRODUCTS  A  BINARY 
(1  BIT  PER  PIXHpi OUTPUT  IMAGE.  THE  INPUT  FILE  CONTAINS  l  9-BIT  PIXEL 
PER  WORD  IF  THE'  IMAGE  IS  IN  THE  256  X  256  DE  ANZA  FORMAT.  OTMERVIZE, 
THE  IMAGE  CONTAINS  2  PIXELS  PER  WORD.  THE  OUTPUT  FILE  CONTAINS 
16  I -BIT  PIXELS  PER  WORD.  BOTH  THE  INPUT  AND  OUTPUT  FILS  HAVE 
256  WORDS  PER  RECORD. 

MODIFIED  9/24/90  TO  ACCEPT  COMMAND  STRING  INPUT 
COMMAND  SYNTAX  IS: 

IOUTFILI -INF1L/LN : * : */TH:i 
ROB  HANNON 

LOGICAL* 1  NULL,LFL0,NANFIL(34) ,A(8I92) ,B(512) 

INTEGERS  I  BYTE 
LOGICAL* I  BYTE 

INTEGER  lFLG.NL.KPJTHR.EXSTAT 
INTEGER  SVTAB(9> ,THVAL(S> ,LNVAL(5> 

REAL  PROMPT 

DATA  PROMPT/' THR'/, IEXT/3R£XT/,EXSTAT/1/ 

DATA  NL/-1/.NP/-1/.ITHR/-1/ 

EQUIVALENCE  (BYTE, IBYTE) 

DATA  LUNIN/3/,  LUNOUT/4/,  NULL/9/ 


— GET  DATA  FROM  TERMINAL 
WRITE  (5,16) 

16  FORMAT) ’  FOR  INPUT  FILE:  TYPE  ALINES,  fPOINTS/LINE,  AND  FILENAME’) 
RF.AD  (5.15)  NL,  NP,  NCHIN,  (NAN1N(  1) ,  1-1  ,NCHIN) 

16  F0RNAT(2IS,0,36A1) 

NAMIN(NCHIN* l)  •  NULL 

iMitialis*  GETCMD  »»4  iwitoh  doieriptor  litbltl 

CALL  ASSIGN  <G,’TI:’,3>  tAfiifi  LUN  for  GETCND 

CALL  INICND  (6, I EXT, PROMPT)  ISol  LUN  tad  Til*  oit.proBpt 

CALL  CSISW  (SVTABI I  > ,  'LN '  ,2,  ,LNVAL)  IDitlari  AM  *wl  teh 

CALL  CSISW  (SVTAB(S) , 'TH' , I , ,THVAL>  IDooloro  /TH  (witch 

CALL  CSISV  (LNVALdl , 'D' ,NL)  IDieliri  flrit  >*l  for  AMit 

CALL  CSISV  (LNVALO) ,  'D'  ,NP>  IS.co.d 

CALL  CSISV  (THVAL(1),'D',ITHR)  ID.ol.r.  »»1  for  /TH:» 


C  Gat  tOMtU  string  u<  proim  it 
16  CALL  GETCND  (B,.,.,LFLG> 

IF  (LFLO  NE.  e>  GO  TO  1000  (Error 
EXSTAT  *  I  (Sot  successful 

CALL  CSII <B, , , 1FL0)  ICoapress  string 

IFIIFLO  .CT.  0)  00  TO  1020  (Error 

IFtlFLG  .EQ.  Oi  LFLO  •  'I'  (Equal  sign  is  string 

IFUFLC  .LT.  0)  LFLO  •  ’O'  (No  equal  sign 

CALL  CSID£F<LFLG, , , , *H8  ')  (S.l  default  sst 

NL  •  -I  (Sst  Tor  chock 

NP  •  -1 
I TUB  •  -1 

CALL  CSI2(LPL0,NAMFIL,,8VTAB)  (Pars,  string 
IF<I<AWML<II  .BO.  0)  00  TO  1040  (Error 

I F  <  <  NL  .EO.  -I)  .0«.  <NP  .EO.  -*»>  00  TO  1030  (Error 
IFtlTHR  .BQ.  -1>  GOTO  1030  (Error 


(Parse  string 


OPEN  <  UN  I T*LUN 1 N , MAME*NAMF I L , TYPE* 1 OLD ' .ACCESS* 'DIRECT' .READONLY, 
RECORDS 1ZE* 128,ERR*90l> 


TOTPIX  «  FLOAT <NL>  •  NP 
NRECIN  •  TOTPIX  /  SI2 
NBLOUT  •  TOTPIX  /  4096 
L1NC  •  l 
rroT  •  a 

IF  <NL  .NE.  256  .OR.  NP  .NE.  256)  00  TO  18 


NRECIN  • 
LINC  *  2 
KTOT  •  16 
CONTINUE 


NRECIN  •  2 


C  Gat  iaforaatioa  ead  open  oatpat  file 

IF  (LFLG  .BQ.  *  I ’ >  GO  TO  50  (Eqael  sign  in  string 
NANFIL<25>  •  'B'  (Sat  default  elteasioa 

NANFIL<26)  •  •  1  * 

NA«FIL<27>  •  'N' 

GO  TO  66 

50  CALL  CSIDEFt 'O' , , .NANFILC IS) , 'BIN')  tSet  default  ext 
CALL  CSI2< 'O' .NAJIFIL)  (Parse  string 

lFlNANFILM)  .EO.  0)  GO  TO  1040 
C  WRITE  (S  20) 

C  20  FORNATt ' ’TYPE  THRESHOLD,  AND  NAME  FOR  BINARY  OUTPUT  FILE') 
C  READ  <5,25)  ITHR, NCHOUT , (NAHOUTt I) , I* 1 .NCHOUT) 

C  2S  FORNAT <  14,0,3041 ) 

C  NAMOUT (NCHOUT- I )  -  NULL 

60  0PENIUN1T*LUN0UT .NANE-NAMFIL, TYPE* 'NEW .ACCESS* 'DIRECT' , 

I  RECORDS I ZE- 128,ESR*902, 1NIT1ALSIZE*NBL0UT) 


(Set  default  ext 
(Parse  string 


TIH  -  SECNDSC0.) 

DO  100  J  •  I .NRECIN, KTOT 

C 

C - READ  IN  4096  8-BIT  PIXELS 

C 

DO  200  X  •  I .KTOT 

XOFF  •  <  X- 1 ) *5 1 2 

READ  tLUNIN' 1*K- l)  <A<KOFF* I ) , 1*1 ,512) 

200  CONTINUE 

C 
C 

LOFT  •  1 

DO  300  L  •  1,512 
IBYTE  ■  0 
C 

C - PERFORM  THE  THRESHOLD  IN  GROUPS  OF  8  PIXELS 

C 

DO  400  I  *  1,8 

l NT  •  A<LOFF) 

LOFF  *  LOFF  ♦  LINC 

IF  < I NT  .LT.  0)  INT  •  INT  ♦  256 
ibyte  ■  ishf t < ibj te .  —  1 ) 

if  list  ,ge.  ithr)  ibyte  •  iordby te,  128) 

400  con  I  inns 

b 1 1 )  ■  byte 

300  ceatiaas 

C 

e-- •rite  ant  one  record  of  data  (4096  binary  pixels) 
a 

•rite  < laaoa t ’ j/k tot*  1)  (b < I ) , i* 1 ,512) 

100  aaatiaaa 
a 
a 

CLOSE  (UNIT  •  LUNIN) 

CLOSE  (UNIT  •  LUNOUT) 

TINE  •  SECNOS(TIN) 

MUTE  (5,991)  TINE 


a 


999  00  TO  10 


10*  t  Mlt  OOMIld 


e 

e 

*901  writ*  (5,911)  (NAMFILI 1) , i« 1,341 

911  formate  orror  opoolit  iopot  rilo  ’,30«1) 

•o  to  999 

902  «ri  to  (5,912)  (NANFILU ) ,  i- 1 ,34) 

912  formate  orror  oroatiap  ootpot  filo  ',30al> 

§o  to  999 

C 

1000  IF1LFL0  .EX).  "366)  CALL  EXST (EXSTAT)  lEiit  with  ttataa 

1020  TYPE  1021 

1021  FORMAT  e  •••  command  (fatal  orror  •••') 

00  TO  9000 

1030  TYPE  1031 

1031  FORMAT  e  •••  orror  to  iapit  filotpoo  ••••) 

00  TO  9600 

1040  TYPE  1041 

1041  FORMAT  e  •••  Error  ia  oatpit  fitoapoo  •••’) 

00  TO  9000 

C 

9000  EXSTAT  ■  4  ISot  for  aotor  orror 

—  CALL  RESCMD  tRoiot  cad  iapat  lo  top  lovol 

00  TO  10  tOot  aoxt  ooammnd 

a 
e 

and 


PROGRAM  LINEBYLINECTTY) ; 


(•NSP*) 


<•  VERSION  46C  10/16/80 

THIS  VERSION  PRODUCES  A  PERIMETER  LIST  OF  ONE  LARGE  BLOB,  WHICH  IS  SUPPOSED 
TO  BE  A  HAND.  FOR  ROB  HAN  MON  AND  THE  PALMPRINT  STUDY. 

WRITTEN  BY  GREG  MYERS 

SRI  INTERNATIONAL,  4-3006 
333  RAVENSVOOD  AVENUE 
MENLO  PARK,  CA.  9402S 


THIS  ALGORITHM  IS  TAKEN  FROM  A  REPORT  BY  OERRY  AGIN  AT  SRI  ENTITLES 
IMAGE  PROCESSING  ALGORITHMS  FOR  INDUSTRIAL  VISION".  THIS  ALGORITHM 
IS  REFERRED  TO  AS  CONNECTIVITY  ANALYSIS".  IT  SEGMENTS  A  BINARY 
IMAGE  INTO  "BLOBS"  (CONNECTED  AREAS)  OF  THE  SAME  COLOR"  (OREY  LEVEL). 

ONLY  ONE  PASS  IS  MADE  THROUGH  THE  IMAGE.  AND  ONLY  ONE  LINE  IS  ACCESSED 
AT  A  TIME  (HENCE,  THE  PROGRAM  NAME  ’LINE-BY-LINE*  >  .  FEATURES  OF  EACH 
BLOB  ARE  COMPUTED,  SUCH  AS  ITS  AREA,  CENTER  OF  GRAVITY,  BOUNDING  RECTANGLE, 
AND  A  LIST  OF  ITS  PERIMETER  POINTS. 


THE  LETTERS  NSf '  WITHIN  A  COMMENT  MEAN  'NON-STANDARD  PASCAL*. 

THE  LINES  MARKED  *NSP'  MAY  REQUIRE  MODIFICATION. 

THE  LETTERS  <MD*  WITHIN  A  COMMENT  MEAN  'MACHINE-DEPENDENT'. 

THESE  LINES  MAY  ALSO  REQUIRE  MODIFICATION. 

BECAUSE  THE  COMMAND  'DISPOSE'  IS  NOT  IN  THE  "SWEDISH"  VERSION  OP  PASCAL 
AT  SRI,  'DISPOSE'  STATEMENTS  ARE  ENCLOSED  IN  COMMENT  STATEMENTS  AND 
ARE  IGNORED. 

READ  AND  WRITE  STATEMENTS  FROM  THE  USER’S  TERMINAL  DO  NOT  REQUIRE 
FILE  SPECIFICATION  IN  THIS  VERSION  OF  PASCAL.  THE  FILE  'TTY'  IS 
ASSUMED. 

EOLN  IS  AT  THE  BEGINNING  OF  AN  INPUT  LINE  WHEN  READING  FROM  THE  TTY. 

•) 


CONST  NPIXELSPERLINE  •  512;  (•  THIS  PROCRAM  HANDLES  512  X  512  IMAGES  ONLY  •) 

NUNES  •  512; 

NBITSPERCHAR  •  8;  (*MD* ) 

NCIURSPKRUNE  •64;  (•  •  NPIXELSPERLINE  /  NBITSPERCHAR  •)  (*MD*) 

NCHARPERBLOCK  •  SI2;  (*MD*> 

BKGND  *0;  INK  •  1 ; 

(•  INPUT  COMMANDS  •) 

ON  •  ’T’;  OFF  •  'F'; 

LINEIN  •  'L'l 
LISTBLOBS  •  ‘B’: 

LISTPER IMS  •  'P1; 

LISTRUNUNOTHS  ■  'R': 

LISTACTIVELINESEOMENTS  •  'A'; 

DIAGNOSTICS  •  'D'; 

STOP  •  'S’; 


A-6 


TYPE  BINARY  •  BJCGND. .  INK 
CWJTYPE  •  SET  OF  'A 


VMV 1  lib  -  Obi  WT 

ARRAYS I NT  •  ARRAY ( I . .5]  OF  INTEGER; 


BLOBPTR  •  (BLOBS; 

PTRPER  INSECT  ION  •  (PER  INSECT  ION; 

PERINPTR  •  (FERIM; 

BLOBS  •  RECORD 

COLOR;  BINARY; 

COMP:  INTEGER; 

PERIN:  PTRPER INSECT  I ON; 
NPERINPTS:  INTEGER: 

AREA , XNEAN , YNEAN ;  REAL; 
XNIN.XNAX.YNIN.YNAX:  INTEGER; 
PARENT,  NEXT:  BLOBPTR  END; 


PER INSECT  ION  •  RECORD 

LEFT, RIGHT:  PERINPTR; 

PREV .NEXT :  PTRPER INSECT  ION  END; 


PERI*  •  RECORD 

LINE, COL:  INTEGER; 
NEXT:  PERINPTR 


END; 


SEGPTR  •  (SEGMENT; 

SEGMENT  •  RECORD 

ST ARTCOL , ENDCOL :  INTEGER; 
BLOB:  BLOBPTR; 

NEXT:  SEGPTR  END; 


PTRRUNLENOTH  •  (RUNLENGTHS; 

RUNLENGTHS  •  RECORD 

ST  ARTCOL, ENDCOL:  INTEGER; 
NEXT:  PTRRUNLENOTH  END; 


(•  GLOBAL  VARIABLES  •> 

VAR  BLOB .BLOBSDONE :  BLOBPTR; 

RECYCLEDPTR :  PERINPTR; 

ARRAYP:  ARRAY (1.. I II  OF  PERINPTR j, 

ACT  I VEL I NE ,  PREVSEG , CURRSEG , SEG :  SEGPTR; 

NEWLINE, NEVSBG.LASTNEVSEG:  PTRRUNLENGTH ; 

CND.ONOFF:  CHAR; 

CNDSET :  CNDTYPE; 

I NT KR ACTIVE:  BOOLEAN* 

L I NENUN , COL i NEWCONPNUN , NLREADJfT I NES :  INTEGER: 

TH  ACEBLOBS ,  TRACEPER I  NS ,  TRACEACT I  VEL  I NESEGNENTS ,  THACERUNLENOTHS:  CHAR 
TRACED UGNOST ICS:  CHAR: 

INPUTINAGE. DATAFILE. PERiNFILE:  TEXT; 

INPUTNANE, DATANANE , PER I NNANE :  ARRAYI1..30)  OF  CHAR;  (»ND»> 

BIT:  ARRAY! I • . NPIXELSPERLINEI  OF  BINARY; 

I.CHARTOT.L1NEGROUP:  INTEGER; 

TYPELINENUNS:  BOOLEAN;  _ 

HAN DON LY , BLOBWR ITTEN :  BOOLEAN;  AREATH;  REAL;  IAREATH:  INTEGER; 
TRACEPARENT, RECURS:  BOOLEAN; 

NSKIP:  INTEGER; 

(••P»  NEW  PAGE  •> 

PROCEDURE  PAUSE; 

VAR  DUNNYCHAR:  CHAR; 

BEGIN  _ 

VRITEUM 'TYPE  ANY  CHARACTER  TO  CONTINUE'); 

READLN; 

READ l DUNNYCHAR I 
END;  (•  PAUSE  •> 


PROCEDURE  TYPERUNLENGTHS; 

VAJ  jILPTRi  PTRRUNLENGTH ; 

END;  (•  TTPERUNLENOTBS  •) 

PROCEDURE  TYPEACT I VEL I NESEGNENTS ; 

VAR  ALSPtR:  SEGPTR; 

BEGIN 

END;  (•  TYPEACT I VEL I NESEGNENTS  •) 


A- 7 


««8P*  HEW  PAGE  *1 

PROCEDURE  VK1TEBLG6(VAR  DEVICE:  TEXT;  VAX  BLOB:  BLOBPTR) 


BEGIN 

END;  <•  VRITEBLOB  •) 


PROCEDURE  TYPEBLOBSCVAR  DEVICE-TEXT) ; 

VAR  ALSPTR:  SEGPTR; 

BEGIN 

END;  «•  TTPEBLOBS  •) 

(•tP»  NEW  PAGE  •) 

PROCEDURE  DIRPOINTIVAR  DIR.Xl ,Y1 ,X2, Y2:  INTEGER); 


«• 


-1 


DIRECTIONS: 

K 

e  1 
2 


-i 

3 


Y  0  4  . 

1  S  6 

0 

7 

BEGIN 

CASE  DIR  OP 

3,2,1:  Y2 

■ 

Yl-l; 

4.0  :  Y2 

■ 

Yl; 

S ,6,7:  Y2 

■ 

Yl*l 

END; 

CASE  DIR  OF 

3.4.S:  X2 

■ 

X 1  —  l ; 

2.6  :  X2 

■ 

XI; 

0,1.7:  K2 

• 

Xl«l 

END; 

END;  <*  DIRP01NT 

•) 

•) 


PROCEDURE  UNPACKS (VAR  XX:  INTEGER;  VAR  DIR:  ARXAVStNT) ; 


BEG  I A 


,X,XDIV4096,XDIVS12,XDIV64,XD1V8:  INTEGER; 


X  :*  -XX; 

XDIV4096 

XDIVSI2 

XD1V64 
XDIV8 
DIR  I  II  ;• 
DJRI2) 


•  X  DIV  4096; 

•  X  DIV  S12; 

-  X  DIV  64; 

•  X  DIV  8: 

X  -  XDiV8*8: 

XDIV8  -  XD1V64*8: 


vini4j  nv  —  auitwo 

D I R I 3 1  :•  XDIV64  -  X01VS12»8: 

DIR  14)  :•  XDIVSI2  -  XDIV4096*8; 

DIR  I  SI  :•  XDIV4096; 

END;  <•  UNPACKS  •) 

(•8P»  NEW  PAGE  *> 

PROCEDURE  UNPACKPERINS4V4R  PTR.UPPTRR, UPPTRL:  PERINPTR); 

VAR  DIR:  ARRAVSINT; 

1,1 1:  INTEGER; 

BEGIN 

ARRAYPlIlf .LINE 
ARRAYPI I 1 » .COL 


•  PTR* -LINE; 

•  PTRt.COL  -  NPIXELSPERLINE; 


UNPACK5(PTR* .NEXT* .LINE, DIR) ; 

FOR  I  :•  I  TO  S  DO  BEGIN 
1 1  :■  I  *  Is 

DIRPOINTfDIRI  II  .ARRATPIll*  .COL,  ARRAYPI  (I  t.LINE.ARRA  TPIIt)*.COL,  AURA  TPIIDt.LINE); 
END;  <•  FOR  •) 


UNPACK54PTR* .NEXT* .COL,  DIR) 

IF  TRACEDI AGNOSTICS  •  ON  THEN  FOR  1  :•  !  TO  S  DO  WR1TELNIDIRI II :4) ; 

FOR  I  :•  6  TO  10  DO  BEGIN 
11  :•  1  ♦  l; 

D1RP0INT(DIRI 1-S), ARRAYPI 11*. COL, ARRAYPI 1 1*.  LINE, ARRATPI 111*. COL, ARRAYPI I II*. LINE) 
END;  (•  FOR  •) 


UPPTHR  :•  ARRAYPI II; 

UPPTRL  :■  ARRAYPI II); 

END;  4*  UNPACKP ERIKS  •) 

<*8P»  NEW  PACE  •) 

FUNCTION  CONNECT ED (VAX  POINT! .P0INT2:  PERINPTR):  BOOLEAN; 


BEGIN 

IF  *ABS(P0INT1*  .LINB-POINT2*  .LINE)  <•  l>  AND 

*ABS( POINT  1 * . COt  -POl NT2* .COL  )  <■  1)  THEN  CONNECTED  :«  TRUE 
ELSE  CONNECTED  :•  FALSE; 

END; 


PROCEDURE  WRITEPOINTS (VAR  DEVICE:  TEXT; 

VAR  RIGHT.LEFT:  PERINPTR; 

VAR  NPERINPTS:  INTEGER); 

VAR  PTR. PTR2, UPPTRR, UPPTRL:  PERINPTR; 

BEGIN 

(•  IF  i TRACED I AGNOSTICS  •  ON)  AND  RECURS  THEN  BEGIN 
URITELNI  ’RECl)RS»TRU£' ) ; 

PTR2  :•  RIGHT; 

WHILE  PTR2  <>  LEFT  DO  BEGIN 

WRITELN (PTR2*  ,LINE:4,PTR2*  .COL: 4)  ; 

PTR2  :•  PTR2t .NEXT;  END; 

WRITELN <PTR2t . LINE:4,PTR2* .COL: 4) ;  END;  •) 

PTR  :*  RIGHT; 

NPERINPTS  :•  NPERINPTS  *  1; 

WITH  PTR*  DO  WRITELN<DEVICE,LINE:4,’  ’.C0L:4); 

WHILE  PTR  <>  LEFT  DO  BEGIN 

IF  PTR* .NEXT*  COL  >  NPIXELSFERLINE  THEN  BEGIN 

IF  RECURS  THEN  BEGIN  WR I TELN ( ' 2ND  RECURSION') ;  PAUSE  END; 
UNPACKPER I  NS  t  PTR* .  NEXT ,  UPPTRR ,  IIPPTRL)  ; 

I*  IF  TRACED I AGNOSTICS  ■  ON  THEN  WRITELN* ’WRITEPOINTS  —  A’); 

IF  TRACED 1AGNOST ICS  >  ON  THEN 
RECURS  :•  TRUE;  •) 

(•  PTR2  :«  UPPTRR; 

WHILE  PTR2  < >  UPPTRL  DO  BEGIN 

WR ITELN < DEV ICE . PTR2* . L I NE : 4, PTR2* . COL : 4) ; 

PTR2  :«  PTR2* .NEXT ;  END; 

WRITELN  <  DEV  ICE ,  PTR2*  ,LINE;4,PTR2*  .C0L-.4)  ;  •) 

WR  ITEPO I  NTS  (DEVICE,  UPPTRR ,  UPPTRL .  NPER I  NPTS>  s 
(•  IF  TRACED  {AGNOSTICS  «  ON  THEN  WRITELN* ’WRITEPOINTS  --  B'); 

IF  TRACED 1 AGNOSTICS  •  ON  THEN  RECURS  :*  FALSE;  •> 

PTR  ;«  PTR*. NEXT*. NEXT;  END 
ELSE  BEGIN 

*•  IF  NOT  CONNECTED* PTR, PTR*. NEXT)  THEN  BEGIN 
WRITELN* ’POINTS  NOT  CONNECTED’); 

WRITELN (DEVICE. ’POINTS  NOT  CONNECTED’);  PAUSE;  END;*) 

<•  SKIP  UNNECESSARY  PERIMETER  POINTS  •) 

IF  *PTR« .NEXT  O  LEFT*  AND  CONNECTED (PTR. PTR* .NEXT* .NEXT) 

THEN  PTR  :•  PTR*. NEXT; 

PTR  PTR*. NEXT; 

NPERINPTS  :•  NPERINPTS  ♦  I; 

(•  IF  TRACED  I  AGNOSTICS  -  ON  THEN 

WITH  PTR*  DO  WRITELN(LINE;4,  ’  \C0L:4>;  •) 


END; 


WITH  PTR*  DO  WR!TELN<DEVICE,LINE:4,' 
END;  I*  ELSE  •) 

END;  <•  WHILE  *1 
*•  WRITEPOINTS  •» 


,COL:4> ; 


(«*?♦  NEW  PAGE  •) 

PROCEDURE  WRITER ERINS (VAR  DEVICE:  TEXT;  VAR  BLOB: 


BLOBPTR) ; 


•BLoi  COMPONENT  *  \C0NP:4> 


VAR  PSPTR:  PTRPER INSECT  ION; 

BEGIN 

WITH  BLOB*  DO 

IF  FERIN  <>  NIL  THEN  BEGIN 
NPERINPTS  :*  0; 

IF  NOT  HANDONLY  THEN  WRITELN (DEV ICE. ' 

IF  NOT  HANDONLY  THEN  WRITELN(DEVICE. 

PSPTR  :•  PERIN; 

REPEAT 

IF  NOT  HANDONLY  THEN  WRITELN (DEV ICE, ' 

IF  NOT  HANDONLY  THEN  WRITELN(DEV ICE, ’LINE  COL'); 

WITH  PSPTR*  DO  WRITEPOINTS (DEV ICE, RIGHT, LEFT, NPERINPTS) ; 
PSPTR  :■  PSPTR*. NEXT 
UNTIL  PSPTR  •  PERIR; 

END  (•  IF  •) 

END;  (*  WRITEPERINS  •) 


PROCEDURE  TYPEP ERINS (VAR  DEV ICE: TEXT) ; 

VAR  ALSPTR :  SEGPTR; 

BEGIN 

WRITELNIDEVICE, ' 

VR1TELN(DEV1CE, 'BLOB  PERIMETERS’); 
ALSPTR  :•  ACTIVELINE: 

WHILE  ALSPTR  <>  NIL  DO  BEGIN 

WRITEPERINS (DEV  ICE, ALSPTR* .BLOB) ; 
ALSPTR  :•  ALSPTR*. NEXT  END; 

END;  <•  TYPEPERIMS  •) 


(•BP*  NEW  FACE  •> 

PROCEDURE  NEW (V AH  NEWPO INTER:  PERINPTR) ; 


(•NSP»> 


(•  T«tS  ROUTINE  IS  NEEDED  ONLY  IP  THE  COMUND  'DISPOSE'  IS  NOT  AVAILABLE 
IN  THIS  VERSION  OF  PASCAL  •) 

BEGIN 

IP  RECYCLEDPTR  <>  NIL  THEN  BEGIN 
NEWPO INTER  :«  RECYCLEDPTR: 

RECYCLEDPTR  :•  RECYCLEDPTR* .NEXT; 

NEWPO INTER*  .NEXT  :•  NIL  END. 

ELSE  NEW (NEWPO INTER) ; 

END;  **  NEW  •> 


PROCEDURE  DELETEPEXINS4VAR  RIGHT, LEFT:  PERINPTR); 

VAR  POINT:  PERIITTR; 

BEG  IV 

(•  ADD  THE  LIST  OF  PERIMETER  POINTS  OF  THE  BLOB  TO  THE  LIST  OF 
RECYCLED  PERIMETER  POINTS.  THIS  SECTION  OF  CODE  SHOULD 
BE  USED  IF  THE  COMMAND  'DISPOSE'  IS  NOT  IMPLEMENTED  IN  THIS  VERSIOH 
OF  PASCAL  •) 

LEFT*. NEXT  :•  RECYCLEDPTR; 

RECYCLEDPTR  :•  RIGHT; 

(•  DISPOSE  OF  PERIMETER  POINTS  AND  PERIMETER  SECTION.  THIS  SECTION  OF  CODE 
CAN  BE  USED  ONLY  IF  THE  COMMAND  'DISPOSE'  IS  IMPLEMENTED.  •> 

<•  WHILE  BIGHT  <>  LEFT  DO  BEGIN 
POINT  :•  RIGHT; 

RIGHT  :•  RIGHT*. NEXT; 

DISPOSE! POINT)  END; 

DISPOSEIRIGHT)  •> 

END;  (•  DELFT EPERIMS  •> 

(•BP*  NEW  PAGE  •)  _  _ 

FUNCTION  DIRECT  ION (VAR  XI.Yl.XE.YR:  INTEGER) :  INTEGER; 


<•  DIRECTIONS: 
X 

-1  0  1 

-13  2  1 

Y  0  4  .  0 

I  5  6  7 


•) 


VAR  DELTAY:  INTEGER; 

BEGIN 

DELTAY  :•  Y2  -  Yl; 

CASE  X2  -  XI  OF 

-I:  CASE  DELTAY  OF 
-1:  DIRECTION 
0:  DIRECTION 
I:  DIRECTION 


0: 


CASE  DELTAY  OF 
-1:  DIRECTION 
I:  DIRECTION 


I:  CASE  DELTAY  OP 
-I:  DIRECTION 
0:  DIRECTION 
1:  DIRECTION 


3; 

4; 

5 


2; 

6 


I; 

V 


END; 


END; 


END 


END; 


END;  (•  CASE  DELTAY  •) 
<•  DIRECTION  •) 


<*BP*  NEW  PAGE  *)  _ 

PROCEDURE  C0HPACTP0INTS4VAR  PTR.PTRB:  PERINPTR;  VAR  I:  INTEGER); 
Vy^PREVLINE.PREVOOL. POSITION. DIR:  INTEGER; 

1I1TH  PTR*  DO  BEGIN 
IF  I  •  0  THEN  BEGIN 
PTRB* .LINE  :•  0; 

PTRBt.COL  »•  0; 

PTRB*. NEXT  :•  NIL  END 
“  IF  I  <  6  THEN  BEGIN 

IF  I  ■  I  THEN  POSITION  I*  I  ELSE  POSITION  :•  POSITION  •  8: 

PTRB*. LINE  :■  PTRB* .LINE  -  POSITION  •  DIRECT  ION (PREVOOL,PREVLINE,OOL, LINE) ; 
END 


BEOtN 

IF  I  •  6  THBI  POSITION  :•  I  ELSE  POSITION  :■  POSITION  •  ft 

PTJSt.OX  t.  PTRB* .COL  -  POSITION  •  DIRECTIONIPREVCOL.PtiVLIHE, COL, LINE)  s 

END; 

PREVLINE  :•  LINE; 

PEEVCOL  :•  COL; 

END;  (•  WITH  •) 

END;  <•  CO', />  ACT  POINTS  •  ) 

<*«r*  NEW  PAGE  •) 

PDOCEMJHE  PKPERIRStVAR  PEEIH:  PTEP EE INSECT ION I ; 

VAR  PTE ,  PTEA ,  PTEB ;  PEEINPTE; 

FIRSTOROUPOF 10:  BOOLEAN;  <•  NOT  USED  EIGHT  NOV  •) 

NPEX INPTS.lt  INTEGER, 

BEGIN 

PTE  PEEINf .EIGHT ; 

I  :•  0; 

F1ICTGROUPOF10  :•  TRUE; 

IF  PTH  <>  PEEIH! .LEFT  THEN 
REPEAT  BEGIN 

IF  I  PTH*  .LINE  <  UNQHm-1)  AND  (PTE*  .COL  >  0) 

_  _  AND  <PTH*.OOL  <-  NPIXELSPERLINE)  THEN  BEGIN 

IF  I  •  0  THEN  BEGIN 
PTEA  PTH; 

NEWlPTRB)  END; 

OONPACTPO I NTS < PTR . PTRB , I ) ; 

It*!*  1; _ 

IF  I  •  II  THEN  BEGIN 

(•  TO  IDENTIFY  THAT  COUP ACTED  RECORDS  FOLLOW  *) 

PTR A* .COL  : «  PTRA* .COL  ♦  NPIXELSPERLINE; 

«»  INSERT  COMPACTED  PERINETER  RECORD  INTO  THE  LIST  •> 

PTRB* . NEXT  :»  PTR*. NEXT; 

DELETEPF.RINSIPTRA*  .NEXT, PTH) ; 

PTRA*. NEXT  :•  PTRB; 

PTR  :*  PTRB; 

I  :»  0; 

END;  (•  IF  I  •  II  •) 

END  (•  IF  3  CONDITIONS  THEN  •) 

ELSE  IF  I  >  0  THEN  BEGIN 

<•  AN  INCONPLETE  CONPACTION  OOCURED  --  RESET  VARIABLES  •) 

I  :•  0; 

FIRSTGROUPOFJ0  TRUE; 

DEUtTEPERlNS(PTRB.PTRB)  END; 

PTH  t-  PTR* .NEXT; 

END  (•  REPEAT  BEGIN  •) 

UNTIL  PTR  -  PERIH* .LEFT t 
END;  <•  PAPER  INS  •) 


PROCEDURE  PACKFERIXS; 

VAR  NPERIHPTS:  INTEGER; 

BEGIN 

SEO  :•  ACTIVELINE; 

WHILE  SEG  <>  NIL  DO  BEGIN 
WITH  SEG*. BLOB*  DO 

IF  PERIH  <>  NIL  THEN  BEGIN 
PKPERIMS(PERIN) ; 

PEKIH  :•  PERIH*. NEXT  END; 
SEG  :•  SEG* -NEXT  END; 

END;  !•  PACK PER INS  •) 

l»BP»  NEW  PAGE  •) 

FUNCTION  END0F8L0CX: BOOLEAN; 

BEGIN 

IF  CHARTOT  •  NCHARPERBLOCK  THEN  BEGIN 
CHARTOT  :■  0; 

ENDOFBLOCK  :•  TRUE  END 
ELSE  ENDOFBLOCt  :•  FALSE; 

CHARTOT  :•  CHARTOT  »  l; 

BID;  (•  ENDOFBLOCK  •) 


PROCEDURE  OETLINE; 

(•  GET  A  LINE  OF  B INART  HUGE  DATA  FROM  AN  INPUT  FILE.  •> 

VAR  INPUTBTTE.EOBL:  CHAR: 

I . J, INDX.BITS:  INTEGER; 

BEGIN 


<•  READ  IN  ONE  LINE  OP  DATA  FROM  THE  INPUT  PILE  •> 

<•  UNPACK  THE  BINARY  DATA:  IT  IS  ASSUMED  THAT  THE  DATA  IS  INTERPRETED  AS 
A  SERIES  OF  8-BIT  CHARACTERS  •) 

INDX  :•  0; 

FOR  I  :•  1  TO  NCHARSPERLINE  DO  BEGIN 

IF  ENDOFBLOCK  THEN  READ* INPUT! MAGE, EOBL) i  (*MD*> 

READ* INPUT I MAGE , I NPUTBYTE ) ; 

BITS  ORDtlNPUTBYTE) ; 

FOR  J  ;•  1  TO  NB I TSPERCHAR  DO  BEGIN 
INDX  :•  INDX  •  I; 

IF  ODD(BITS)  THEN  BITIINDXI  :•  INK 

ELSE  BITIINDXI  :«  BKGND; 

BITS  :«  BITS  DIV  2; 

END. 

END; 

END;  <•  GETLINE  •» 

PROCEDURE  SETLINE; 

(•  SET  UP  STARTING  RUNLENGTH  RECORD  •> 

BEGIN 

IF  LINENUN  •  1  THEN  BEGIN 
NEV(NEVLINE) ; 

NEWLINE* .NEXT  :•  NIL  END; 

NEWLINE* .STARTCOL  :•  MININT* 1 ; 

NEWLINE* -ENDCOL  :•  MAXINT-1; 

NEWSEG  :•  NEWLINE; 

LASTNEWSEG  :•  NEWLINE; 

END;  <•  SETLINE  •> 


PROCEDURE  ADDRUNLENGTH (COLUMN:  INTEGER!; 

VAR  NEWESTSEG:  PTRRUHLENGTH ; 

BEGIN 

IF  NEWSEG*. NEXT  <>  NIL 

THEN  NEWESTSEG  :•  NEWSEG*. NEXT  (•  ADVANCE  NEW  SEGMENT  POINTER  IP 

A  SEGMENT  ALREADY  EXISTS  •> 

ELSE  BEGIN 

NEW I NEWESTSEG);  _ (•  CREATE  A  NEW  SEGMENT  IF  IT  DOESN’T  EXIST  •) 

NEWSEG*. NEXT  :■  NEWESTSEG;  (•  CONNECT  NEWEST  SEGMENT  TO  LIST  •) 

NEWESTSEG*. NEXT  :•  NIL  END; 

NEWSEG*. ENDCOL  :-  COLUMN  -  I; 

NEWESTSEG*. ST ARTCOL  :-  COLUMN; 

NEWESTSEG*  . ENDCOL  :«  NAXINT-1:  _ 

NEWSEG  :•  NEWESTSEG  (•  UPDATE  SEGMENT  POINTER  TO  THE  NEWEST  SEGMENT  •) 

END;  <•  ADDRUNLENGTH  •) 


PROCEDURE  CREATERUNLENGTHS; 

<•  FOR  EACH  TRANSITION,  CREATE  A  NEW  RECORD  FOR  THE  RUNLENGTH  SEGMENT  •) 

VAR  I:  INTEGER; 

BEGIN 

IF  BIT! II  -  INK  THEN  ADDRUNLENGTH 1 1 ) ; 

FOR  I  :■  2  TO  NP I XELSPERL I NE  DO  _ 

IF  BIT  1 1- 1 1  <>  HITCH  THEN  ADDRUNLENGTH  ( I )  ; 

IF  BIT INPIXELSPEKLINE)  •  INK  THEN  ADDRUNLENGTH(NPIXELSPERLINE»I) ; 
LASTNEWSEG  :«  NEWSEG; 

IF  TRACERUNLENGTHS  •  ON  THEN  TYPERUNLENOTHS ; 

END;  (•  CREATERUNLENGTHS  •) 

(••P*  MEW  PAGE  •> 

PROCEDURE  ADDRIGHTPER INPOINT (VAR  PER IMrPTRPER INSECT  ION ;  LINE, COL: INTEGER) ; 

(•  THIS  PROCEDURE  ADDS  ONE  POINT  TO  THE  PERIMETER  LIST.  THE  NEW  POINT  IS 
INSERTED  BEFORE  THE  POINT  DESIGNATED  BY  PER  IN*. RIGHT  •> 

VAR  NEVPOINTER:  PERINPTR ; 

BEGIN 

NEW  (NEVPOINTER) ;  (•NSP«> 

NEVPOINTER*. LINE  :•  LINE; 

NEVPOINTER*. COL  :•  COL; 

NEVPOINTER*  NEXT  :•  PERIMt .RIGHT; 

PERIM* .RIGHT  :•  NEVPOINTER: 

END;  (•  AD0R1GHTPER INPOINT  •) 


PROCEDURE  ADOLEFTPER INPOINT (VAR  PERIM:PTRPERiRSECTION;  LINE, COL: INTEGER) ; 

(•  THIS  PROCEDURE  ADDS  ONE  POINT  TO  THE  PERIMETER  LIST.  THE  NEW  POINT  IS 
INSERTED  AFTER  THE  POINT  DESIGNATED  BY  PER IN* .LEFT  •> 
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VAR  NEVPO INTER:  PERINPTR; 

BEGIN 

NEVVlNEVPO  INTER);  (»KSP«> 

NEVPO INTER*  .LINE  :•  LINE; 

NEV POINTER* .COL  :•  COL; 

NEVPO INTER*  .NEXT  :•  NIL; 

PER  IN*. LEFT*. NEXT  :*  NEN POINTER; 

PEN  IN*  .  LEFT  NEVPOINTER; 

END;  <•  ADDLEFTPERINPOINT  *) 

(••P.  NEW  PAGE  •» 

PROCEDURE  I NSERTSEONENT (VAR  ST ARTCOL , ENDCOL : INTEGER) ; 

VAR  SURRBLOB:  BLOBPTR ; 

SUKRCONP :  INTEGER ; 

SUKNCOLOR:  BINARY; 

A  ,B:  SEGPTR; 

COL:  INTEGER; 

NEWPERIN:  PTRPER INSECT  ION; 

BEGIN 

IF  TRACED I AGNOSTICS  «  ON  THEN  WRITELNI ’BEGIN  INSERTION  •>; 

SURRBLOB  :•  PREVSEG* . BLOB ;  <•  SURROUNDING  BLOB  •> 

SURRCONP  :•  SURRBLOB* .CONP: 

SURRCOLOR  :•  SURRBLOB*  .COLOR; 

I*  CREATE  A  NEV  BLOB  FOR  THE  NEW  SEGMENT  •> 

NEW(BLOB)  ; 

WITH  BLOB*  DO  BEGIN 

IF  SURRCOLOR  •  BKGND  THEN  COLOR  :•  INK 

ELSE  COLOR  :  •  BKGND; 

NEVCONPNUN  :•  NEWCONPNUN  ♦  I; 

CONP  :•  NFWCOMPNIIN; 

AREA  :•  0.0; 

KMEAN  ;•  0.0; 

YNEAN  :-  0.0; 

XNIN  :•  ST  ARTCOL; 

XNAX  -  ENOCOL; 

YNIN  .  -  LINF.NUN; 

YNAX  LINENUH; 

PARENT  :•  SURRBLOB; 

NEXT  :«  NIL; 

IF  COLOR  •  BKGND  THEN  PERIM  ;•  NIL  (•  NO  PERIMETER  LIST  FOB  BACKGROUND  BLOB  •) 
ELSE  BEGIN 
NEVIPERIH); 

PER  IN* .PREV  ;•  PERIN; 

PER IN* .NEXT  ;•  PERIN; 

PERIN*. RIGHT  :•  NIL; 

ADDHIGHTPERINPOINT(PERIN,LINENUN. ST ARTCOL) ,  <«  ADD  FIRST  POINT  •> 

PERIN*. LEFT  :•  PERIN*  .RIGHT: 

FOR  COL  :*  ST ARTCOL* 1  TO  ENDCOL  DO 

ADDRIGHTPERINPOINT(PERIN.LINENUN.COL)  <•  ADD  REMAINING  POINTS  •> 

END  (•  ELSE  •> 

END;  C  WITH  BLOB*  •) 

IP  TRACED I AGNOSTICS  .  ON  THEN  VRITELNI 'DBUG  [NSERTSEONENT  2'); 

(•  FORM  A  NEV  PERIMETER  SECTION  FOR  THE  SURROUNDING  BLOB  •> 

WITH  SURRBLOB*  DO 

IF  COLOR  <  >  BKGND  THEN  BEGIN  _ 

(•  CREATE  A  NEW  PER  I  NET ER  SECTION  AND  INSERT  IT  BEFORE  THE 
CURRENT  PER I NET ER  SECTION  •) 

NEV I NEWPERIN); 

NEWPERIN*. NEXT  :•  PERIM; 

NEWPLRIM* .PREV  :•  PERIN*. PREV; 

PER  IN*. PREV*. NEXT  :«  NEWPERIN; 

PERIN*. PREV  :•  NEWPERIN; 


(•  CONNECT  PREVIOUS  POINT  AND  NEW  POINT  •> 
<•  RESET  LEFT  POINTER  •) 


(•  MOVE  PART  OF  THE  RIGHT  DID  OF  THE  CURRENT  PERIMETER  SECTION  TO  THE 
NEW  PERIMETER  SECTION  •> 


WITH  PERIH*  DO  BEGIN 

NEWPERIN*. RIGHT  :■  RIGHT;  (•  SET  RIGHT  POINTER  OF  MEW  SECTION  •) 

IF  TRACEDI AGNOSTICS  •  ON  THEN  VRITELNI 'DBUG  INSERT  SEGMENT  2A*>; 

WHILE  IRIGHT*. LINE  <>  LINENUH- I )  OR  (RIGHT* .COL  <>  ENDCOL)  _  _  _  .  , 

DO  BEGIN  IF  TRACEDI AGNOSTICS  •  ON  THEN  VRITELNIRIOHT* .LINE, RIGHT* .COL, LINENUM-I .ENDCOL) ; 
RIGHT  :•  RIGHT*. NEXT  END; 

IF  TRACED I AGNOSTICS  •  ON  THEN  VRITELNI 'DBUG  INSERTSBGHENT  2B*>; 

NEWPERIN*. LEFT  :•  RIGHT;  I •  SET  LEFT  POINTER  OF  NEW  SECTION  •) _  _ _ _  _ 

IF  (RIGHT*. NEXT*. LINE  ■  LINENUM-I)  AND  (RIGHT* .NEXT* .COL  >  RIGHT* .COL)  THEN  RIGHT  :•  RIGHT* 
WHILE  (RIGHT* .NEXT*. LINE  •  LINENUM-I)  AND  (RIGHT* .NEXT* .COL  >  RIGHT* .COL)  DO  BEGIN 
IF  TRACEDI AGNOSTICS  •  ON  THEN  VRITELNI RIGHT* .LINE.RIGHT* .COLjLINENUN-l) ; 

RIGHT  :•  RIGHT*. NEXT;  (•  SET  THE  RIGHT  POINTER  OF  THE  CURRENT 

SECTION  TO  THE  ENDPOINT  OF  THE  LINE 
SEGMENT  ON  THE  PREVIOUS  LINE  THAT  IS 
PART  OF  THE  SURROUNDING  BLOB  •) 


r 


IF  RIGHT  •  NEWPERINt . LEFT  THEN  BEGIN 
COL  :•  BIGHT*. COL: 

RIGHT  :■  RIGHT*.  NEXT: 

ADDRIGHTPERINPOIMT (PERIN.LINENUN- I ,OOL> ;  END; 

NEWPER IN*  LEFT* .NEXT  :■  NIL; 

END;  (•  N1TH  PER IN*  •> 

PERIN  :•  NEVPERIN 
END;  <*  WITH  SURRBLOB*  •) 

(•  INSERT  2  SEGNENTS  IN  THE  ACTIVE  LINE  I NNED I ATELY  BEFORE  THE  CURRENT 
SEGNENT  POINTER  •* 

IF  TRACED  I AGNOSTICS  ■  ON  THEN  WRITELN! ’DBUG  INSERTSEGNENT  3’*; 

NEW i A  > ;  NEW ( B  > ; 

At  .STARTCOL  STARTCOL; 

B* .STARTCOL  :•  ENDCOL* I ; 

A* . ENDCOL  :•  ENDCOL; 

Bt.ENDCOL  CURRSEG* . STARTCOL  -  1; 

At. BI.OB  :•  BLOB; 

Bt.BLOB  :•  SURRBLOB; 

PREVSEG* . NEXT  :■  A; 

At. NEXT  ;•  B; 

Bt . NEXT  ;•  CURRSEG ; 

CURRSEO  ;•  A;  (•  RESET  CURRENT  SEGNENT  POINTER  •) 

IF  TRACEACT 1 VEL I NESEGNENTS  ■  ON  THEN  TYPEACT I VEL I NESEOMBfTS ; 

IF  TRACEPERINS  -  ON  THEN  BEGIN 

WRITEPERINSCTTY,  SURRBLOB) ;  WR I TEP  ERINS  (TTY,  BLOB)  END; 

END;  (•  INSERTSEGNENT  •» 

<•«*♦  NEW  PAGE  •) 

PROCEDURE  ADDBLOB  <  VAR  BLOB , TERNBLOB :  BLOBPTR ) ; 

<•  THIS  IS  A  RECURSIVE  PROCEDURE  FOR  ADDING  A  BLOB  TO  THE  END  OP 
THE  LIST  OF  CONPLETED  BLOBS  •) 

BEGIN 

IF  BLOB  •  NIL  THEN  BLOB  ;•  TERNBLOB 

ELSE  ADDBLOB (BLOBt  .NEXT .TERNBLOB) 

END;  <•  ADDBLOB  •> 


PROCEDURE  RECORDBLOB ( VAR  TERNBLOB:  BLOBPTR); 


IF  (HANDONLY  AND  (TERNBLOB* .AREA  <  AREATH) )  OR  BLOBWR ITTEN  THEN  ELSE 
WRITEPERINSCPERINFILE, TERNBLOB) ; 

IF  HANDONLY  AND  (NOT  BLOBWR ITTEN)  AND  (TERNBLOB* .AREA  >•  AREATH)  THEN  BEGIN 
RIjOBWR  ITTEN  ;•  TNUE; 

WRITELN ( 'BLOB  OF  HAND  IS  DONE');  END; 

WITH  TERNBLOB*  .PERIN*  DO  DELETEPER INS(R IGHT .LEFT) ; 

(•  DISPOSE(TERNBLOB*. PERIN);  •)  (»NSP») 

WR I TEBLOB< DATAFILE. TERNBLOB) ; 

ADDBLOB (BLOBSDONE, TERNBLOB) ; 

END;  <•  RECORDBLOB  •) 


(••P»  NEW  PAGE  •) 

PROCEDURE  DELETESEGNENT ; 

VAR  TERNBLOB, R I GHTBLOB . LETTBLOB  REPLACBLOB , ABSORBBIjOB :  BLOBPTR; 

TEHMCONP ,  R1GHTCONP ,  LEFTCONP ,  REPUCCONP ,  ABSORBCONP ;  I NTEGER ; 

RIGHT COLOR:  BINARY; 

NEW AREA:  REAL; 

OLUPERIN:  PTRPER INSECT  ION ; 

PERINLPOINT.PPUNRPOINT,  I HTE1IORPOINT,  NEXT  POINT:  PERINPTR; 

COL, ABSORBRCOL:  INTEGER; 

BEGIN 

IF  TRACED I AGNOSTICS  •  ON  THEN  WRITELN* 'BEGIN  DELETING’); 

TERNBLOB  :*  CURRSEG) .BLOB;  (•  TERMINATING  BLOB  *) 

TKRNCONP  ;•  TERNBLOB*  .COUP; 

RIGHTBLOB  :•  CURRSEG* .NEXT* .BLOB; 

RIGHTCONP  :•  K IGHT Rl. OB*  .COMP ; 

RIGHTCOLOR  :*  RIGHTBLOB* .COLOR; 

LEFTBLOB  :*  PREVSEG* . BLOB: 

LEFTCONP  :*  LEFTBLOB* .COMP; 

IF  RIGHTCOMF  ■  0 
THEN  BEGIN 

REPLACBLOB  :•  RIGHTBLOB;  ABSORBBLOB  :•  LEFTBLOB; 

REPLACCONP  :•  RIGHTCONP;  ABSORBCONP  :•  LEFTCONP  END 

ELSE  BEGIN 

REPLACBLOB  :■  LEFTBLOB;  ABSORBBLOB  :•  RIGHTBLOB; 

REPLACCONP  :•  LEFTCONP:  ABSORBCONP  ;•  RIGHTCONP  END; 

IF  TRACEDIAGNOSTICS  •  ON  THEN  BEGIN 

WRITELN* 'REPLACING  COMPONENT  REPLACCONP , 

'ABSORBED  CONPONENT  •’ .ABSORBCONP) ; 

WRITELN! 'RIGHTCOLOR* ’.RIGHTCOLOR,’  RIGHT  CONP  •  ’.RIGHTCONP)  END; 


(•  CONNECT  THE  ENDS  OF  TVO  PERIMETER  SECTIONS  OF  THE  TERMINATING  BLOB  IP  ITS 
COLOR  IS  NOT  BACKGROUND  •> 

WITH  TERMBLOBt  DO 

IF  COLOR  <  >  BfcGMD  THEN  BEGIN 

(•  ADO  INTERIOR  POINTS  ON  THE  BOTTOM  LINE  OF  THE  TERMINATING  BLOB  TO  THE 
LIST  OF  PERIMETER  POINTS  •* 

FOR  COL  :•  OJRRSEU* .  STARTCOL. 1  TO  CURRSEGt .ENDCOL-t  DO 
ADDLEFTPERINPOINT  < PER  I M , L I NENUN-  1  ,COL>  ; 

<•  IF  THERE  IS  MORE  THAN  ONE  PERIMETER  POINT  IN  THE  BLOB  AND 

IF  THERE  IS  ONLY  ONE  POINT  ON  THE  BOTTOM  LINE  OF  THE  TERMINATING  BLOB, 

IT  IS  AT  BOTH  THE  LEFT  END  CF  THE  CURRENT  PERIMETER  SECTION  AND  THE 
RIGHT  END  OF  THE  NEXT  PERIMETER  SECTION;  DELETE  THE  POINT  AT  THE 
RIGHT  END  OF  THE  NEXT  PERIMETER  SECTION  •> 

IF  (PER IN'  LEFT  <>  PER  IN' . RIGHT)  AND  (CURRSEGt .STARTCOL  •  CURRSEGt .ENDCOL)  THEN 
WITH  PERIN*. NEXT*  DO  RIGHT  :*  RIGHT'. NEXT; 

IF  TRACED  I  AGNOSTICS  •  ON  THEN  WRITELN( 'DISPOSE  CHECK  I  M| 

IF  TRACEFER I MS  •  ON  THEN  WR I TEPER I  MS ( TTY , TERMBLOBt i 

(•  CONNECT  THE  LEFT  END  OF  THE  CURRENT  PERIMETER  SECTION  WITH  THE  RIGHT 
END  OF  THE  NEXT  PERIMETER  SECTION  •> 

WITH  PERIH'  DO  BEGIN 

LEFT*. NEXT  :•  NEXT* .RIGHT; 

NEXT'. RIGHT  :•  RIGHT;  (•  RESET  THE  RIGHT  END  POINTER  OF  THE 
NEXT  PER I NET ER  SECTION  •) 

PREV'.NEXT  :■  NEXT;  (•  CONNECT  THE  PREVIOUS  SECTION  •) 

NEXT t . PREV  :•  PREV  END;  (•  WITH  THE  NEXT  SECTION  •) 

OLDPERIN  :•  PERIH; 

PERIN  :•  PERU*. NEXT;  <•  ADVANCE  THE  PERIMETER  SECTION  POINTER  •) 

(•  IF  THERE  IS  MORE  THAN  ONE  PERIMETER  SECTION,  THEN  DISPOSE  OF  THE  OLD  ONE; 
OTHERWISE,  ADD  THE  BLOB  TO  THE  LIST  OF  BLOBS  DONE  •> 

IF  PERIN  <>  OLDPERIN 

THEN  (•  DISPOSE  (OLDPERIN)  •)  CUSP*) 

ELSE  RECORDBLOB(TERMBLOB) ; 

END;  (•  IF.  WITH  •) 

IF  TRACEDI AGNOSTICS  *  ON  THEN  WRITELN( 'DISPOSE  CHECK  2'); 

IF  TRACEPERIMS  -  ON  THEN  WR  I  TBfER  IMS!  TTY ,  TERM8LOB  >  ; 


(•  CONNECT  THE  ENDS  OF  THE  TWO  PERIMETER  SECTIONS  OF  THE  REPLACING  BLOB  IF  ITS 
COLOR  IS  NOT  BACKGROUND  •> 


WITH  REPUCBLOB*  DO  BEGIN 

IF  COLOR  <  >  BKGND  THEN  BEGIN 

(•  CONNECT  LEFT  END  POINT  OF  ABSORBED  SECTION  TO  A  PERIMETER  POINT  IN  THE 
REPLACING  SECTION  •) 


PER  I NLPO I  NT  :■  PERU* .  RIGHT ;  „ 

WHILE  PER  I  NLPO  I  NT*  .COL  <>  ABSORBBLOB'  .FOUR*.  LEFT'  .COL- 1  DO 
PEKIMLPOINT  ;■  PERIMLPOINT*  . NEXT : 

ABSORBBLOBt. PER  IN'. LEFT'. NEXT  :•  PERIMLPOINT; 

IF  TRACED  I  AGNOSTICS  -  ON  THEN  _ _ 

WR ITELN ( 'PERIMLPOINT  •' .PERIMLPOINT' ,LINE:4, ’ PERIMLPOINT' .C0L:4> ; 
IF  TRACEPERIMS  •  ON  THEN  WRITEPERINSOTY  .REPLACBLOB)  ; 


PER  I MR POINT* .NEXT 


PERIMRPOINT  ;■  PERIN' .RIGHT; 

ABSOR8RCOL  :«  ABSORBBLOB*  .PERIH*  .NEXT*  .RIGHT*  .COL; 

IF  PERIH*  .RIGHT*  .COL  >  ABSORBRCOL*  I  THEN 

WHILE  PERIMRPOINT*  .COL  <)  ABSORBRCOL*  I  DO  PERIMRPOINT 
ELSE  IF  PER  IN*.  RIGHT*.  COL  <  ABSORBRCOL- I  THEN 

FOR  COL  :■  ABSORBRCOL- I  DOWNTO  PERIN*  .RIGHT*  .COL*  I  DO 

ADDR I  GHTI’ER  I MPO I  NT  ( ABSORBBLOB*  .PERIH*.  NEXT ,  L INENDN- 1 ,  COL)  ; 
IF  TRACEDI  AGNOSTICS  •  ON  THEN  WRITELN  ( 'DISPOSE  DBUO  2B’ )  I 
IF  TRACEPERIMS  •  ON  THEN  WR I TEPER IMS (TTY, REPLACBLOB) ; 


(•  DELETE  POINTS  IN  THE  MIDDLE  OF  THE  LINE  SEGMENT  CONNECTING  THE  TVO 
BLOBS  WHICH  ARE  NOT  NOW  PERIMETER  POINTS  *) 


IF  TRACEDI  AGNOSTICS  ■  ON  THEN  _ _ 

WRITELN! 'PERIMRPOINT  ■ ' ■ PERIMRPOINT* .L1NE;4. ' ■ ' .FERIHBPOINT* -C0t:4) ; 

IF  (PERIMLPOINT  <>  PERIMRPOINT I  AND  (PERIMLPOINT  <>  PERIMRPOINT* .NEXT)  THEN  BEGIN 
NEXTPOINT  :*  PEMIMRPOJNT* .NEXT; 

WHILE  NEXTPOINT  <>  PERIMLPOINT  DO  BEGIN 
INTERIORPOINT  ;*  NEXTPOINT; 

NEXTPOINT  :•  NEXTPOINT*  .NEXT f 

(•  DISPOSE! INTERIORPOINT) ;  •)  END  END;  (»NSP«) 

IF  TRACEDI AGNOSTICS  ■  ON  THEN  WR ITELN ( 'DISPOSE  DRUG  4’); 

IF  TRACEPERIMS  •  ON  THEN  WR ITEPER I MS (TTY .REPLACBLOB) ; 


(•  CONNECT  THE  RIGHT  END  POINTS  OF  THE  ABSORBED  AND  REPLACING  SECTIONS  •) 
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IF  PERINLPOINT  •  PERINRPOINT  THEM 

ADbR  IGHTPER I RPO 1  NT  ( ABSORBS  LOB  i  .PERINt .NEXT .LINBIUN, PERINRPOINT t .COL > 

ELSE  BEGIN 

PERINRPOINT*  .NEXT  :•  ABSORBBLOBt . PEKIN'  .NEXT*  RIGHT ; 

ABSORBBLOBt .PERINt .NEXT* .RIGHT  :■  PERINt .RIGHT;  END; 

PERINt  Tight  =-  absorbblob t .per ini .right: 

IF  TRACED 1 AGNOSTICS  •  ON  THEN  URITELNl  'DISPOSE  DRUG  4B’>; 

<•  HEDGE  PER  I NET  EH  SECTION  RINGS  OF  THE  ABSORBED  BLOB  AND  THE  REPLACING  BLOB  •> 

OLDPER1N  :•  ABSORBBLOBt .PERIN; 

IF  (OLDPER 1R  •  OLDPENINt .NEXT)  AND  (OLDPERIN  •  OLDPERINt .PREV) 

THEN  BEGIN  (•  DISPOSElOLDPER IN)  *1  END  (*NSP*> 

ELSE 

IF  REPLACBLOB  <>  ABSORBBLOfi 
THEN  BEGIN 

<•  INSERT  PER  I  NET ER  SECTIONS  OF  ABSORBED  BLOB  BEFORE  PERINt  •) 

IF  TRACED! AGNOSTICS  •  ON  THEN  WRITELNI ’DISPOSE  DBUG  7'); 

OLDPERIN*  .PREVt .NEXT  :■  PERIN; 

OLDPERINt .NEXT) .PREV  ;•  TEKINt .PREV ; 

PERINt. PREVt. NEXT  :•  OLDPENINt .NEXT; 

PERINt. PREV  OLDPER INt . PREV ; 

PERIN  ;•  OLDPER  INt.  NEXT; 

<•  DISPOSElOLDPER  IN)  ;  •)  (•NSP»> 

END 

ELSE  BEGIN 

(•  ADD  PER  I  NET  ER  LIST  TO  THE  BLOB  DESCRIPTOR  OF  THE  TERMINATING  BLOB  (MICH 
HAS  THE  BACKGROUND  COLOR  IF  REPLACING  BLOB  •  ABSORBED  BLOB)  •> 

(•  INSERT  PERI NET ER  SECTIONS  OF  ABSORBED  BLOB  AFTER  PERINt  •) 

IF  TRACED I AGNOSTICS  •  ON  THEN  URITELNl ’DISPOSE  DBUG  8’>; 

PER INt. NEXT t. PREV  :•  OLDPERINt .PREV ; 

OLDPERINt .PREVt .NEXT  :■  PERINt. NEXT; 

PERIN  PERINt. NEXT; 

TERNBLOBt. PERIN  ;«  OLDPERIN: 

OLDPER INt. RIGHT  :•  PER  I APPOINT ; 

OLDPER INt. PREV  :•  OLDPERIN; 

OLDPERINt .NEXT  :»  OLDPERIN; 

REOOHDBLOB I TERNBLOB  >  END; 

IF  REPLACBLOB  <>  ABSORB BLOB  THEN  AREA  :•  AREA  ♦  ABSORBBLOBt .AREA; 

END;  (•  IF  COLOR  <>  BKGND  •) 

END;  (•  WITH  REPLACBLOB!  •) 

IF  ABSORBBLOB  <  >  REPLACBLOB  THEN  BEGIN 

(•  NOTE  —  THE  PARENT  IS  NOT  CONFUTED  CORRECTLY  •> 

IF  TRACEPARENT  THEN  WRITELNI ’LINE  ’ ,LINENUN:3, ’  ABSORBBLOB*’ , 

ABSORBBLOBt ,CONP:3. ’  REPLACBLOB* ’ , REPLACBLOBt .COMP) ; 


(•  CHANGE  ALL  INSTANCES  OF  THE  ABSORBED  COMPONENT  IN  THE  ACTIVE  LINE  TO 
THE  REPLACING  COMPONENT  •) 

SEG  :•  ACTIVELINE; 

WHILE  SEG  O  NIL  DO  BEGIN 

IF  SEGt.BLOB  •  ABSORBBLOB  THEN  SEGt.BLOB  :•  REPLACBLOB; 

WITH  SEGt.BLOBt  DO  IF  PARENT  •  ABSORBBLOB  THEN  BEGIN; 

PARENT  :«  REPLACBLOB; 

IF  TRACEPARENT  THEN  WRITELNI ’REPLACING  PARENT  IN  SEGMENT  FOR  BLOB  *’ .COMP: 4) j  END; 
SEG  :•  SEG t . NEXT  END; 

(•  CHANGE  ALL  INSTANCES  OF  THE  ABSORBED  COMPONENT  NUMBER  IN  THE 
LIST  OF  BLOBS  DONE  TO  THE  REPUCING  COMPONENT  NUMBER  *) 

END;  (•  IF  ABSORBBLOB  <>  REPLACBLOB  •) 

<•  DISPOSE I ABSORBBLOB)  •>  (*NSP») 

IP  TRACED! AGNOSTICS  *  ON  THEN  WRITELNI ’DISPOSE  DBUG  6’); 

<•  DELETE  THE  CURRENT  SEGMENT  AND  THE  SEGMENT  FOLLOWING  THE  CURRENT  ONE  PROM  THE 
ACTIVE  LINE:  SET  THE  CURRENT  SEGMENT  POINTER  TO  THE  FIRST  SEGMENT  AFTER  THE 
DELETED  ONES  •> 

PREVSEGt .NEXT  :•  CTJRRSBGt . NEXT t .NEXT; 

<•  DISPOSE i CURRSEGt .NEXT) ;  •)  <*NSF»> 

(•  DISPOSE  I CURRSEG)  ;  •)  (»NSP»> 

CURRSBG  :•  PREVSEGt . NEXT ; 

IF  TRACEACTIVELINESEGMENTS  •  ON  THEN  TTPEACT I VELINESEONEMTS ; 

IF  TRACEBLOBS  •  ON  THEN  TYPEBLOBSITTY) ; 

IF  TRACEPERINS  •  ON  THEN  TYPEPERIMSITTY) 

BID;  (•  DELETESEOMENT  •» 


I 


(••P*  NEV  PAGE  •* 

PROCEDURE  UPDATE? ERIMtVAR  PERIN:PTRPERINSECTION;  CURRST ARTCOL , CURRENDCOL , 
NEWSTART COL ,  SEVEN DCOL :  INTEGER) ; 

(•  FOR  CASE  3  ONLY  •) 

VAR  COL:  INTEGER; 

BEGIN 

(•  LEFT  SIDE  OF  PERIMETER  SECTION  •> 

IF  CURRST ARTCOL*  I  <  NEVST ARTCOL  THEN 

FOR  COL  CURRST  ARTCOL*  I  TO  N  EVST  ARTCOL  - 1  DO 
ADDLEFTFERIMPOtNTIPERIM.LINENUM- I ,COL> ; 

IF  CURRST  ARTCOL- 1  >  NEVST ARTCOL  THEN 

FOR  COL  :■  CURRST ARTCOL- I  DOWNTO  NEVST ARTCOL* I  DO 
ADDIEFTPER I MPO 1 NT  1  PER  I N . L I NENUN . COL ) ; 

ADOLEFTPER INPOINT IFERIM.LINENUM.NEVST ARTCOL) ; 

PERIN  :•  PER  IN* .NEXT;  <•  ADVANCE  PERIMETER  SECTION  POINTER  •> 

(•  RIGHT  SIDE  •) 

IF  CURRENDCOL  >  NEVENDCOL  *  I  THEN 

FOR  COI.  :•  CURRENDCOL- 1  DOVNTO  NEVENDCOL*  l  DO 
ADDKIGNTrER INPOINT (PER  IN, LI NENUN- 1  ,COL)  ; 

IF  CURRENDCOL  <  NEVENDCOL  -  I  THEN 

FOR  COL  :•  CURRENDCOL*  I  TO  NEVENDCOL- 1  DO 
ADORIGHTPER  INPOINT  (PERIN, L1NENUN, COL)  ; 

ADDR I GHTPER  INPOINT  (PERIN ,  LI  NENUN,  NEVENDCOL) ; 

END,  <•  UPDATEPERIN  •> 


(••P*  NEV  PAGE  •) 

PROCEDURE  CASI3PROCESSING(VAR  CURRSEG : SEGPTR ,  ST ARTCOL, ENDCOL: INTEGER) ; 

VAR  NEVAREA.XLINEAVE:  REAL; 

BEGIN 

IP  TRACED I AGNOSTICS  •  ON  THEN  VRITELNI 'BEGIN  CASE  3  PROCESSING’); 

VITH  CURRSEG'  .  BLOB*  DO 

IF  CONP  <>  6  THEN  AREA  AREA  ♦  ENDCOL  -  ST ART COL  ♦  1; 

IF  (CURRSEG'. BLOB*. COLOR  O  BKGND)  AND  (CURRSEG* -BLOB' .TWIN  <>  LI NENUN)  THEN _ 

UPDATEPERIN  (CURRSEG*  .BLOB*  .  PERIN,  CURRSEG*  .ST  ARTCOL,  CURRSEG*  .ENDCOL.  ST  ARTCOL.  ENDCOL) ; 

IF  TRACEBLOBS  •  (Ml  THEN  WR I TEBLOB( TTY .CURRSEG* .BLOB) ; 

IF  TRACEPEKINS  *  ON  THEN  HR  ITEPER  INS  (TTY  .CURRSEG*  .BLOB) 

END;  (•  CASE3PROCESS I NG  •) 

(••P*  NEV  PAGE  •> 

PROCEDURE  PROCESSLINE; 

VAR  ONKZERO:  INTEGER; 

BEGIN 

IF  LINENUN  >  1  THEN  BEGIN  _ 

NEV ( BLOB > ;  (•  CREATE  A  BLOB  DESCRIPTOR  FOR  THE  BACKGROUND  •) 

WITH  BLOB*  DO  BEGIN 
COLOR  :•  BKGND. 

CONP  0; 

PERIN  :•  NIL; 

AREA  :*  0.0; 

XNEAN  :•  0.0; 

YNF.AN  0.0; 

XNIN  :«  0; 

XNAX  :>  0; 

YNIN  :•  0; 

YNAX  0; 

PARENT  :•  NIL; 

NEXT  :•  NIL  END; 

(•  INITIALIZE  THE  ACTIVE  LINE  •) 

NEV (ACTIVELINE) ;  (•  THE  FIRST  SEGMENT  CORRESPONDS  TO  THE  BACKGROUND 

WHICH  IS  TO  THE  LEFT  OF  THE  FIRST  COLUMN  •) 

ACTIVELINE* .STARTCOL  :■  NINIHT*1; 

ACTIVELINE*. ENDCOL  :•  NAXINT-2; 

ACTIVELINE*  .BLOB  :«  BLOB:  _  _  _ _ _  __  „„  . 

NEV (CURRSEG) ;  <•  THIS  SEGMENT  DESCRIPTOR  CORRESPONDS  TO  THE  BACKGROUND 

WHICH  IS  TO  THE  EIGHT  OF  THE  LAST  COLUMN  •) 

ACTIVELINE*. NEXT  :•  CURRSEG; 

CURRSEG* .STARTCOL  :•  NAX1NT-I; 

CURRSEG) .ENDCOL  :•  NAX1NT-I; 

CURRSEG*  BLOB  ;•  BLOB; 

CURRSEG* . NEXT  :■  NIL 
END;  (•  IF  LINENUN  •  I  •) 


-V 


IF  LINENUN  >  KSKIF  THEN  BEGIN 


CURRSBG  :•  ACTIVELINE;  <•  INITIALIZE  THE  CURRENT  SEGMENT  POINTER  •> 

NEWSEG  ;•  NEWLINE;  (•  INITIALIZE  THE  NEW  SEGMENT  POINTER  •> 

REPEAT 

WITH  NEWSEGt  DO  BEGIN 

IF  TRACED I AGNOSTICS  •  ON  THEN  BEGIN 

VRITELNI  ‘ClIRRSEG*  .STARTCOL  •  '  .CURRSEG*  .STARTCOL, 1  CURRSBG t  .ENDCOL  •  ’CURRSEG* .ENDCOL) ; 
WR 1TELN ( ’ NEWSEGt . STARTCOL  •  ’.STARTCOL,'  NEWSEGt . ENDCOL  •  ’.EKDCOU  END; 

IF  CURRSEGt .BLOB* .COLOR  •  INK  THEN  ONEZERO  ;•  I 

ELSE  ONEZERO  :•  6: 

WHILE  STARTCOL  >  CURRSEGt .ENDCOL* ONEZERO  DO  DELETESEGMENT ;  (•  CASE  1  •> 

IF  TRACED I AGNOSTICS  •  ON  THEN  WRITELN ( ' DBUG  PROCESSLINE  i!>; 

IF  ENDCOL  <  CURRSEGt . STARTCOL -ONEZERO  THEN  I NSERTSEGMENT (STARTCOL, ENDCOL) ;  (•  CASE  2  •) 
CASK3PROCESS I NG  ( CURRSEG ,  STARTCOL ,  ENDCOL  >  ; 

CURRSEGt  .STARTCOL  :•  STARTCOL; 

CURRSEGt .  ENDCOL  ;•  ENDCOL; 

END;  (•  WITH  NEWSEGt  •> 

PREVSEG  :*  CURRSEG;  <•  ADVANCE  THE  PREVIOUS  SEGMENT  POINTER  •) 

CURRSEG  :•  CURRSEGt . NEXT ;  <•  ADVANCE  THE  CURRENT  SEGMENT  POINTER  •) 

NEWSEG  NEWSEGt .NEXT  (•  ADVANCE  THE  NEW  SEGMENT  POINTER  *> 

UNTIL  NEWSEG  >  LASTNEWSEGt . NEXT ; 

IF  TRACEACTIVEL1NESEGNENTS  •  ON  THEN  TYPEACT1VELINESEGNENTS; 

IF  CURRSEG  «  NIL  THEN  BEGIN  WRITELN (TTY, 'ERROR:  CURRSEG-NIL1 > (PAUSE  END; 

WHILE  CURRSEGt. NEXT  <>  NIL  DO  DELETESEGMENT;  (•  PROCESS  ALL  SEGMENTS  REMAINING 

IN  THE  ACTIVE  LINE  EXCEPT  THE 
LAST  SEGMENT  •> 


END;  <•  IF  LINENUN  >  NSXIP  •) 
EMD;  (•  PROCESSLINE  •> 


NEW  PAGE  •» 

BEGIN 

MANDONLY 


(•  MAIN  PROGRAM  •> 


;•  TRUE;  BUWWKITTEN  :•  FALSE; 
WRITELNCTTFE  INPUT  IMAGE  NAME:'); 
KEAOUf ;  READ!  INPUTXAME)  ; 

RESET!  INPUT  IMAGE.  INPUTNANE)  ; 


NEWS  I TEC DATAFILE, 'DUMMY. DAT'): 
-  OUTPUT  PERIMETER 


FILE  NAME;'); 


(•NSP*) 

(•NSP»> 

!»NSP»> 


WR  (TELN  l  'TYPE _  _ 

READLN;  READ (PERI MM AME ) ;  (•MSP*) 

REWRITECPKRIMFILE.PERIMNAME) :  C*NST*) 

UR  I TELN C 'TYPE  THRESHOLD  FOR  BLOB  AREA:'); 

READLN:  READ! IAREATH) ;  AREATH  :•  IAREATH; 

WR ITELN ( 'PERIMETER  LIST  COMPACTED  EVERY  X  LINES;  TYPE  X:’>! 

READLN;  READ(LINEGROUP) ; 

TRACEPARENT  :■  FALSE: 

WR ITELN ( 'TYPE  *  OF  LINES  TO  BE  SKIPPED  AT  THE  TOP  OF  THE  IMAGE:'); 
READLN;  READ(NSXIP) ; 

RECURS  :•  FALSE; 

WRITELN!  'TYPE  "I"  FOR  INTERACTIVE  DEBUGGING’) ; 

WR  ITELN  ('  OR  "U"  FOR  UNINTERRUPTED  PROCESSING:’); 

READLN;  READ(CND) ; 

IF  CMD  •  ’ I '  THEN  INTERACTIVE  :»  TRUE 
ELSE  INTERACTIVE  :•  FALSE; 

VRITELNI  TYPE  ”T”  TO  TYPE  OUT  LINE  NUMBERS:’); 

READLN;  READ(ONOFF) ; 

IF  ONOFF  ■  'T'  THEN  TYPELINENUMS  :•  TRUE; 

LINENUN  :■  0;  NEVCONPNUN  :•  0;  BLOBSDONE  :•  NIL;  RECYCLEDFTR  :•  NIL; 
CHARTOT  0; 

NEWCARRAYP  III): 

FOR  I  :-  I  TO  10  DO  BEGIN 
NEWCARRAYP  11*11); 

ARRAYP lilt. NEXT  :«  ARRAY? 1 1 • 1 1 1  END; 

ARRAYF (lilt. NEXT  :■  NIL; 

IF  NOT  INTERACTIVE  THEN  NLREAD  :•  I; 

TRACEBLOBS  :«  OFF;  TRACEFEXINS  :•  OFF;  TRACERUNLENGTHS  c»  OFF; 
TRACEACTJVELINESEGMENTS  :•  OFF;  TRACED1 AGNOST ICS  :•  OFF; 

CMDSET  :•  I’L’,  'B',  ’P\  'A',  'R' ,  'D' .  'S’l; 


REPEAT 

IF  INTERACTIVE  THEN  BEGIN 

WRITELN ( 'TYPE  OOMNAND:  ”L”  ,  ”B”.  "F” ,  "R" .  "A” ,  ”D” , 
READLN:  READ(CMD) ; 

IF  (CMD  •  LINEIN)  THEN  BEGIN 

WRITELNCTTFE  NUMBER  OF  LINES  TO  BE  READ:  ’ )  I 
READLN;  READ (NLREAD)  DID 
ELSE  NLREAD  :•  I  END 


OR  "S”'); 


IF  CRD  •  LINE  IN  THEN  FOR  NT  I  RES  :•  t  TO  NLREAD  DO  BEGIN 
LINENUM  :•  LINENUN  ♦  I; 

IF  TYPELIHENUMS  THEN  WRITELNI 'LINE  ' ,LINENUM:3> I 
GETLINE; 

IF  (LINENUN  •  II  OR  (LINENUM  >  NSKIP)  THEN  BEGIN 

IF  LINENUM  NOD  LINEGROUP  •  0  THEN  PACKPERIRSt 
SETLINE, 

CREATERUNLENOTRSi 
PROCESSLINE ; 

END) 

END 

ELSE  IF  I  CUD  <>  STOP)  AND  (CND  IN  CMDSET)  THEN  BEGIN 
READ(ONOFF) ; 

CASE  CND  OF 

DIAGNOSTICS:  TRACED 1 AGNOSTICS  !•  ONOFFi 
LIST  BLOBS:  TRACEBLOBS  :•  ONOFF; 

LIST  I*  ERINS:  TRACEPER  IMS  :*  ONOFF; 

LISTRUNLENCTHS:  TRACERUNLENGTHS  :•  ONOFF; 
LISTACriVELINESEGMENTS:  TRACEACTIVELINESECMENTS  :•  ONOFF 
END;  (•  CASE  •) 

END;  (•  IF  •) 

UNTIL  (CND  •  STOP)  OR  (LINENUN  •  NL1NES)  OR  BL0BVR1TTEN; 

LINENUN  :•  LINENUN  ♦  1; 

SET1.INE;  (•  PROCESS  A  LINE  WITH  NO  RUN  LENGTHS  TO  •) 

PROCESSLINE;  <•  ENSURE  COMPLETION  OF  BLOB  PROCESSING  •> 

WRITELNI 'LINEBYLIN  —  ALL  DONE'); 


PROGRAM  SCAN 

C  -  Proirn  to  raad  laromltil  aaqaaatial  data  filag  of  hul  data 
C  -  that  war#  gaaaratad  with  prog  LINEBYLIN  ea  E1KONIX  gaaag, 

C  -  aid  attract  position  ii  ipaoa  of  IKa  fiigors  aid  tKa  oorraspoadiig 
C  “  iatardigital  spaoas. 

1  FORMAT  ( Id, IX, 14) 

2  FORMAT  </'  Tipi  is  filaaaaa  of  iapat  data  fila  »> 

4  FORMAT  (Q.34A1) 

5  FORMAT  (Id) 

6  FORMAT  ('  FINGER  •  INCLUDES  INTERDIGITAL  SPACES  >  I  a  LITTLE 
♦  FINGER') 

7  FORMAT  (/IX. 'FINGER  *\3X,’ARRAY  POSITION'  3X. 'XVAL' ,SX, 'TTAL' . 
♦SX, 'LENOTH’.dX, 'WIDTH' ,2X, 'RATIO  lyW’ ,3X, 'SLOPE' .dX. 'dIst’ I 

8  FORMAT  (2X, 12, 10X, 15,8X,F8.2,X,F8.2,X,F8.2, 
♦X,F8.2.X,F8.2.X.F8.2.2X,F8.3) 

9  FORMAT  (/•'  Laagth  of  leaa  wi.dow  (12*6aa>  •  ') 

FORMAT  ( 12) 

FORMAT  (/IX,'  Laagth  of  goat  aiado*  •  ’,14,’  (12  •  8  aa>') 
FORMAT  (/»'  Aaglo  for  flagon  ■  ’) 

FORMAT  (F6.0) 

FORMAT  (/•'  Aaglo  Tor  thaab  •  ') 

FORMAT  (/•'  Fingar  aaglo  •  *,F6.8,'  Thaab  aagla  •  '.F6.0I 
FORMAT  (/////IX, 'Iapat  fila  1 , I4A1.2X, 'Ootpal  fill  '.I4A1) 
FORMAT  (2X.I6.6(2X,F8.2>> 

FORMAT  (/’  At  Thaab  Crotoh:  Handwidth  •  ’ ,FB.2,2X, ’Paria  • 

♦F13. 2,2X,  'Aral  •  •  .FI3.2.2X,  'Ratio  (P**2/A>  ■  '  ,H.3) 

FORMAT  (/•’  Tjpi  1  to  haaa  data  priatad  oat  at  aad  of  toft  ’> 
FORMAT  (/•'  Tap#  1  to  aark  aoaaara  poiata  oa  pala  plot  ’> 

FORMAT  (7(X, Is) ) 

FORMAT  (/’  Entire  Hoad:  Faria  •  *  .FI3.2.2X, 'Ana  •  ’.F13.2.2X, 
♦’Ratio  (P**2/A)  ■  '.F6.3) 

FORMAT  ('  Type  I  if  yoa  aaat  to  ehaago  paraaotora  ’> 

FORMAT  (3l3,X,7FI0.3,/,( I0X.7F10.3) ) 

FORMAT  ('  Tip*  aaaorio  eodoa  for  SUBJECT  ID,  SESSION,  TRIAL  ’> 
FORMAT  (14,14,14) 

FORMAT  (/'  Typa  ia  filaaaaa  for  oatpat  data  fila  ’) 

FORMAT  (/'  Typo  DIST  aad  R  critarioa  for  oiraatara  aaal’) 
FORMAT  (F8.2.FB.3) 

FORMAT  (/'  Typo  aaabar  (odd)  of  poiata  to  gkippod  at  atari') 

DIMENSION  END(2> ,NID(2> .HEAD (2) .TEMP (2) 

DIMENSION  FNGT IP ( 18,2) jFNGLENl 18) ,FNGVID(I8) ,VALUE(30) 

DIMENSION  SLOPE! 10) ,ZYER0< 10) ,LINE(30) 

DIMENSION  TANGLE( 100) ,DLENG( 100) , CURVE) 10) 

INTEGER  IDATA(6000) , FINGER, FNGNUR( 18) , IVIDTH( 18) , ID!N(6O00) 
LOGICAL*)  IFILE(34),IFILE2(34) 


ANGLE)  •  90. 
ANGLE2  •  90. 
IDIS  •  18 
DISTC  •  40.0 
CHIT  •  0.90 
I  START  •  201 
I FLAG  •  0 
IF  •  0 
TYPE  2 
ACCEPT  4, 
TYPE  28 
ACCEPT  4. 


IDafaall  eritarioa  aagla  ahift  for  fiagtipa 
tDofaalt  oritarioa  aagla  ahift  for  thaab 
tDofaalt  dial  botvaoa  poiata 
'Storting  dial  f roa  fiagtip  for  oarao 
ICritoriof  fit  for  oaraotaro  fa. 
tSkip  thaaa  firat  poiata 


151 ,  ( IFILE( J) ,  J«t,ISI> 

152,  ( IFILE2( J) ,  J-1.1S2) 


A-19 


TYPE  26 

ACCEPT  27,  lAUTH, ISESS, I TRIAL 
TYPE  20 

ACCEPT  S,  I FLAG 
TYPE  21 

ACCEPT  5,  HPLOTF 

TYPE  24 

ACCEPT  5,  IF 

IF  (IF  .EQ.  0>  GO  TO  49 

TYPE  9 

ACCEPT  5,  ID1S 
IDIS  -  ID1S/2 
TYPE  12 

ACCEPT  13,  ANGLE I 
TYPE  14 

ACCEPT  13.  ANGLE2 
TYPE  29 

ACCEPT  30,  DISTC.CNIT 
TYPE  31 

ACCEPT  S.  I START 
CALL  ASSIGN! 1 , IFILE, IS1) 

DO  50  £* l  32000  2 

READ  d’l, END-51. ER*«5I>  IPIN(K) ,  IDIN(K*1> 
STOP  •?  NO  EOF  7  1 
IF  (IDIN<K>  .ME.  0)  00  TO  52 
X  •  K  -  2  IDisreftrd  otorraa  poiats 

00  TO  SI 
NPTS  -  X 
CALL  CLOSE!  I) 

CALL  SORT! IDATA.IOIN, NPTS* 1) 


9.  '  4*f  la  altar j tha 

C  -  Iaitialite  cooators  -  fill  eatl*  at  kata  throstH  data,  eoeateetljr 
C  -  apdatint  tha  average  aa|ta.  Moira  throifh  data  with  a  wiedow 
C  -  of  tixe  IDIS,  in  slept  of  IDIS/2.  Paiats  etefd.  with  triaatla  fa. 

C  -  Vhaa  fiad  a  large  shift  is  tha  average  aagla  (ANGLE  for  fitters, 

C  -  ANGLE2  f  or  thaab)  store  looatioa  of  tip,  raiaitialixa 
C  -  average  eatte.  Store  tha  looatioa  of  tha  tip  ia  tha 
C  -  data  arrajr,  as  vail  as  X  asd  Y  positioas  ia  tpaaa. 

C  -  Calcalate  froa  these  both  laattn  aad  width  of  fiatart 
END!21  •  AVE!  1DATA,  1ST  ANT*  1 ) 

END! 1 )  •  AVE! IDATA, 1ST ART) 

NID(2I  •  AVE! IDATA, IDIS* 1START* I > 

NIDI  I >  *  AVE! IDATA. IDIS*ISTART> 

HEAD (2)  •  AVE(lDATA.2*IDIS*ISTAirr*l) 

HEAD! 1 >  •  AVE! IDATA, 2*I0IS*ISTART> 

ANGLE  -  ANGLE  I 
FINGER  -  I 
1CNT  -  0 
N  •  I 

ALPHA  «  F ALPHA (END, HEAD I 
SVNALP  •  ALPHA 
AVEALP  •  ALPHA 

DO  208  lOaiDISMSTART, NPTS,  IDIS 
END! I)  •  MID! I > 

END(2>  •  HID(2) 

NID(I)  •  HEAD ( 1 ) 

N1D(2I  •  HEADI2) 

HEADI2)  •  AVE! IDATA, 1*1) 

HEAD! I >  •  AVE! IDATA, I > 

IF  (I  .LT.  N-l)  00  TO  200 
ALPHA  •  FALPHA! END, HEAD) 

IF  (ABS! ALPHA- AVEALP)  .LE.  100.)  00  TO  79 
IF  (AVEALP  .or.  0.0)  ALPHA  •  ALPHA  •  360. 

IF  (AVEALP  .LT.  0.0)  ALPHA  ■  ALPHA  -  360. 

79  DALPHA  »  ABS I ALPHA  -  AVEALP) 

C  IF  (FINOER  .LE.  4)  PRINT  17. I ,EH0(2> ,END( l) .HEAD (2) .HEAD! I ) . 

C  ♦ ALPHA, AVEALP 

IF  ( ICNT  .LT.  6)  00  TO  99  tNove  at  least  lOsa  before  teat 
IF  (DALPHA  ,OT.  ANGLE)  00  TO  100 
99  N  •  N  *  1 

ICNT  •  ICNT  ♦  1 
SVNALP  •  SVNALP  *  ALPHA 
AVEALP  •  SVNALP /N 
00  TO  200 

-  foaad  a  fitter 

199  N  •  I 

ICNT  >  a 
SVNALP  •  ALPHA 
AVEALP  «  ALPHA 

IF  (FINGER  .EO.  9)  AINU  ■  AN0LE2 
Ilf  IF  (FINOS  OT.  7)  00  Tp  165 
C  *  Far  fitters  fist  test  ths  loaatios  af  tha  tip 

C  *  ti  dratri at  •  !•••  throsth  tha  aaatar  af  the  fiatar  aad  fled  tha  poiat 
C  -  or  itoraootiaa  with  tho  tip. 

L  •  IDELTS! IDATA, I— IDIS— 1 ,52. ,-2) 

N  •  IDELTS! IDATA, I -IDIS* l,S2. ,2) 
t  ■  IDELTS! IDATA, l-IDIS-l, 72., -2) 

X  •  IDELTS! IDATA, l-IDIS-1, 72.. 2) 


LINE(F1NGKR*4-3)«J-1 
LINE(FINGFJl*4-2)  »L- 1 
LINE(FINGER*4-I )  «N- 1 
LIN£(F1NGER*4>  -K-l 
XAVEI  •  AVEWDATA.L) 

XAVE2  •  AVE( I DATA, Ml 
YAVF.I  •  AVE( IDATA.L-l) 

VAVE2  •  AVE( 1DATA.N- I > 

XLM  •  (XAVEI  ♦  XAVE2)/2 
YLN  ■  (YAVE1  ♦  YAVE2)/2 
XAVEI  •  AVE< I DATA, J) 

XAVE2  •  AVF.(IDATA.K) 

YAVEI  •  AVE< (DATA. 2-0 
YAVE2  ■  AVE(IDATA.K-l) 

XJK  •  (XAVEI  ♦  XAVE21/2 
YJK  -  (YAVEI  ♦  YAVE2>/2 
SLOPE (FINGER)  -  (YUI  -  YJK>/(XLM  -  XJK) 

Y2ER0IFINGER)  •  YJK  -  SLOPE(FINGER)*XJK 

0  •  2. 

LI  -  L 
LI  ■  II  ♦  2 

APR  I  ME  »  FLOAT  ( IDATAILl- 1 ) )  -  SLOFEIFINGER)  •  FLOAT  UDATAILl  >  > 
IF  (ABSlYZERO (FINGER)  -  APR I HE)  .LT.  D)  GO  TO  166 
IF  (LI  .LT.  M>  GO  TO  ISO 
LI  •  L 
D  ■  2*D 
GO  TO  ISO 
TYPE  13.  D 

FNGNl/Ml  FINGER)  •  LI- 1 
FNCTIP(FINGER,2>  •  IDATAILl) 

FNGT IP (FINGER, 1)  •  IDATA(Ll-l) 

GO  TO  190 

FNGNUN(FINCER)  •  1-IDIS 
FNGTIP(FINGER,2)  •  MID(2) 

FNGT IP (FINGER, I )  «  MIDI 1> 

FINGER  ■  FINGER  ♦  1 

IF  (FINGER  .EX).  18)  GO  TO  281  fFoaad  all  fiogsrs 
CONTINUE 
IDIS  •  IDIS*2 

Calealats  leagth  of  fiager 
DO  S18  1m  ,9,2 
J  -  I 

IF  (1  .GT.  6)  J  •  -I 

DELTAX  •  FNGTIPI I ,2)  -  FNGTIPI I*J,2I 

DELTAY  *  FNGTIPI  1,1)  -  FNGTIPI I«J,  1>' 

FNGLEN(I)  •  SORT (DELTAX* *2  ♦  DELTaY**2) 

FNGLENl  I )  •  DSTHCEI I  DATA ,  FNGNUM  ( I  > »  1  ,FNGND)f  (I*J>*1) 

CONTINUE 
TYPE  18.1 

Caloslats  aidlh  af  fiagers  at  a  dlataaoa  DIST  back  fra*  Up. 

DIST  -  92. 

DO  521  I  •  1,9,2 

L  •  IDELTSI IDATA.FNGNUMI I)* 1 ,D1ST ,-2) 

M  •  IDELTSI IDATA.FNGNUMI I )♦!, DIST, 2) 

DELTAX  •  FLOAT! IDATAIL)  -  IDATAIR) ) 

DELTAY  -  FLOAT  1 1  DATA  I L- 1 )  -  IDATa(R-D) 

FNGVIDII)  -  S0RT(DELTAX**2  ♦  DELTAY**2) 

FNGWIDII)  •  DSTNCEI IDATA.L.R) 

DIST  •  123. 

IF  (I  .EO.  7)  DIST  •  92. 

IF  (I  .EO.  9)  1  «  10 
IWIDTHI I»2)  -  L-l 
IW1DTH(1»3)  •  M-l 
CONTINUE 
TYPE  18.2 

Calaalata  earaatara  of  Uagsrtips  fro*  ♦  to  -  DIST  from  tip. 

First  Had  dal  la  angla  for  taafaatg  arooad  tho  tips,  that  fit 
tha  aagla  function  with  a  straight  lino  aad  gtora  slops  sad  fit. 


DIST  •  DISTC 
DO  530  I* 1, 9,2 

CALL  TANANGl IDATA,TANGLE,DLENG,D1ST,FNGNURII)» 1 ,6, 12, NP) 
CALL  LINFIT(DLENG,TANGLE,NF,SL0PE2,R) 

IF  (ABS(R)  .GT.  CRIT)  GO  TO  529 

DIST  •  DIST  -  2.0 

IF  (DIST  .LT.  4.8)  GO  TO  529 

GO  TO  525 

CURVE! I)  •  SL0PE2 

CURVE) 1*1)  »  DIST 

DIST  •  DISTC 

TYPE  531.  R 

FORMAT  (re. 3) 

CONTINUE 
TYPE  10,3 


C 


C  -  Caloalata  Hendvidth,  pariaatar  aid  arc*  9A 1 1  oa  ia  froa  tha  tipa  af 
C  -  tha  liltla  and  iadei  fiaaers,  raapao  t  i  *al  jr . 

L  *  IDELTS(IDATA,FNGNUR( I) ♦ I , 185. ,-2> 

M  •  . OELTS< IDATA ,FNGNUM(7) ♦ 1 ,226. ,2) 

DELTAK  •  FLOAT ( IDATA (L>  -  IDATA(M) ) 

DELTAY  ■  FLOAT ( IDATA (L-l)  -  1 DATA (A- 1 ) ) 

HNDWID  ■  SORT (DELTAX-*2  ♦  DELTAY—2) 

HNDVID  •  DSTNCEl IDATA.L.M) 

CALL  ARFRIM! IDATA.L.M, AREA, PERIM) 

FINGAR  •  AREA 
FNGPRM  >  PERIN 
IWIDTHI2)  •  L-l 
IVIDTHt 1 1 )  •  M-l 
TYPE  10. 4 


C  - 
C  - 
C  - 


C  - 

c  - 


C  -  Fiad  araa  aad  pariaatar  —  taraiaata  haad  fraa  paiat  I0aa 
C  •  praaiaal  to  tha  tip  of  tha  thaab  ta  a  poiat  I3aa  praaiaal  to  tha 
C  -  tip  of  the  little  fineer. 

L«  IDELTSt I  DAT  A ,  FNGNUM  < 1) »1 ,267. ,-2> 

M  •  IDELTSt IDATA, FNGNUN!9) *  1 ,20S. ,21 
IVIDTHt l>  •  L-l 
IVIDTHt  Ml  •  M-l 

CALL  AJIPR  IN  (IDATA.L.M,  AREA.  PER  IM> 

TYPE  10.5 
C 

C  -  Priat  iafo  aad  write  faatara  aalaaa  ta  fita.  Faatiroa  writtaa  ia  ardart 
C  -  loaeth,  width,  ratie<L/V> .haadwidth.f iapar  pariaatar  aad  araa, 

C  -  ratio  (fiaaar  A/PI,  haad  pariaatar  aad  area,  ratio  (A/P>. 

RATIO  •  !(FNGPRM/20.S)-*2)/(FINGAR/420.  > 

RAT 102  •  ( (PERIM/20.S)**2. >/< AREA/420. > 

IF  tIFLAG  .NE.  I)  GO  TO  551 
PRINT  22, t IVIDTHt I ) ,  1-1.71 
PRINT  22,(IVIUTH(II,  1*8, 141- 

PRINT  16, t IFILEt 1 1 ,  I* 1 , 14) , ( 1FILE2C I > ,  1-1,14) 

PRINT  IS.  ANGLE!, ANGLE2 
PRINT  11,  1DIS 
PRINT  6 
PRINT  7 

DO  550  I- I, FINGER- 1,2 

PRINT  8,  ( I ♦ 1 >/2,FNGNUM( I) , FNGTIPt  1 ,2)  .PWTIPtl ,  1) , 

♦FNGLEN 1 1 ) /20 . S , FNOV t D ( 1 1 /20 . 5 , FNGLEN ( I ) /FNOV I D ( I ) , 

-CURVE! I). CURVE! l-l) 

PRINT  19,  HNDVID/20.S,FNGPRN/20.S,FINGAR/420.,RATIO 
PRINT  23.  PER  1 R/20 . 5 , AREA/420. , RAT 1 02 
CALL  ASS >GN  (1 , 1FILE2, IS2) 

DO  G00  l-l, S 

VALUE!  I)  ■  FNGLEN!  I-2-IV20.5 
VALUE! I-S)  •  FNGVID! 1-2-1 1/20.5 
VALUE! I ♦ 10)  ■  FNGLEN! 1-2- 1 ) /FNGVID ( 1*2-1 ) 

CONTINUE 

HNDVID/20.5 
FNOPRN/20.5 
FINGAR/420. 

RATIO 
PER I  A/ 20. 5 
AREA/420. 

RAT  102 
1.9,2 


550 


SSI 


VALUE!  16) 
VALUE!  17) 
VALUE!  18) 
VALUE!  19) 
VALUE! 17) 
VALUE!  18) 
VALUE!  19) 
DO  601 


CURVE! I ) 
CURVE! 1-1) 


S01 


VALUE!  I- 19) 

VALUE!  1-20) 

CONTINUE 
NUMFEA  •  29 

URITE  (1,25)  (IAUTH, ISESS, ITRIAL, (VALUE! J) ,  J-l.NUMPEA)) 


C  -  Plot  haad  pariaatar  aaaberiaa  peiata  far  fiapar  tipa  aad  iatardlplal 
C  -  apaoas. 

1 YOFST  •  600 
IXOFST  •  100 
J  -  t 

CALL  INITT<96@) 

CALL  MOVABSI IXOFST- 1 DATA (2) , I YOPST- IDATA! I >  > 

DO  700  1-3.NPTS.2 

CALL  DRVABS! IXOFST* (DATA! 1*1), 1 YOFST- IDATA! I > > 

IF  (NFLOTF  .NE.  1)  00  TO  710 
IF  (I  .NE.  IVlimili))  00  TO  700 
GO  TO  711 

IF  ! I  .NE.  LINE(J>)  00  TO  700 
IF  (I  .NE.  FNONUM ( J) )  00  TO  700 
ENCODE! 2, 10, JJ) J 
CALL  A0UTST(2, JJ> 

J  -  J  ♦  1 
700  CONTINUE 

C  -  Plot  lieaa  that  iataraapt  fiaaar  tipa 
00  TO  1002 

DO  1000  X- 1,7  _  _ 

CALL  NOV ABS( IXOFST, I  YOFST*TZEMO(X>> 

DO  1001  1-100,900,100 

IY  •  IFIX(YzERO(X>  *  SLOPE! K)«FLOAT( 1) ) 

CALL  DRVABS! IXOFST* I, I YOFST* IY) 

CONTINUE 
CONTINUE 

CALL  FIN1TT(0,79O> 

STOP 
DID 


C  - 

710 

711 


1000 

1002 


A-22 


nn 


r- 

i 


c 

c 


FUNCTION  AVE(IDATA.I) 

Faaatioa  ie  do  •  fi»a  poiat  lrtti|itir  mri|«  u4  ratara  itl«  u 
AVE. 

DIMENSION  I DATA (6000) 

FAMT1  •  FLOAT! IDATA! 1-4) )•«. 1 
PART2  •  FLOAT ( I DATA ( 1-2) ) *0.2 
PART3  •  FLOAT ( I DATA ( I ) >*0.4 
PART A  •  FLOAT ( I DATA ( I#2>>»0.2 
PANTS  •  FLOAT! I DATA ( 1*4) >*0. 1 
AVE  •  PANT  I  ♦PART2*PART3#PART4#PART5 
RETURN 
END 


FUNCTION  IXYGETdX,  7,  I  DATA) 

C  -  Faaatioa  to  fi*4  i  poiat  with  n  X<V>  oaiiiotaat  to  aaothar  poiit  (IX) 
C  -  J  il  tho  itirtiii  poiit  for  tho  nareh  thraagh  tho  array. 

C  -  If  J>IO  tho  narch  wilt  bo  backward!  throigh  tho  orroy. 

C  -  Rotoroo  (I)  for  tho  oo«  poiot. 

DIMENSION  I DAT A (6000) 


10 

20 


ID* 


M 

L  ■ 
IF 
1 
I 


.OT.  10)  L  •  -2 
*  L 


IF  (ABSIIDATA(I)  -  IX)  .LT.  M>  00  TO  100 


IF  (I  .OT. 
N  ■  2*N 
00  TO  10 
IF  <1  .LT. 
M  •  2*N 
00  TO  10 
1XYGET  •  I 
RETURN 
END 


I)  00  TO  30 


6000)  00  TO  20 


SUBROUTINE  ANTRIM! IDATA.L.N, AREA, PERU) 

C  *  Rot  till  to  oolooloto  oroo  ootf  poriaotor  of  a  portioa  of  tho  hood  with 
C  -  boaodrios  IDATA(L)  to  IDATA(N) . 

DIMENSION  1DATA(6000) 

AREA  •  0.0 

PERIM  -  0.0 

X2  •  FLOAT  < I DATA (L) > 

¥2  •  FLOAT! I DATA (L- 1 ) ) 

DO  20  l*L«2,N,2 
XI  -  X2 
¥1  •  ¥2 

X2  •  FLOAT! I DATA! I)) 

¥2  -  FLOAT ( I  DATA  < I - 1 ) ) 

BELT AX  -  X2  -  XI 
DELTAY  .  Y2  -  Yl 

PERU  <  PERU  ♦  SORT (DELTAX**2  ♦  DELTA¥**2> 

AREA  •  AREA  •  DOT  AX  •  I  (Yl#Y2)/2.0> 

20  CONTINUE 

DELTAX  »  FLOAT!  1  DAT  AIL)  -  IDATAU)) 

AREA  •  AREA  «  DELTAX  •  ((FLOAT! I DATA (N- !)♦ IDATA(L-I) ) )/2.0) 

RETURN 

END 


SUBROUTINE  SORTdDATA.  IDU.XPTS) 

•  tin  to  oort  data  oo  that  it  atvaya  olarta  at  (ha  lowor  wrial. 


10 

ion 

IIR 


2  ID 


DIMENSION  IDATA<6000) , IDIN(OOOO) 

00  10  I«NFTS,2j-2 

IF(IDIN(I)  .EQ.  1)  00  TO  100 
IF((NFTS-I)  .(71.  10)  00  TO  200 
DO  110  l« I ,NPTS 
I DATA! I )  •  ID1N(I) 

RETURN 
L  •  0 

00  210  E*I-l ,NPTS  (Start  alth  T 

L  •  L  ♦  1 
IDATA(L)  •  IDIN(K) 

CONTINUE 
DO  220  X-1,1-2 
L  •  L  *  I 
IDATA(L)  ■  IDIN(X) 

CONTINUE 

RETURN 


22D 


SUBROUTINE  TANANOdDATA, TANGLE, DLENG, DIST, 

♦JSTART , JSTEP , JRANGE.NP) 

C  *  loitix  to  atloiltlt  tha  u|U  fanotioa  arosad  •  iadtioa  of  tho 
C  “  poriaotor.  Tho  taolioa  it  -  to  ♦  DIST  froa  JSTART.  Tho  laafth  of 
C  -  tho  poriaotor  for  oooh  taagant  it  J RANGE,  tad  tanfanta  tro  ooloolotod 
C  -  for  otoh  taction  it  dttt  itlorotlt  JSTEP .  NP  it  tho  atabor  of  t«|lot 
C  "  ooloolotod  Analat  tro  ttorod  ia  TANGLE  tad  tho  oorroapoadiat 
C  -  loot th  it  DLENG. 

DIMENSION  IDATA(600O)  , TANGLE!  100)  ,DLENG(  100) 

DATA  PI/3. 1416/, RADIAN/S7. 296/ 

L  •  IDELTSdDATA, JSTART, DIST, -2) 

N  •  IDELTStIDATA, JSTART, DIST, 2) 

NP  •  0 
TLENG  •  0. 

DO  100  l.L.M, JSTEP 
NP  •  NP  ♦  I 

dx  •  av£( i data, i  orange)  -  aveiidata. i> 

DV  •  AVE< IDATA, loJRANQE- I >  -  AVEIIDATA, l-l) 

C  DX  -  FLOAT  ( IDATAd  ♦  JRANGE)  -  I  DATA!  I ) ) 

C  DY  •  FLOAT  1 1 DATA! 1 ♦JRANGE- l >  -  I DATA < I - 1 ) ) 

TLENG  -  (SORT (DX**2  ♦  DY**2>)/2.0  ♦  TLENG 

DLENG (NP)  •  TLENG 

IF  (OX  . NE.  0.)  GO  TO  20 

ANGLE  •  PI/2.  •  RADIAN 

IF  (DY  .LT.  0.0)  ANGLE  •  ANGLE  •  *1.0 

GO  TO  30 

20  ANGLE  •  AT AN2 ( DY , DX )  •  RADIAN 

30  IF  (NP  -BO.  I)  ANGLE!  •  ANGLE 

IF  (ABS(ANGLE-ANGLEl)  .LE.  100.)  GO  TO  40 
IF  (ANGLE!  .Of.  0.0)  ANGLE  •  ANGLE  *  300. 

IF  ( ANGLE 1  .LT.  0.0)  ANGLE  •  ANGLE  -  360. 
ao  tangle (Np>  -  angle 

ANGLE I  •  ANGLE 
100  CONTINUE 

RETURN 
END 


SUBROUTINE  LINFITIDLENG, TANGLE, NP,B,R) 

-  Roitiao  to  girt  a  loot)  sqoaroi  fit  to  dttt  with  a  litt  Y*A*BX. 

*  X  array  of  dttt  for  iadapaadaat  atritblt 

-  Y  array  of  data  for  doptadaat  tariabla 

-  NP  naaber  of  ptira  of  data  poiata 

-A  Y  intarcapt 

-  SIGNAA  standard  dariatioa  of  A 

“  B  tlopa 

-  SI GRAB  standard  daviatioa  of  B 

-  R  liaaar  corralation  ooafficiaat 

DIMENSION  TANGLE! 100) .DLENG! 100) 

SUM  -  FLOAT (NP) 

SUNK  •  0. 

SUMY  •  0. 

SUMX2  •  0. 

SUMY2  •  0. 

SUMXY  •  0. 

DO  100  1-1 ,MP 
XI  -  DLENG! 1) 

Yl  ■  TANGLEd) 

simx  -  sunx  -  xi 


SUMY  •  SUMY 
SUNX2  •  SUHX2 
SUMY2  •  SUMY2 
SUMXY  •  SUMXY 
CONTINUE 


Yl 


XI  «2 
Yl**2 
XI*Y1 


-  Calaalata  aoaffieiaait  tad  ataadard  dasialioaa 

DELTA  •  SUM  •  SUMX2-SUMX»»2 
A  •  (SUMX2*SUMY-SUMX*SUMXY) /DELTA 
t  •  (SUMXY*SUM -SUMX»SUMY ) /BELT A 
C  •  NP  -  2 

VARNCE  ■  iSU«T2-A,*2*SUM-B»*2«SUMX2-24 
I ( A*SUNY •B*SUNXY-A*B*SUNX> )/C 

SIGNAA  •  SQRT(VARNCE*SUNX2/DELTA) 

SIGMAB  •  SORT(VARNCE*SUM  /DELTA) 

R  •  (SUM*SUMXY-S  . 

RETURN 
END 


-SUNX*  SI' MY )  /SORT  ( DELTA*  (SUM*SUMT2-SUWT**2)  ) 
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DOCUMENTATION  OF  LINEBYLINE 


Gregory  K.  Myers 
SRI  International 
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I  INTRODUCTION 


This  algorithm  is  taken  from  "Image  Processing  Algorithms  for  In¬ 
dustrial  Vision,"  a  report  by  Gerry  Agin  of  SRI  International.  The 
algorithm  is  referred  to  as  "connectivity  analysis."  It  segments  a 
binary  image  into  "blobs"  (connected  regions)  of  the  same  "color" 

(gray  level).  Only  one  pass  is  made  through  the  image,  and  only  one 
line  is  accessed  at  a  time  (hence  the  name  "LINEBYLINE") .  Features  of 
each  blob  are  computed,  such  as  the  area,  center  of  gravity,  bounding 
rectangle,  and  perimeter  points.  This  document  explains  the  concepts 
used  in  LINEBYLINE  and  provides  an  example  of  processing. 


II  RUN-LENGTH  CODING 


The  first  step  in  processing  is  the  run-length  coding  of  the  binary 
image.  Each  line  of  the  binary  image  is  converted  into  a  series  of  run- 
length  segments.  A  "1"  denotes  the  color  of  "object"  and  "0"  denotes 
the  color  of  "background." 

To  illustrate  how  run-length  coding  is  done,  let  us  consider  the 


following  8x8  binary  image: 


Column  number: 


Row  number :  1 

2 

3 

4 

5 

6 

7 

8 


(Note:  "b"  denotes  an  implicit 


12345678 

bbbbbbbbbb 
b  0  0  0  1  1  1  0  0  b 

b  1  1  0  0  1  1  1  0  b 

b  0  1  0  1  1  0  1  1  b 

b  0  1  1  1  1  0  0  1  b 

b  0  1  1  1  0  0  0  1  b 

bOlOllOllb 
b  0  1  0  0  1  1  1  0  b 
bOOOOOlOOb 
bbbbbbbbbb 

0"  beyond  the  margins  of  the  image.) 
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FHECXOUO  not  1UNK-NOT  FILMED 


The  algorithm  for  run-length  coding  examines  each  pixel  of  the  row 
in  turn  and  records  the  column  number  of  any  element  that  differs  from 
its  predecessor.  A  0  is  assumed  to  precede  the  first  element  and  to 
follow  the  last  one.  For  this  reason,  if  the  last  column  of  some  par¬ 
ticular  row  contains  a  1,  a  transition  will  be  recorded  in  column  9. 

In  general,  there  will  be  twice  as  many  column  numbers  recorded  as 
there  are  segments  of  contiguous  Is  in  the  row. 

The  result  of  run-length  coding  our  example  is  as  follows: 


Row 

1: 

Col. 

4,7 

Row 

2: 

Col. 

1,3, 5, 8 

Row 

3: 

Col. 

2, 3, 4, 6, 7, 9 

Row 

4: 

Col. 

2, 6, 8, 9 

Row 

5: 

Col. 

2, 5, 8, 9 

Row 

6: 

Col. 

2, 3,4, 6, 7, 9 

Row 

7: 

Col. 

2, 3, 5, 8 

Row 

8: 

Col. 

6,7 

Later  in  this  paper,  we  will  refer  to  the  run-length  segments  that 
make  up  a  line  or  row.  These  are  contiguous  sequences  of  pixels  of 
either  color  (0  or  1).  For  example,  row  1  has  exactly  three  run-length 
segments.  The  first  segment  is  all  Os  and  extends  from  minus  infinity 
up  to  (but  not  including)  column  4.  The  second  segment  its  Is  runs  from 
column  4  to  (but  not  including)  column  7.  The  third  and  last  segment 
consists  of  Os  from  column  7  to  plus  infinity.  Because  every  row  must 
start  and  end  with  0,  there  is  always  an  odd  number  of  run-length  seg¬ 
ments  in  any  line. 

The  starting  and  ending  column  numbers  of  each  run-length  segment 
are  stored  in  a  linked  list  for  later  processing  by  the  connectivity- 
analysis  procedures.  When  the  connectivity  analysis  for  all  the  run- 
length  segments  in  a  line  has  been  completed,  the  run-length  segments 
are  deleted.  The  creation  of  run-length  segments  from  a  binary  image 
line  could  be  easily  performed  in  hardware. 


Ill  INTRODUCTION  TO  CONNECTIVITY  ANALYSIS 


The  purpose  of  connectivity  analysis  is  to  separate  a  binary  image 
into  connected  components.  We  call  these  connected  components  blobs . 

A  hole  is  a  special  case  of  a  blob  entirely  surrounded  by  another  blob 
of  the  contrasting  color.  When  connectivity  analysis  is  applied  to  the 
8x8  example  presented  in  the  previous  section,  the  result  contains 
three  components:  the  background  (component  0),  an  "object"  (component  1), 
and  a  hole  (component  2) .  The  three  components  are  shown  below. 

000  +  +  +00  ...111..  ...+  +  +.. 

+  +00  +  +  +0  11.  .111.  +  +..+  +  +. 

0  +  0  +  +  .  +  +  .1.11.11  .  +  .+  +2+  + 

0  +  +  +  +..  +  .1111.  .1  .+  +  +  +22  + 

0  +  +  +...+  .111.  ..1  .+  +  +222  + 

0  +  0  +  +.++  .1.11.11  .  +  .+  +2+  + 

0  +  00  +  +  +  0  .1..111.  .  +  ..+  +  +. 

00000  +  00  . 1..  . +  .. 

All  pixels  in  the  background  are  assigned  the  same  region  number,  even 
though  in  the  8x8  figure  the  background  appears  to  consist  of  four 
separate  regions.  This  is  because  we  assume  the  image  is  embedded  in 
an  infinite  field  of  0s.  Any  region  of  0s  that  touches  the  margin  of 
the  picture  will  therefore  be  classified  as  background. 

The  algorithm  extracts  and  identifies  connected  components,  com¬ 
putes  useful  feature  information,  and  requires  buffering  of  no  more 
than  one  line  of  data.  This  section  describes  the  general  procedure 
for  connectivity  analysis.  Section  IV  gives  some  details  that  ..re 
skipped  in  this  section  for  the  purpose  of  explaining  the  process  more 
easily. 

Connectivity  analysis,  as  implemented  here,  makes  use  of  a  data 
structure,  which  we  call  the  active  line,  to  keep  track  of  the  process¬ 
ing.  The  active  line  consist  of 


A  linked  list  of  segment  descriptors.  Each  segment  descriptor 
contain  the  starting  and  ending  column  numbers  of  a  run-length 
segment,  the  component  number  to  which  the  segment  belongs,  and 


a  pointer  to  the  next  segment  descriptor.  A  dummy  segment 
descriptor  gives  the  ending  column  number  of  the  final  run- 
length  segment. 

•  A  pointer  to  the  current  segment  being  processed. 

Before  processing  any  data,  the  active  line  is  initialized  to  con¬ 
tain  a  single  segment  descriptor,  representing  the  background.  Before 
processing  any  row  of  run-length  data,  the  pointer  to  the  current  segment 
pointer  must  be  initialized  to  point  to  the  first  segment.  As  we  start 
to  process  our  example  image,  the  state  of  affairs  is  thus: 

Image  so  far:  bbbbbbbbbb 

Active  line: 

Start  column  -infinity  +infinity 

Component  number  0  dummy 

Current  segment:  I 

In  the  first  row  of  example  data,  transitions  were  found  in  columns 
4  and  7.  This  implies  that  the  row  consists  of  three  run-length  segments 
Os  from  -infinity  to  (but  not  including)  column  4,  Is  from  column  4  to 
(but  not  including)  column  7,  and  Os  from  column  7  to  +infinity.  After 
processing  the  first  row  we  will  want  the  data  structure  to  represent 
those  three  segments. 

One  procedure  will  examine  each  run-length  segment  in  turn.  Each 
segment  will  take  its  turn  as  the  new  segment,  represented  by  its  start¬ 
ing  and  ending  column  numbers.  The  new  segment  must  be  matched  against 
the  partially  completed  analysis  embodied  in  the  current  line.  In  par¬ 
ticular,  we  must  determine  whether  the  new  segment  overlaps  the  current 
segment ,  that  is,  the  segment  of  the  active  line  pointed  to  by  the 
current  segment  pointer.  Three  cases  are  possible: 


The  two  segments  do  not  overlap  because  the  new  segment  Is  to  the 
right  of  the  current  one  (that  is,  the  starting  column  number  of 
the  new  segment  is  larger  than  the  ending  column  number  of  the 
current  segment) . 

Current  segment:  XXXXXXXX 

New  segment :  XXXXXXXX 

Case  2: 

The  two  segments  do  not  overlap  because  the  new  segment  is  to  the 
left  of  the  current  one  (that  is,  the  ending  column  number  of  the 
new  segment  is  smaller  than  the  starting  column  number  of  the 
current  segment). 

Current  segment :  XXXXXXXX 

New  segment :  XXXXXXXX 

Case  3: 

Neither  Case  1  nor  Case  2  obtain,  and  the  segments  overlap. 

Current  segment :  XXXXXXXXXXXX 

New  segment :  XXXXXXXXXXX 

Different  actions  must  be  taken  in  each  of  the  three  cases.  We  will 
explain  each  case  as  it  occurs  in  the  course  of  analyzing  our  example 
image. 

The  first  segment  of  the  new  row  goes  from  -infinity  to  4.  The 
current  segment  pointer  points  to  the  first  segment  in  the  active  line, 
which  goes  from  -infinity  to  +infinity.  When  these  two  segments  are 
compared,  we  find  that  Case  3  applies:  the  two  segments  overlap. 

Action  on  Case  3: 

Copy  the  start  column  number  from  the  new  segment  to  the  current 
one.  Then  advance  the  current  segment  pointer  to  the  next  segment 
in  the  active  line. 


Copying  the  starting  number  of  the  new  segment  (-infinity)  to  the 
current  segment  results  in  no  change  in  the  data  structure,  because 
-infinity  was  there  to  begin  with.  The  current  segment  pointer  is  ad¬ 
vanced  to  the  next  segment  descriptor  in  the  active-line  list.  After 
processing  the  first  segment  in  the  first  row,  the  active-line  data 
structure  looks  like  this: 


Image  so  far: 
New  segment: 


bbbbbbbbbb 
b  0  0  0 


Active  line: 

Start  -inf  +inf 

Component  0  dummy 


Current  segment: 


I 


We  now  analyze  the  next  segment.  This  new  segment  runs  from 
columns  4  up  to  (but  not  including)  column  7.  The  current  segment 
pointer  now  points  to  the  dummy  segment  descriptor,  which  starts  at 
+infinity.  Therefore  Case  2  applies:  The  new  segment  does  not  match 
any  existing  segment  in  the  active  line.  Room  must  be  made  in  the  data 
structure  for  it. 

Action  on  Case  2: 

Insert  two  new  segments  descriptors  (call  them  A  and  B)  in  the 
active  line  before  the  current  segment.  Let  the  current  segment 
pointer  point  to  A.  Copy  the  ending  column  number  of  the  new 
segment  to  the  start  column  number  of  B.  Choose  an  unused  number 
for  a  new  component  and  let  it  be  the  component  number  for  A.  Copy 
the  component  number  from  the' segment  preceding  segment  A,  to  B. 

Now  proceed  as  in  Case  3. 

The  situation  after  processing  the  second  segment  in  the  first  row 
is  thus: 


Image  so  far: 

b  b 

b  b  b  b  b 

b  b  b 

b  0 

0  0 

New  segment: 

111 

Active  line: 

(A) 

(B) 

Start 

-inf 

4 

7 

+inf 

Component 

0 

1 

0 

dummy 

Current  segment : 

1 

Two  new  segment  descriptors  were  inserted  in  the  active  line.  Descriptor 
A  represents  part  of  a  newly  discovered  blob  in  the  image;  it  is  assigned 
component  number  1.  Descriptor  B  extends  the  background  down  the  right 
side  of  the  new  blob. 


Descriptor  B  takes  its  starting  column  number  from  the  ending 
column  number  of  the  new  segment  (7)  and  its  component  number  from  the 
segment  before  A  (0).  Having  created  a  new  pair  of  segment  descriptors, 
we  can  proceed  to  match  the  new  segment  to  segment  A  by  performing  the 
Case  3  action.  The  starting  column  number  of  the  new  segment  (4)  is 
copied  to  descriptor  A,  and  the  current  segment  pointer  is  advanced  to 
point  to  segment  B. 

The  third  (and  last)  segment  in  the  first  row  runs  from  column  7 
to  +infinity.  This  will  be  matched  against  the  current  segment,  which 
also  runs  from  7  to  +infinity.  This  is  another  Case  3  overlap,  and  the 
data  structure  after  processing  that  new  segment  will  be: 

Image  so  far :  bbbbbbbbbb 

b  0  0  0  1  1  1 


New  segment : 

0  0b 

Active  line: 

Start 

-inf 

4 

7 

+inf 

Component 

Current  segment: 

0 

1 

0 

dummy 

1 

This  concludes  processing  of  the  first  row. 

The  second  row  has  transitions  at  1,  3,  5,  and  8,  or  five  segments 
to  be  processed.  The  current  segment  pointer  is  reset  to  point  to  the 
first  active  segment. 

The  first  segment  of  this  row  runs  from  -infinity  to  1.  When  it 
is  matched  to  the  first  segment  in  the  active  line,  a  Case  3  overlap  is 
discovered.  The  appropriate  action  (extending  the  segment  representing 
the  background)  will  be  performed. 

The  next  segment  runs  from  column  1  to  (but  not  including)  column  3. 
When  it  is  matched  against  the  current  segment  starting  at  column  4, 

Case  2  applies.  The  newly  discovered  blob  is  assigned  component  number  2, 
producing  the  following  situation  after  processing  the  second  segment: 


Image  so  far: 


bbbbbbbbbb 
bOOOlllOOb 
b 

New  segment:  2  2 

Active  line:  (A)  (B) 

Start  -inf  1  3  4  7  +inf 

Component  02010  dummy 

Current  segment:  ! 


The  last  three  segments  on  this  second  row  all  produce  Case  3  over 
laps  when  they  are  matched  against  the  active  line.  At  the  conclusion 
of  processing  the  second  row,  the  following  data  structure  obtains: 


Image  so  far:  bbbbbbbbbb 

bOOOlllOOb 
b  2  2  0  0  1  1  1 

New  segment:  0  b 

Active  line: 

Start  -inf  1  3  5  8  +inf 

Component  0  2  0  1  0  dummy 

Current  segment :  ® 


Row  3  has  transitions  at  columns  2,  3,  4,  6,  7,  and  9.  Matching 
the  seven  segments  of  this  row  will  result  in  four  Case  3s,  a  Case  2, 
and  two  more  Case  3s.  After  processing  the  third  row  the  following 
will  be  the  state  of  the  data  structure: 


Current  segment: 


In  the  fourth  row  there  are  transitions  at  columns  2,  6,  8,  and  9, 
for  a  total  of  five  segments.  The  first  two  of  these  result  in  Case  3 
overlaps.  But  when  we  match  the  third  segment  of  the  fifth  row  (which 
runs  from  column  6  to  8)  with  the  next  segment  of  the  active  line 
(running  from  column  3  to  column  4)  we  discover  that  Case  1  holds.  Here 
is  the  situation  as  we  discover  the  fact: 


Image  so  far:  bbbbbbbbbb 

b  0  0  0  1  1  1  0  0  b 
b22001110b 
b  0  2<0>1  1(3)1  1  b 
b  0(2  2  2  2) 


New  segment: 

[0  0] 

Active  line: 

(A) 

(B) 

Start  -inf 

2 

3 

4  6 

7 

9 

+inf 

Component  0 

2 

0 

1  3 

1 

0 

dummy 

Current  segment : 


The  current  segment  in  the  active  line,  marked  by  "<  >"  above,  is 
not  matched  by  any  segment  in  the  new  row.  The  new  segment,  marked  by 
"[  ]"  must  eventually  match  the  segment  denoted  by  "(  )".  We  must 
delete  the  unmatched  segment  "<  >"  from  the  active  line.  Furthermore, 
the  previous  segment,  marked  by  "(  )",  forms  a  "bridge"  between  com¬ 
ponents  1  and  2.  We  must  merge  the  two  components. 

Action  on  Case  1: 

Consider  the  segment  in  the  active  line  before  the  current  one 
(call  it  "A")  and  the  segment  after  the  current  one  (call  it 
"B").  If  the  component  numbers  of  A  and  B  are  different,  then 
merge  the  two  components  by  changing  all  instances  of  B's  component 
number  in  the  active  line  to  A's  component  number.  Delete  the 
current  segment  and  segment  B  from  the  active  line.  Let  the  cur¬ 
rent  segment  pointer  point  to  the  segment  after  B.  Now  go  back 
to  the  beginning  of  the  entire  matching  procedure,  matching  the 
new  segment  against  the  updated  active  line. 

By  this  procedure,  we  merge  components  1  and  2.  The  merging  changes 
the  component  numbers  of  the  fourth  and  sixth  segment  descriptors  in  the 
active  line  from  1  to  2.  The  current  segment  and  segment  B  get  deleted, 
and  the  current  segment  pointer  is  advanced  to  point  to  the  segment 


starting  at  column  6.  After  ve  have  taken  these  steps,  the  data  struc¬ 
ture  looks  as  follows: 


b  b  b  b 
2  0  0  b 

2  2  0  b 

3  2  2  b 

3  3 


Image  so  far: 

New  segment: 

Active  line: 

Start 

Component 

Current  segment : 


b  b  b  b  b  b 
b  0  0  0  2  2 
b  2  2  0  0  2 
b  0  2  0  2  2 
b  0  2  2  2  2 


(A) 

-inf  2  6 

0  2  3 

I 


7  9  +inf 

2  0  dummy 


Now  the  new  segment  may  be  matched  anew  against  the  active  line  and 


current  segment.  This  time  around,  the  match  is  Case  3,  as  are  the  re¬ 
mainder  of  matches  in  this  fourth  row.  After  finishing  this  row,  the 
state  of  the  data  structure  is  as  shown  below: 


Image  so  far: 


Active  line: 

Start 

Component ; 
Current  segment: 


b  b 
b  0 
b  2 
b'  0 
b  0 


-inf 

0 


b  b  b  b 
0  0  2  2 
2  0  0  2 
2  0  2  2 
2  2  2  2 


2 

2 


b  b  b  b 
2  0  0  b 

2  2  0  b 

3  2  2  b 
3  3  2  b 


6 

3 


8 

2 


9 

0 


+inf 

dummy 

I 


In  the  fifth  row,  there  are  five  segments  with  transitions  at 


columns  2,  5,  8,  and  9.  Only  Case  3  overlaps  occur. 


The  "image  so  far"  is  not  part  of  the  data  structure;  it  only  is  pre¬ 
sented  as  an  aid  in  visualization  of  the  connectivity-analysis  process. 
We  have  changed  all  Is  in  the  image  to  2s  to  indicate,  conceptually,  the 
merging  of  the  two  components. 


The  sixth  row  has  transitions  at  columns  2,  3,  4,  6,  7,  and  9.  The 
third  segment  (from  columns  3  to  4)  will  generate  a  Case  2  match,  the 
others  will  all  be  Case  3.  The  new  component  created  by  the  Case  2 
match  is  assigned  number  4.  Here  is  what  the  data  structure  looks  like 
after  processing  that  row: 


Image  so  far: 


bbbbbbbbbb 
b  0  0  0  2  2  2  0  0  b 
b  2  2  0  0  2  2  2  0  b 
b02022322b 
b  0  2  2  2  2  3  3  2  b 
b  0  2  2  2  3  3  3  2  b 
b  0  2  4  2  2  3  2  2  b 


Active  line: 
Start  -inf 
Component  0 

Current  segment : 


+inf 

dummy 


In  processing  row  7,  with  transitions  at  2,  3,  5,  and  8  a  Case  1 
situation  occurs.  The  analysis  of  component  3  is  complete:  Its  segment 
descriptor  disappears  from  the  active-line  data  structure.  The  two 
segments  on  either  side  of  component  3  already  have  the  same  component 
numbers,  so  no  merging  is  necessary: 


Image  so  far: 


bbbbbbbbbb 
b  0  0  0  2  2  2  0  0  b 
b  2  2  0  0  2  2  2  Ob 
b  0  2  0  2  2  3  2  2  b 
b  0  2  2  2  2  3  3  2  b 
b  0  2  2  2  3  3  3  2  b 
b  0  2  4  2  2  3  2  2  b 
b  0  2  4  4  2  2  2  0  b 


Active  line: 

Start 

Component 


+inf 

dummy 


Current  segment: 


In  row  8,  there  are  transitions  at  columns  6  and  7.  Another  Case  1 
situation  occurs  here,  merging  component  4  with  component  0  (the  back¬ 
ground).  Here  is  the  situation  after  we  have  processed  the  eighth  row: 

Image  so  far:  bbbbbbbbbb 

b  0  0  0  2  2  2  0  0  b 

b  2  2  Q  O  2  2  2  0  b 

b  0  2  0  2  2  3  2  2  b 

b  0  2  2  2  2  3  3  2  b 

b  0  2  2  2  3  3  3  2  b 

b  0  2  0  2  2  3  2  2  b 

b  0  2  0  0  2  2  2  0  b 

b  0  0  0  0  0  2  0  0  b 


Active  line: 

Start 

-inf 

6 

7 

+inf 

Component 

0 

2 

0 

dummy 

Current  segment :  1 

After  the  last  row  of  image  data,  another  row  of  all  Os  will  com¬ 
plete  the  analysis,  removing  component  2  from  the  active  line. 

We  have  engaged  in  a  little  oversimplification  in  order  to  present 
the  basic  ideas  behind  connectivity  analysis  as  plainly  as  possible. 

In  the  next  section  we  will  fix  the  algorithm  to  operate  correctly  in 
all  cases.  But  before  reading  on,  turn  to  the  beginning  of  this  section, 
where  we  presented  the  three  components  we  wish  the  algorithm  to  extract. 
You  may  verify  that  the  results  we  have  obtained  are  what  was  intended 
except  that  different  numbers  were  assigned. 

IV  CONNECTIVITY  ANALYSIS  IN  DETAIL 

The  simplified  algorithm  presented  in  the  previous  section  is  de¬ 
ficient.  It  has  errors,  and  it  doesn't  compute  any  features  of  the  blobs 
or  any  perimeter  lists.  In  addition,  processing  of  some  borderline 
cases  of  overlapping  segments  was  not  adequately  described.  At  the  end 
of  this  section  the  complete  algorithm  with  its  amendments  and  corrections 
will  be  presented. 
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A. 


Handling  Some  Special  Cases 


Let  us  finish  the  processing  of  our  test  image.  To  complete  the 
analysis,  a  last  row  of  Os  running  from  -infinity  to  +infinity  is  added 
to  the  active-line  data  structure,  leaving  us  with  the  following  situa¬ 
tion: 

Image  so  far: 


New  segment: 


Active  line: 

Start 

-inf 

6 

7 

+inf 

Component 

0 

2 

0 

dummy 

Current  segment:  I 


We  have  just  processed  a  line  with  nothing  in  it,  but  the  active 
line  incorrectly  points  to  a  segment  there.  Tt  takes  a  Case  1  match 
between  segments  to  delete  component  2  from  the  data  structure,  but 
there  is  nothing  left  to  process.  What  we  must  do  is  to  check  at  the 
end  of  each  line  whether  the  current  segment  pointer  points  to  the  dummy 
segment  at  the  end.  If  it  does  not,  then  the  Case  1  processing  must  be 
performed  until  the  current  segment  pointer  does  point  to  the  last  seg¬ 
ment. 

A  second  error  concerns  the  merging  of  two  regions  when  a  "bridge" 
connects  them.  The  rule  for  Case  1  states  that  if  the  two  component 
numbers  are  different,  then  the  number  on  the  right  replaces  the  number 
on  the  left  everywhere  else  in  the  active  line.  However,  if  the  com¬ 
ponent  number  on  the  right  happens  to  be  the  background,  the  replacement 
would  renumber  the  background.  We  wish  the  background  to  be  always 
component  number  0  and  must  modify  the  renumbering  rule  accordingly. 


bbbbbbbbbb 
b  0  0  0  2  2  2  0  0  b 
b  2  2  0  0  2  2  2  Ob 
b  0  2  0  2  2  3  2  2  b 
b  0  2  2  2  2  3  3  2  b 
b  0  2  2  2  3  3  3  2  b 
b  0  2  0  2  2  3  2  2  b 
b  0  2  0  0  2  2  2  Ob 
b000002Q0b 
bbbbbbbbbb 


B.  Blob  Descriptors 

The  main  reason  for  doing  connectivity  analysis  is  to  obtain  infor¬ 
mation  about  the  components  we  find  thereby.  The  information  we  have 
derived  related  to  color,  location,  area,  perimeter,  and  points  to  other 
components  that  surround,  neighbor,  or  form  holes  in  the  component  we 
are  examining.  The  algorithms  for  deriving  such  information  will  be 
described  later  in  this  report.  Such  information  is  stored  in  a  blob 
descriptor  for  each  component. 

We  must  modify  our  algorithm  description  to  include  the  use  of  blob 
descriptors.  When  a  new  component  is  started  as  a  result  of  Case  2 
processing,  a  new  blob  descriptor  must  be  created.  When  segments  are 
removed  from  the  active  line  in  Case  1  processing,  we  must  keep  track 
of  the  blob  descriptors  that  are  removed.  Several  segment  descriptors 
in  the  active-line  data  structure  can  have  the  same  component  number; 
hence  they  can  point  to  the  same  blob  descriptor.  Therefore,  when  a 
segment  descriptor  is  removed  as  the  result  of  Case  1  processing,  the 
algorithm  must  check  to  see  if  other  segment  descriptors  in  the  active 
line  refer  to  the  same  component.  If  this  is  the  last  reference  to 
that  component,  its  processing  is  finished. 

At  this  point,  the  blob  descriptor  describes  a  connected  component 
whose  connectivity  has  been  completely  determined.  Furthermore,  the 
extraction  of  blob  features  and  perimeter  lists  is  complete.  We  have 
more  than  one  option  as  to  what  to  do  with  the  blob  descriptor.  We  can 
pass  the  address  of  the  completed  blob  descriptor  to  an  application- 
dependent  subroutine  that  will  do  further  processing,  take  some  action 
depending  on  what  was  found,  or  delete  the  blob  descriptor  to  make  its 
memory  space  available  for  other  blobs  or  other  purposes.  Alternatively, 
the  address  of  the  blob  descriptor  can  simply  be  added  to  a  list  of 
isolated  blobs  in  the  image  and  further  processing  will  be  deferred 
until  the  entire  picture  is  analyzed.  In  the  current  implementation  the 
features  and  perimeter  lists  of  the  blob  are  written  to  output  files, 
and  the  blob  descriptor  is  recycled  to  save  memory  space. 


B-16 


C.  Treatment  of  Diagonally  Adjacent  Pixels 

Two  pixels  belong  to  the  same  component  if  and  only  if  there  exists 
a  path,  along  adjacent  elements  of  the  same  color,  from  one  pixel  to  the 
other.  When  dealing  with  points  on  a  rectangular  grid,  it  is  reasonable 
to  ask  whether  diagonal  points  are  adjacent.  If  diagonal  points  are  to 
be  considered  adjacent,  then  the  components  of  an  image  are  said  to  be 
8-connected ,  since  each  pixel  has  eight  connected  neighbors.  If  diagonal 
points  are  not  considered  adjacent,  then  the  components  are  called  4- 
connected .  Our  connectivity  algorithm,  as  presented  so  far,  finds  the 
4-connected  components  of  an  image.  If  we  were  to  analyze  the  image  of 
a  checkerboard,  then  each  square  of  the  image,  black  or  white,  would  be 
a  separate  component. 

It  would  be  desirable  to  modify  the  connectivity  algorithm  so  it 
has  the  property  that  every  component  (except  the  background)  is  en¬ 
tirely  enclosed  by  one  and  only  one  component,  which  is  of  the  opposite 
color.  Then  every  boundary  between  two  regions  is  a  closed  curve  that 
separates  an  enclosing  region  from  an  enclosed  region.  The  components 
of  a  scene  thus  analyzed  may  be  arranged  hierarchically,  each  component 
having  a  single  superior  (encloser)  and  zero  or  more  inferiors  (holes 
or  enclosed  regions) . 

A  way  of  achieving  this  is  to  specify  that  all  white  cells  (Is) 
are  8-connected  and  that  all  black  cells  (Os)  are  4-connected.  With 
such  a  convention,  the  image  of  a  checkerboard  would  be  analyzed  as  a 
number  of  individual  white  squares,  all  embedded  in  a  single  black 
region  (the  background) .  To  implement  this  requires  modifying  the 
tests  for  overlap  that  determine  Case  1,  2,  or  3.  The  modification  has 
to  be  such  that  one  set  of  tests  applies  when  a  segment  of  Is  is  being 
considered  and  a  different  set  of  tests  applies  for  a  segment  of  Os. 

D.  The  Connectivity  Algorithm 

We  shall  now  consolidate  what  has  gone  before.  The  algorithm 
will  be  presented  in  an  informal  fashion.  Operations  such  as  creating 
and  modifying  lists,  or  obtaining  storage  for  a  descriptor  are  imple¬ 
mentation  details  that  will  not  be  mentioned  here. 
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First,  some  definitions. 


•  Segment  descriptor  is  an  array  or  a  block  of  data  that  contains 
four  items:  starting  and  ending  column  numbers,  a  pointer  to  a 
blob  descriptor,  and  a  pointer  to  the  next  segment  descriptor 
in  the  active  line. 

•  The  active  line  is  a  linked  list  of  segment  descriptors.  The 
current  segment  pointer  points  to  a  segment  descriptor  on  this 
list.  The  active  segment  is  the  segment  descriptor  pointed  to 
by  the  current  segment  pointer. 

•  A  blob  descriptor  is  a  record  of  data  that  contains  at  least  one 
item:  a  color,  which  may  be  either  0  or  1.  Additional  itmes 

in  a  blob  can  be  used  for  feature  analysis,  as  will  be  described 
later.  A  component  number  is  a  number  that  identifies  a  blob. 

To  process  an  image: 

•  Obtain  a  blob  descriptor  to  represent  the  background.  Set  the 
color  word  of  the  background  blob  to  0. 

•  Initialize  the  active  line  to  contain  two  segment  descriptors. 

The  first  segment  descriptor  should  have  a  column  number  smaller 
than  zero  and  a  component  number  pointing  to  the  background 
blob.  The  column  number  of  the  second  segment  descriptor  should 
be  a  large  positive  number,  a  number  greater  than  the  number 

of  columns  in  the  image.  The  component  number  of  the  second 
descriptor  is  irrelevant. 

•  Process  each  row  of  the  image,  as  described  below. 

•  Finish  by  processing  an  extra  row  consisting  of  all  zeros. 

To  process  a  row: 

•  Initialize  the  current  segment  pointer  to  point  to  the  first 
segment  descriptor  in  the  active  line. 

•  Obtain  the  run-length  representation  of  the  row  (see  Section 
II).  The  run-length  data  should  start  with  a  negative  number 
and  end  with  a  large  positive  number. 

For  every  pair  of  adjacent  numbers  in  the  run-length  data,  in 
turn,  call  the  segment-processing  operation  (defined  below). 

While  the  current  segment  pointer  does  not  point  to  the  last 
segment  in  the  active  line,  perform  the  deletion  operation  (de¬ 
fined  below).  Repeat  this  step  until  the  current  segment  pointer 
does  point  to  the  last  segment. 


In  this  context,  the  word  while  refers  to  a  condition  to  check  and  an 
action  to  perform.  It  means  to  repeat  the  action  zero  or  more  times 
until  the  condition  is  not  true. 


To  process  a  new  segment,  given  a  starting  column  number  and  an 
ending  column  number : 

•  While  the  starting  column  number  is  greater  than  the  ending 
column  number  of  the  current  segment  in  the  active  line,  do  the 
deletion  operation  (defined  below).  (This  is  Case  1.) 

•  If  the  ending  column  number  is  less  than  the  starting  column 
number  of  the  current  segment,  perform  the  insertion  operation 
(defined  below),  passing  on  the  starting  and  ending  numbers  of 
the  segment.  (This  is  Case  2.) 

•  Do  feature  extraction  for  Case  3  processing. 

•  Copy  the  starting  and  ending  column  numbers  of  the  new  segment 

to  the  starting  and  ending  column  numbers  of  the  current  segment. 

•  Advance  the  current  segment  pointer  to  point  to  the  next  segment 
descriptor  in  the  active  line. 

To  perform  the  insertion  operation,  given  starting  and  ending 
column  numbers: 

•  Obtain  the  component  number  of  the  segment  descriptor  preceding 
the  current  segment.  Call  that  component  the  surrounding  com¬ 
ponent  . 

•  Obtain  a  new  blob  descriptor.  Call  it  the  new  component.  Set 
the  color  word  of  the  new  component  to  the  opposite  of  the  color 
of  the  surrounding  component. 

•  Obtain  two  new  segment  descriptors  and  insert  them  in  the  active 
line  immediately  before  the  current  segment,  calling  the  first 
segment  "A"  and  the  second  "B."  Segment  A  receives  the  new  com¬ 
ponent  number  and  the  starting  column  number.  Segment  B  re¬ 
ceives  the  surrounding  component  number  and  the  ending  column 
number . 

To  perform  the  deletion  operation: 


•  Call  the  component  number  of  the  current  segment  the  terminated 
component .  Call  the  component  number  of  the  segment  preceding 
the  current  one  the  left  component  and  the  component  number  of 
the  segment  following  the  current  one  the  right  component.  Call 
the  left  component  the  replacing  component  and  the  right  component 
the  replaced  component.  If  the  right  component  points  to  the 
background,  then  call  it  the  replacing  component  and  the  left 
component  the  replaced  component;  otherwise,  call  the  left 
component  replacing  and  the  right  component  replaced. 

•  Do  feature  extraction  for  Case  1  processing. 

•  If  the  replacing  component  and  the  replaced  component  are  dif¬ 
ferent,  then  find  all  instances  of  the  replaced  component  number 
in  the  active  line,  and  change  them  to  the  replacing  number. 


•  Delete  the  current  segment  and  the  segment  following  the  cur¬ 
rent  one  from  the  active  line.  Let  the  current  segment  pointer 
point  to  the  first  segment  after  the  deleted  one. 

•  Search  for  instances  of  the  terminated  component  number  in  the 
active  line.  If  there  are  no  remaining  instances,  call 
application-dependent  subroutines,  as  appropriate,  passing  the 
address  of  the  terminated  component's  blob  descriptor. 


V  EXTRACTION  OF  FEATURES 

The  preceding  sections  outlined  a  procedure  for  isolating  connected 
components  (blobs)  from  one  another  in  an  image.  The  results  of  the 
analysis  are  a  number  of  blob  descriptor  records  that  contain  information 
about  the  blob.  Each  blob  descriptor  consists  of  several  features  or 
characteristics  of  interest. 

A.  Color 

The  simplest  feature  is  color.  The  color  of  a  blob  will  be  0 
if  the  blob  is  black  and  1  if  it  is  white.  If  the  color  of  the  background 
is  appropriately  set  up  at  initialization  time,  then  the  following  will 
ensure  that  each  blob  descriptor  subsequently  created  has  the  correct 
color:  whenever  a  new  blob  descriptor  is  created  (Case  2)  obtain  the 
color  feature  of  the  surrounding  blob  and  store  the  opposite  color  in 
the  appropriate  place  in  the  new  blob  descriptor. 

B.  Parent,  Child,  and  Sibling 

An  important  class  of  features  describes  inclusion  relationships 
among  blobs.  Because  we  have  been  careful  in  the  way  connections  are 
made  between  diagonally  adjacent  pixels,  we  can  guarantee  that  every 
component  (except  the  background)  has  one  and  only  one  surrounding 
component  blob.  We  call  any  blob's  surrounder  its  parent ,  any  holes 
(or  blobs  enclosed  by  this  one)  children,  and  other  blobs  enclosed  by 
the  same  parent  siblings. 

At  blob  creation  time  (Case  2)  a  pointer  to  the  surrounding  blob 
descriptor  can  be  stored  as  the  new  blob's  parent  feature.  The  child 


feature  of  the  new  blob  is  set  to  0.  We  must  also  record  the  fact  that 
the  parent  blob  has  an  additional  child.  There  is  only  a  single  entry 
in  the  blob  descriptor  for  the  child  feature.  The  child  feature  of  any 
blob  points  to  the  most  recently  created  surrounded  blob;  the  sibling 
feature  of  that  most  recently  created  child  points  to  the  next  most 
recently  created  child,  and  so  forth.  The  procedure  to  follow  in  Case  2 
processing  is  as  follows:  Copy  the  child  feature  of  the  surrounding 
blob  to  the  sibling  feature  of  the  new  blob.  Set  the  child  feature  of 
the  surrounding  blob  to  point  to  the  new  blob. 

Case  1  processing  frequently  causes  two  blobs  to  be  merged  into 
one.  Topological  considerations  guarantee  that  the  two  blobs  to  be 
merged  will  have  the  same  color  feature  and  the  same  parent.  But  to 
preserve  the  integrity  of  the  parent-child-sibling  links,  several 
actions  must  be  performed.  In  each  of  the  children  of  the  replaced  blob, 
the  parent  feature  must  be  changed  to  point  to  the  replacing  blob.  All 
the  children  of  the  replaced  blob  must  be  concatenated  onto  the  sibling 
list  of  the  children  of  the  replacing  blob.  The  pointer  to  the  replaced 
blob  itself  must  be  removed  from  the  sibling  list  that  includes  both 
the  replaced  and  the  replacing  blobs.  Finally,  if  the  common  parent  of 
the  two  merged  blobs  points  to  the  replacing  blob  as  its  child  feature, 
that  child  feature  must  be  replaced  by  the  contents  of  the  sibling 
feature  of  the  replaced  blob. 

A  greal  deal  of  bookkeeping  overhead  is  eliminated  by  recording  only 
parent  relationships.  The  children  and  siblings  can  be  regenerated 
from  the  list  of  parent  relationships.  In  the  current  implementation  of 
LINEBYLINE,  the  parent  relationship  is  not  computed  because  it  is  not 
used  by  subsequent  programs. 


C.  The  Bounding  Rectangle 

The  bounding  rectangle  of  a  blob  is  specified  by  the  minimum  and 
maximum  values  of  x  and  y  over  the  blob.  It  is  easily  extracted  as 
follows:  when  a  blob  is  created  (Case  2),  copy  the  current  value  of  y 
(which  is  constant  across  an  entire  line  of  input)  to  the  new  blob's 
ymin  feature.  At  the  same  time  place  the  starting  column  number  (of  the 


contrasting  run-length  segment)  in  the  new  blob's  xrom  feature  and  the 
ending  column  number  in  the  blob’s  ymin  feature. 


For  each  run-length  segment  added  to  a  blob  (Case  3) ,  set  the  xmin 
feature  to  the  smaller  of  (1)  the  starting  column  number  and  (2)  the 
previous  value  of  the  xmin  feature.  Do  the  analogous  operation  for  the 
xmax  feature  and  the  ending  column  number.  The  ymax  feature  is  set  to 
the  current  line  number. 

When  two  blobs  are  merged  as  a  result  of  Case  1  processing,  the 
ymin,  xmin,  and  xmax  features  of  the  replacing  blob  must  take  the  ex¬ 
treme  values  of  the  two  blobs  to  be  merged. 

D.  Area  and  Center  of  Gravity 

The  area  of  a  blob  is  set  to  0  when  the  blob  descriptor  is  created. 

When  a  new  line  segment  of  the  blob  is  processed,  the  number  of  pixels 

in  the  segment  is  added  to  the  area  of  the  blob.  When  two  blob:,  are 

merged,  their  areas  are  added  together. 

The  center  of  gravity  (CG)  is  recomputed  after  each  new  line  segment 
is  added  to  the  blob.  The  updated  CG  is  a  weighted  average  of  the  CG 
of  the  new  line  and  the  CG  of  the  previously  processed  portion  of  the 
blob.  The  weights  are  proportional  to  the  areas  of  the  new  line  segment 
and  the  previously  processed  portion,  respectively.  When  two  blobs  are 
merged,  a  weighted  average  of  the  two  CG's  is  computed  in  a  similar 
manner  as  above. 

E.  Perimeter  Lists 

Associated  with  each  inked  blob  is  a  set  of  perimeter  points.  Be¬ 
cause  each  line  of  the  image  is  processed  sequentially,  the  perimeter 
points  are  obtained  in  a  somewhat  haphazard  manner.  During  processing, 
the  perimeter  of  a  simple  convex  blob  has  two  uncompleted  ends:  One  at 
the  starting  column  of  the  active  line  segment  (the  left  end),  and  one 
at  the  ending  column  of  the  active  line  segment  (the  right  end).  The 
perimeter  points  are  stored  as  a  singly  linked  list  called  a  perimeter 
section.  The  right  endpoint  is  at  the  beginning  of  the  list,  and  the 


r 


left  endpoint  is  at  the  end  of  the  list.  When  there  is  more  than  one 
active  line  segment  corresponding  to  the  same  blob,  one  perimeter  section 
is  associated  with  each  active  line  segment.  The  perimeter  sections 
are  arranged  as  a  doubly  linked  ring.  This  structure  is  illustrated 
below,  using  component  2  in  the  previous  example  after  six  lines  of 
processing. 
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The  structure  of  the  perimeter  section  records  and  perimeter  point 
records  in  the  illustration  above  are  shown  here: 


left 

right 

line 

previous  section 

column 

next  section 

r.nup.i.H.ui 

Perimeter  section  record  Perimeter  point  record 

When  each  new  line  of  pixels  is  processed,  the  lists  of  perimeter 
points  are  modified  according  to  which  case  applies.  In  all  cases 
perimeter  lists  are  compiled  for  inked  blobs  only. 

Case  3  requires  the  simplest  processing.  When  a  new  line  segment 
is  added  to  an  inked  blob,  the  left  end  of  the  associated  perimeter 
section  and  the  right  end  of  the  next  perimeter  section  are  extended. 
An  example  is  shown  below: 


Image  so  far: 
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New  segment:  [111] 


Perimeter  lists 
before  new  seg¬ 
ment  is  added: 


Perimeter  lists 
after  new  seg¬ 
ment  is  added: 
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For  Case  2  (insertion  of  a  new  segment)  there  are  two  subcases. 
Subcase  A  occurs  if  the  new  segment  has  the  color  ink;  otherwise  Sub¬ 
case  B  occurs.  In  Subcase  A  a  new  perimeter  section  for  the  new  blob  is 
created  and  will  contain  all  of  the  points  in  the  new  segment.  An 
example  is  shown  below. 


Image  so  far: 


New  segment: 
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Perimeter  list  after 
new  segment  is  added: 


In  Subcase  B  the  perimeter  section  of  the  surrounding  segment  is 
split  into  two  perimeter  sections.  An  example  is  shown  below. 


Image  so  far: 


New  segment: 
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Perimeter  lists 
before  new  seg¬ 
ment  is  added: 
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Perimeter  lists 
after  new  seg¬ 
ment  is  added: 
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For  Case  1  (deletion  of  a  segment)  there  are  also  two  subcases. 
Subcase  A  occurs  if  the  deleted  segment  has  the  color  ink;  otherwise 
Subcase  B  occurs.  In  Subcase  A  two  adjacent  perimeter  sections  of  the 
terminating  blob  are  merged.  If  the  terminating  blob  has  only  one 
perimeter  section,  its  left  and  right  ends  are  joined,  and  processing 
of  this  blob  is  finished.  An  example  is  shown  below. 


Image  so  far: 


New  segment: 


bbbbbbbbbb 
b  0  0  1  1  1  1  0  0  b 
b  0  1  1  1  1  1  0  0  b 
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Perimeter  lists 
before  new  seg¬ 
ment  is  added: 
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Perimeter  lists 
after  new  seg¬ 
ment  is  added: 
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In  Subcase  B  the  perimeter  sections  associated  with  the  replacing 
and  replaced  components  are  merged.  If  the  replacing  and  replaced 
components  belong  to  different  blobs,  the  two  blobs  will  be  merged  and 
the  two  rings  of  perimeter  sections  will  be  merged  into  one  ring.  An 
example  is  shown  below. 


Image  so  far: 


New  segment: 
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Perimeter  lists 
before  new  seg¬ 
ment  is  added: 


Perimeter  lists 
after  new  seg¬ 
ment  is  added : 


2]  1  W  2 


a  re 


•M— t 


To  save  space  during  processing  of  the  image,  and  to  store  the 
completed  perimeter  lists  in  less  disk  space,  the  perimeter  lists  are 
periodically  packed.  The  variable  LINEGROUP  controls  how  often  the 
lists  are  compressed  (the  words  "pack"  and  "compress"  mean  the  same 
thing).  For  example,  when  LINEGROUP  =  25,  compression  occurs  at  lines 
25,  50,  75,  etc.  Those  points  in  the  list  that  previously  have  been 
compressed  are  skipped.  When  a  perimeter  list  has  been  completed  and 
is  ready  to  be  written  to  disk,  it  is  unpacked  and  then  repacked  if 
the  Boolean  variable  FINALPACKING  is  true.  If  FINALPACKING  is  false, 
the  perimeter  list  is  only  unpacked.  Currently,  the  value  of  FINALPACKING 
is  specified  by  the  user  at  the  beginning  of  LINEBYLXNE. 


Each  unpacked  perimeter  point  consists  of  two  16-bit  integers  with 
values  between  1  and  512.  One  Integer  is  for  the  line  number  Y  and  the 


second  integer  is  for  the  column  number  X.  Compression  is  achieved  by 
converting  the  perimeter  information  from  a  series  of  (Y,X)  points  to  a 
starting  (Y,X)  point  and  a  series  of  directions  from  one  point  to  the 
next.  The  starting  point  has  the  same  format  as  an  unpacked  perimeter 
point,  except  that  512  is  added  to  the  column  number.  Therefore,  a 
column  number  greater  than  512  identifies  that  point  as  the  beginning  of 
a  series  of  compressed  points.  Because  the  perimeter  is  8-connected,  there 
are  eight  possible  directions  indicating  the  positions  of  subsequent 
perimeter  points.  They  are  numbered  from  0  to  7  and  are  assigned  as 
follows : 


Change  in  X 
-10  1 

-13  2  1 

Change  in  Y  0  A  0 
15  6  7 

For  example,  if  the  perimeter  list  of  the  last  example  above  (deletion 
of  a  segment,  Subcase  B)  were  to  be  compressed  starting  with  the  point 
(3,1),  the  direction  to  the  next  point  (2,1)  would  be  2;  the  direction 
from  (2,1)  to  the  next  point  (1,1)  would  be  again  2;  and  the  direction 

from  (1,1)  to  the  next  point  (1,2)  would  be  0. 

Each  of  the  eight  directions  can  be  represented  as  a  3-bit  code 
(000  ■  0,  001  *  1,  010  =2,  ...,  Ill  *  7).  Five  direction  codes  fill 
15  bits  of  a  16-bit  word.  The  16th  bit  (the  sign  bit)  is  set  equal  to  1 
to  indicate  that  the  word  contains  compressed  perimeter  points.  Each 
set  of  five  direction  codes  is  stored  alternately  in  the  line  and  column 
locations  of  the  perimeter  point  record.  If  we  continue  the  example, 
the  first  five  directions  would  be  2,  2,  0,  6,  and  7.  The  would  be 
coded  and  placed  in  a  16-bit  word  as: 

1  111  110  000  010  010 

16  15  1A  13  12  11  10  9  8  7  6  5  A  3  2  1 


B-27 


This  word  is  interpreted  as  the  two's  complement  integer  -1756  octal, 
or  -1006  decimal.  The  second  five  directions  are  0,  1,  1,  0,  and  5. 
Their  bit  representation  is: 

1  101  000  001  001  000 

16  15  14  13  12  11  10  987  654  321 


This  is  interpreted  as  -27670  octal,  or  -12216  decimal.  Hence,  the 
first  11  perimeter  points  are  compressed  into  the  space  of  two  perimeter 
point  records  as  follows:  (3,513),  (-1006,-12216). 

Compression  of  groups  of  less  than  10  points  is  not  attempted. 
Therefore,  during  processing,  the  list  of  compressed  points  is  inter¬ 
spersed  with  uncompressed  points.  In  addition,  there  are  usually  a 
few  uncompressed  points  at  the  end  of  the  list.  In  our  example,  the 
12th  point,  (3,8)  remains  uncompressed.  The  completely  compressed 
list  is  therefore: 


The  unpacking  process  is  just  the  reverse  of  the  packing  process. 
The  column  number  of  each  perimeter  point  is  examined.  If  it  is  greater 
than  512,  the  start  of  a  series  of  compressed  points  has  been  detected. 
Subtracting  512  from  the  column  number  results  in  the  (Y,X)  coordinate 
of  the  starting  point.  The  sign  of  the  column  number  of  the  next 
perimeter  point  is  examined.  If  it  is  less  than  or  equal  to  0,  then 
compression  of  10  points  has  been  detected.  The  locations  of  these 
points  are  recomputed  from  the  3-bit  direction  codes  that  are  extracted 
from  the  two  16-bit  words.  This  process  continues  until  a  perimeter 
point  record  with  a  column  number  greater  than  0  is  encountered. 


VI  HIERARCHY  OF  PROCEDURES 


This  chart  shows  the  hierarchy  of  procedures  within  LINEBYLINE. 
All  of  the  procedures  that  are  called  by  each  procedure  are  indented 
and  listed  below. 


LINEBYLINE 

PACKPERIMS 

PKPERIMS 

COMPRESSED 

UNPACKPERIMS 

UNPACK5 

DIRPOINT 

DELETEPERIMS 

NEWPERIM 

DIRECTION 

GETLINE 

ENDOFBLOCK 

SETLINE 

CREATERUNLENGTHS 

ADDRUNLENGTHS 

PROCESSLINE 

DELETESEGMENT 

ADDLEFTPERIMPOINT 

NEWPERIM 

ADDRIGHTPERIMPOINT 

NEWPERIM 

RECORDBLOB 

WRITEPERIMS 

WRITEPOINTS 

PKPERIMS 

(see  above) 
COUNT AND WRITE 
COMPRESSED 
UNPACKPERIMS 

(see  above) 
WRITEPOINTS 

(see  above) 
DELETEPERIMS 

WRITEBLOB 

DELETEBLOB 

INSERTSEGMENT 

NEWBLOB 

ADDRIGHTPERIMPOINT 
(see  above) 

CASE3PR0CESSING 

UPDATEPERIMS 

ADDLEFTPERIMPOINT 
(see  above) 
ADDRIGHTPERIMPOINT 
(see  above) 
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VII  INPUT  IMAGE  FORMAT 


The  binary  image  that  is  read  by  LINEBYLINE  is  512  *  512  pixels. 
There  are  8  pixels  per  byte  (each  pixel  is  one  bit).  Therefore,  one 
line  of  data  occupies  64  bytes,  and  the  entire  image  occupies  512  x  64  = 
32K  bytes  on  disk.  The  data  is  read  one  byte  (one  character  in  PASCAL) 
at  a  time. 

However,  there  is  one  small  difference  due  to  the  PDP  11/40  version 
of  PASCAL.  The  bytes  are  stored  on  disk  in  blocks  of  512  bytes.  The 
PASCAL  I/O  processor  interprets  a  block  boundary  as  an  extra  byte. 

To  compensate  for  this  inaccuracy,  the  LINEBYLINE  software  skips  an 
"imaginary"  byte  once  every  512  times  (when  the  Boolean  function 
ENDOFBLOCK  is  true).  On  most  computer  systems,  this  would  not  be  neces- 
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ANNOTATIONS  FOR  THE  SPSS  RESULT  PACKAGE 
(Submitted  Separately) 


The  printout  of  an  SPSS  analysis  contains  several  parts.  Following 
is  a  description  of  each  part  with  the  page  number  of  the  listing  on 
which  that  part  of  the  printout  begins.  Preceding  the  actual  output, 
some  initial  information  is  listed,  including  the  control  file  used 
to  run  SPSS  and  the  data  that  are  to  be  read  by  SPSS.  These  data  are 
ordered  as  a  sequence  of  trials.  The  numbers  represent:  Subject  code 
(1-30),  session  number,  trial  number,  and  29  feature  values. 

•  Page  1:  SPSS  description;  format  of  the  data  that  are  to  be 
read  in. 

•  Page  3:  A  partial  list  of  the  data  that  were  read  in.  Each 
trial  is  listed,  but  only  five  selected  feature  values  are 
listed. 

•  Page  9:  SPSS  options  used. 

•  Page  10:  Number  of  trials  per  subject  (group). 

•  Page  10 :  Feature-value  means  for  each  subject. 

•  Page  12:  Feature-value  standard  deviations  for  each  subject. 

•  Page  15 :  Wilk's  lambda  and  univariate  F-ratios  for  each  feature. 

Each  of  these  numbers  is  related  to  the  discriminating  power  of 
the  specific  feature.  For  Wilk's  lambda,  the  smaller  the 
number,  the  better  is  the  discriminating  power  of  that  feature. 

The  opposite  is  true  for  F-ratios:  The  larger  the  F-ratio,  the 
greater  is  the  discriminating  power  of  that  feature. 

•  Page  17:  Within-groups  correlation  matrix.  This  is  similar  to 
the  covariance  matrix,  except  that  each  feature  mean  is  first 
normalized  to  mean  *  0  and  standard  deviation  =  1  before  calcu¬ 
lation  of  the  covariance.  A  number  in  this  matrix  that  approaches 
1  indicates  a  high  degree  of  correlation  between  those  two 
features.. 

•  Page  19:  Prior  probability  of  finding  any  particular  subject 
within  the  population.  A  priori,  this  probability  was  set  to 
be  equal  for  each  subject. 
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Page  1^:  'iwenty-nine  discriminant  fuac  taOiu  wc.i ..  determined. 

The  data  in  this  chart  demonstrate  the  declining  discriminating 
power  of  each  additional  discriminant  function. 

Page  20 :  The  standard  discriminant-function  coefficients  are 
listed.  The  product  of  all  normalized  feature  values  and  the 
appropriate  discriminant-function  coefficient  will  be  the  dis¬ 
criminant  score  used  to  classify  each  subject. 

Page  23:  Here  start  listings  of  both  the  unstandardized 
discriminant-function  coefficients  and  the  location  in  29- 
dimensional  space  of  the  centroids  of  each  group.  The  centroid 
is  the  mean  discriminant-function  score  for  each  subject. 

Page  30:  Classification  results:  Actual  group  is  the  subject 
code  number;  "highest  group  D**2"  is  the  classification  result. 

In  all  cases  the  group  was  classified  correctly.  "Probability 
P(G/x)"  is  the  probability  that  the  sample  belongs  to  Group  G 
given  that  the  discriminant  measure  x  was  observed.  "Probability 
P(x/G)"  is  the  probability  that  the  discriminant  measure  x  would 
be  observed  given  that  the  sample  is  from  Group  G.  The  actual 
discriminant  scores  (upon  which  the  classification  is  based)  are 
listed  on  the  far  right. 

Page  78:  Territorial  map  of  first  two  (most  discriminating) 
discriminant  functions.  This  is  a  2D  representation  oi:  rhu  data. 
The  centroids  (*)  are  separated  by  numbers  or  letters  (correspond¬ 
ing  to  that  particular  group,  numbered  0-9  or  A-Z)  that  define 
the  equal  probability  lines  between  groups. 


MISSION 
of 

Rome  Air  Development  Center 

RAVC  plant)  and  executes  research,  development,  tut  and 
selected  acquisition  programs  In  Support  of  Command,  Control 
Communications  and  Intelligence  (C3 I)  activities.  Technical 
and  engineering  support  within  areas  of  technical  competence 
is  provided  to  E SV  Program  Offices  ( POs )  and  other  ESP 
elements.  The  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  and  control,  sur¬ 
veillance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 
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